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Abstract. This paper is a continuation of the study on the stabiUty speed for 
Markov processes. It extends the previous study of the ergodic convergence speed 
to the non-ergodic one, in which the processes are even allowed to be explosive 
or to have general killings. At the beginning stage, this paper is concentrated on 
the birth-death processes. According to the classification of the boundaries, there 
are four cases plus one having general killings. In each case, some dual variational 
formulas for the convergence rate are presented, from which the criterion for the 
positivity of the rate and an approximating procedure of estimating the rate are 
deduced. As the first step of the approximation, the ratio of the resulting bounds 
is usually no more than 2. The criteria as well as basic estimates for more general 
types of stability are also presented. Even though the paper contributes mainly to 
the non-ergodic case, there are some improvements in the ergodic one. To illustrate 
the power of the results, a large number of examples are included. 
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1. Introduction 

Consider a birth-death process on the nonnegativc integers Z_|_ with birth rates 
6„ > (n ^ 0) and death rates a„ > (n ^ 1). Define 

6o • • • bn-i . . 

//q = 1, I2n = , n ^ 1. (1.1) 

oi • • • a„ 

We say that the birth-death process is nonexplosive if the following Dobrushin's 
uniqueness criterion holds: 



fc I II 

" " ' ' oo (1.2) 



oc ^ k - oo oo ^ 

fc=0 1=0 i=u k=i 

(cf. Dobrushin (1952), or Wang and Yang (1992, Corollary 5.2.1), or [10; Corol- 
lary 3.18]). This implies a useful condition that 

When 'Yl^ l^k < oo, each of (1.2) and (1.3) is equivalent to the recurrent condition: 
X^o°(^nMn)~^ = oo. Otherwise, (1.3) cannot imply (1.2) since one can easily 
construct a counterexample so that J2o^ i^k = oo but 

cx) oo J 

1=0 k=i 

Thus, under (1.3), the process may not be unique. 

It is well known that for a birth-death process, the transition probabilities 
{Pij(t)) satisfy 

limp,,(t)=:7r,-^0 (1.4) 

t— >oo 

for all i,j G We are now interested in the exponential convergence rate 
a* = sup |a : \pij{t) — 7rj| = C>(exp[— at]) as t ^ oo for all i,j G E^. (1.5) 

In the ergodic case (i.e., limt^oo Pi jit) > for all i, j), we have Z := Y^Y=oi^3 ^ 
oo and then iXj := tXj/Z > for all j ^ 0. In this case, the problem has been 
well studied, see, for instance, van Doorn (1981; 2002), Zeifman (1991), Kijima 
(1997), [2, 12], and the references therein. The problem becomes trivial in the 
zero-recurrent case for general irreducible Markov chains, since we have on the 
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one hand tTj = for all j, and on the other hand, /q°° pii{t)dt = oo for all i. Hence, 
the exponential decay can only happen in the transient case: 



CXD ^ 



n=0 ^"'^'^ 



< OO. (1.6) 



Since the process is ^u-symmetric: ^iPij{t) = fijPji{t) for all i, j and t, it is natural, 
as we did in the ergodic case, to use the L^-theory. As usual, denote by || • || and 
(•, •) the norm and the inner product on the real Hilbert space L'^ifi), respectively. 
Let 

J(f = {f : f has finite support}. (1.7) 

Define 

with the minimal domain ^™'"(Z)) consisting of the functions in the closure of 
=jr with respect to the norm || • \\d: \\f\\n = + D{f). Next, define 

Ao = inf{D(/) : 11/11 = 1, / e ^} = M{D{f) : \\f\\ = 1, / G ^^'^D)}. 

From now on, we often write /oo or /(oo) as the limit of / at infinity provided it 
exists. In the definition of Aq, it is natural to add the boundary condition /oo = 
but this can be ignored since on the one hand, for each / € ,J^, we have foo = 0, 
and on the other hand Jif is a core of the Dirichlet form (D, iF™'"(L')) (i.e., the 
form is regular) by [10; Proposition 6.59]. For a large part of the paper, we are 
dealing with this minimal Dirichlet form or the minimal process. 

We now make a connection between a* and Aq. The proofs of the next three 
propositions are delayed for a moment. 

Proposition 1.1. For a general non-ergodic symmetric semigroup {Pt}t^o with 
Dirichlet form {D,^(D)) (not necessarily regular) on L'^i/J-), the parameter Xq, 

Xo = M{D{f):\\f\\ = l, fe&{D)}, (1.8) 

is the largest e such that 

\\Ptf\\^\\f\\e-'\ i^0,/GL2(^). (1.9) 



It was proved in [2; Theorem 5.3] that for birth-death processes, under (1.2), 
the exponentially ergodic convergence rate a* coincides with the L^-exponential 
one, denoted by Ai: 

\\Ptf - vr(/)|| ^ 11/ - vr(/)||e-^^* for alH ^ and / G L\p), 



where '/r(/) = / fdp/p{E). For non-ergodic birth-death processes, we have simi- 
larly a* = Ao, as mentioned at the end of [2]. Here is a generalization. 
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Proposition 1.2. For a general non-ergodic ji-symmetric Markov chain with 
Dirichlet form {D, ^(D)), we have a* = Aq defined by (1-8). 

About (1.3), we have the following result. 

Proposition 1.3. Let ^""^(D) = {/ G L'^in) : D{f) < oo}. Then the Dirichlet 

form {D,S^™^^{D)) is regular iff (1.3) holds. In other words, the Dirichlet form 
corresponding to the rates (a^) and (bi) is unique i^f (1.3) holds. 

Proposition 1.2 reduces the study on a* to the first (or principal) eigenvalue 
Aq. This is the starting point of this paper. In the two cases we have discussed so 
far, the state is a reflecting (Neumann) boundary, denoted by code "N". For Ai, 
since the process starting from any point will certainly come back, the infinity may 
be regarded as a reflecting (Neumann) boundary. However, for Aq, the situation 
is different. As we will prove in the next section, the corresponding eigenfunction 
decreases to zero at infinity. Hence, the infinity may be regarded as an absorbing 
(Dirichlet) boundary, denoted by code "D". Thus, for the temporary convenience, 
we rewrite Ai = A^^ and Aq = A^°. Replacing the Neumann boundary at by 
the Dirichlet one (i.e., bo = 0), we obtain two more cases for which we have the 
decay rates (eigenvalues) A^^ and A^^, respectively. The main body of this paper 
is devoted to study these four cases. Now, the rate a* coincides with, case by 
case, one of X™ , A^°, A°^, and A^°. Here are simple examples to show the 
difference in the different cases. 

Examples 1.4. 

(1) Let Qi = Si, bi = pi + j, S > p. Then A^^ = S - ^ if -y > and so does 

A^N if^ = o. 

(2) Let ai = i, h = 2{i + 7). Then X^^ = j if j > and A°° = 1 i/7 = 0. 

The rate in the first example is the difference of the coefficients of leading 
terms, independent of 7. This is somehow natural. Surprisingly, the second one 
is determined by the constant term only except 7 = at which case there is a 
jump from A^° to A^°. Thus, for the convergence rate, the role played by the 
parameters (aj,6j) is mazed and then one may wonder how far we can go (see 
Theorem 1.5 below for a preliminary answer). 

The main body of the paper is devoted to the quantitive study of the con- 
vergence rate. For this, our key result (variational formulas) plays a full power. 
For those readers who arc interested only in the qualitative criteria and basic 
estimates, here is a short statement. 

Theorem 1.5 (Criterion and basic estimates). Let (1.3) hold. Then in spite 
of bo > or bo = 0, the exponential convergence rate a* defined in (1.5) for the 
unique process is positive 

(1) iff S^^-*^ < 00 in the case ofY^-fii < 00; and otherwise, 



(2) 



iffS^^-^^ < 00 



where 




SPEED OF STABILITY FOR BIRTH-DEATH PROCESSES 



5 



More precisely, we have the basic estimate S ^ /A ^ a* ^ S ^ , where the constant 
S is equal to k^^-^^) or k^^-^) according to bo > or bo = 0, respectively: 



- 1 /CO 



(^(^•^))-^= inf [fE^V'+fE^T 



771 — 1 ^ 



j=n •' ■' 



1=1 



-1-1 / m X —1 

E^i 

^ j=n 



Here, the superscript of k^'^'^\ for instance, means that it is in the case studied 
in Section 7 and the constant is given in (7.5). 

The proof of Theorem 1.5 and its extension are given in Section 7. The more 
general quahtative results are presented in Section 8 and in Summary 9.12 for the 
killing case. 

To have an impression about the progress made in the paper, let us have a look 
at the new points made in the well-developed case, Section 6. 

(1) The uniqueness condition (1.2) is replaced by using the maximal process. 

(2) The more complete dual variational formulas are presented in Theorem 6.1. 

(3) Even though the criterion Theorem 6.2 is known before, the upper bound 
in its improvement (Corollary 6.4) is newly added so that the ratio of the 
bounds is now no more than 2 as shown by a group of examples. Moreover, 
a new criterion (Corollary 6.6) which has been expected naturally (in view 
of Theorem 4.2) for a long time, is now presented. 

(4) A more effective sequence for the upper estimate given in Theorem 6.3 is 
introduced to replace the original one. The monotonicity of the approxi- 
mating sequences arc proved here for the first time. 

We now return to prove the propositions above. 

Proof of Proposition 1.1. Replace by emax the largest exponential rate in (1.9). 
Then we have Emax ^ because of the contractivity of the semigroup in every 
L^-space (p ^ 1). We need to show that Aq = emax- The proof of Aq ^ Emax is 
easier since by an elementary property of the Dirichlet form and (1.9), we have 
for every / with ||/|| = 1, 

D{f) = lim t ^(/ - Pt/, /) ^ lim ^(1 - e"^— *) = w, (1.10) 

where lim ^ means an increasing limit. Hence, we have Aq ^ ^max- To prove 
£max ^ Ao, assume that Aq > 0. Otherwise, the assertion is trivial. Noticing that 
D{f) = (-0/, /) for the generator Q of {Pt} on L'^ifJ,) and for every / G ^{^), 
we have 

^jPtff = 2{Ptf, mf) = -2D{Ptf). (1.11) 
Next, since Ptf G ^{D) for each / G L'^in), by the definition of Aq, we have 

~2D{Ptf) < -2Xo\\Ptff. 
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Thus, IIPJII ^ ll/lle"^"* for alH ^ and / G and then for all / G L'^{n) 

since the density of in L'^{n) and the contractivity of the semigroup {Pt}t^o- 
The assertion now follows since ffmax is the largest rate. □ 

Proof of Proposition 1.2. The proof for a* ^ Aq is rather easy. Simply applying 
Proposition 1.1 to the indicator function / = we obtain 

Note that this also provides a non-trivial estimate of the constant in (1.5). 

To prove that Aq ^ ct* , we may assume that a* > 0. One may follow the proof 
of [12; proof of part (4) of Theorem 8.13]. In the last part of the original proof, 
we have 

\\Ptff = if,Pnf) ^ ll/llLe-^"** E 

for every bounded / with compact support. Here, we have used the assumption 
that ^ Cij-e""**. □ 

Proof Proposition 1.3. Since the Q-matrix is conservative, by [10; Lemma 6.52 
and Theorem 6.61], (D , ^"^'^^ (D)) is a Dirichlet form and is indeed the maxi- 
mal one. Note that in the conservative case, every Q-process (in particular, the 
semigroup generated by a Dirichlet form) satisfies the backward Kolmogorov's 
equation by [10; Theorem 1.15(1)]. 

(a) Let (1.3) hold. Then the Dirichlet form should be regular. Otherwise, we 
have two different birth-death semigroups generated by {D,^^^^{D)) and the 
minimal Dirichlet form (D , !^™^'^ (D)) , respectively. They satisfy first the back- 
ward and then also the forward Kolmogorov's equations by [10; Theorem 6.16]. 
This is impossible since condition (1.3) is the uniqueness criterion for the pro- 
cess satisfying the Kolmogorov's equations simultaneously, due to Karlin and 
McGregor (1957a, Theorem 15) (cf. Hou et al. (2000, Theorem 6.4.6 (1); 1994, 
Theorem 12.7.1)). Note that criterion (1.3) is equivalent to the uniqueness for the 
process satisfying one of the Kolmogorov equations since every symmetric process 
as well as the minimal one satisfies both of the equations. This is the reason why 
(1.3) is weaker than (1.2). 

(b) Next, let (1.3) fail. Then we have < °° < 
Moreover, (1.2) fails. Note that the birth-death Q-matrix has at most a single 
exit boundary, and there is precisely one if (1.2) fails. Besides, the non-trivial 
(maximal) exit solution zx is bounded from above by 1. In view of [10; Proposi- 
tion 6.56], there are infinitely many Dirichlet forms. The minimal one is regular 
but not the maximal one {D, Sl™^{D)). □ 

Actually, Proposition 1.3 is a particular case of a result we will study at the 
end of Section 9 (Theorem 9.22). 

The remainder of the paper is organized as follows. In the next two sections, 
we study A'^^. Sections 4, 6 and 7 are devoted to A^^, and A^^, respectively. 
By exchanging N and D, we formally obtain a dual of A'^^ and A^^ (resp. A'^'^ 
and A'^^) which is studied in Section 5 (resp. 7). In each case, we present a group 
of dual variational formulas for the first (non-trivial) eigenvalue. By using the 
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formulas, we then deduce a criterion for the positivity of the eigenvalue and an 
approximating procedure for estimating the eigenvalue. The criteria and basic 
estimates in a quite general setup are given in Section 8. A closely related topic, 
having general killings, is studied in Section 9. In the study of this paper, the 
author has benefited a great deal from our prcvioTis work and from many authors' 
contribution. A part of the contributions is noted in the context. In the ergodic 
case under (1.2), a large number of references are given in [10, 12] and the author 
apologizes for omitting them here. At the end of the paper (Section 10), some 
remarks on the related results, some open problems or open topics, and so on 
are discussed. The analog of Theorem 1.5 for one-dimensional diffusions is also 
included. 

Notation 1.6. To be economical, we use the same notation Xq, S, k, I and 
II and so on, from time to time in different sections with similar hut different 
meaning. To distinguish them if necessary, we write Ag for instance to denote 
the Xq defined hy formula (#). 

2. Absorbing (Dirichlet) boundary at 
infinity: dual variational formulas 

This section begins with the study on the property of eigenfunction of Aq- It 
is fundamental in our analysis and has been studied several times before, see, for 
instance, [3; Lemma 4.2]; [4; proofs of Theorems 3.2 and 3.4]; Chen, Zhang and 
Zhao (2003, Section 2); Shao and Mao (2007, Proposition 3.1). The main body 
of this section is devoted to prove a group of variational formulas (Theorem 2.4 
and Proposition 2.5). Their applications are given in the next section. 

Fix an integer N: 1 < AT ^ oo, and let = {fc G Z+ : ^ A; < AT + 1}. 
Throughout the paper, the infinite case that N = co \s more essential but the 
finite case that A < oo is also included which may be meaningful in matrices 
theory. To avoid the confusion of these two cases in reading the paper, one may 
read the infinite case first and then go back to check the modification for the finite 
case. Besides, note that when N < oo, neither (1.2) nor (1.3) is needed. 

Let us start at a general situation. Consider the operator Q corresponding to 
the birth-death Q-matrix with birth rates bi, death rates Oj, and killing rates 
Ci'^ 0{i e E) as follows. 

nf{i) = bi{fi+,-fi)+ai{fi.^- fi)-Cifi, i e E, = if AT < oo. (2.1) 

In other words, when A < oo, the state A^+1 is an absorbing (Dirichlet) boundary. 
When Ci ^ for 1 < i < A", unless otherwise stated, we assume that Oq = and 
^AT = if A < oo (the other and bi are positive), otherwise, simply replace the 
original cq and cjv by oq + cq and 6jv + cjv, respectively. Now, since oq = 0, /_i 
is free in the last formula. The first eigenvalue Aq is now defined by 

Ao =inf{Z)(/): 11/11 = 1,/G^}, (2.2) 

where 

^(/)= E Mi«'i(/m-/if + /jv+i = if AT < oo. (2.3) 

0^i<N ieE 
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We say that g is an "eigenfunction" of A G M, if ^ satisfies the "eigenequation" : 

^g = -Xg, gN+1 = if < oo. (2.4) 

Note that the "eigenvalue" and "eigenfunction" used in this paper are in a gener- 
ahzed sense rather than the standard ones since here we do not require g E L'^ (/x) . 

Proposition 2.1. 

(1) Every eigenfunction g o/ A € M satisfies 

k 

IJ-khigk - 9k+i) = - Ci)nigi, k e E, gN+i = if N < oo. (2.5) 

(2) // Ao > 0, then q ^ (0 ^ i ^ N) whenever N < oo, and the non-zero 
eigenfunction g of Aq is either positive or negative on E. 

(3) The non-zero eigenfunction g ofX = Ois either positive and nondecreasing, 
or negative and nonincreasing on E. Furthermore, let g > for instance. 
Then gk+i > gk for all k : i ^ k < N whenever Cj > 0. 

Proof, (a) Recall the eigenequation 

ng{i) = bi{gi+i - gi) + ai{gi-i - gi) - Cigi = -Xgi, ieE, (2.6) 
or more generally, the Poisson equation 

bi{gi - gi+i) - ai{gi-i - gi) = hi, ie E, gN+i = if iV < oo, (2.7) 
for a given function h. Multiplying both sides by //j, we get 

i^Aigi - 9i+i) - i^i-ibi-i{gi-i - gi) = i^iK, ieE. (2.8) 

When i = 0, the second term on the left-hand side is set to be zero. Making a 
summation over i, we obtain 

k 

IJ'khigk - 9k+i) = XI '"'^i' ke E, gN+1 = if iV < oo. (2.9) 

i=0 

With hi = (X — Ci)gi, this gives us the first assertion of the proposition. 

(b) To prove the second assertion, note that Aq = if q = (0 ^ i ^ < oo) 
in which case both and N are reflecting and the process is ergodic. Now, since 
Ao > 0, one may assume that 7^ 0, otherwise = by induction. Next, 
replacing g hy g/g^ if necessary, we can assume that go = 1. If 5 is not positive, 
then there would exist a ko £ E, ko 1 such that gi > for i < ko and gk„ ^ 0. 
We then modify g from ko: set gi = gi for i < ko and gi = for i ^ ko + 1. 
By choosing a suitable value £ > at ko, the new function g G gives us 
< which is a contradiction to the definition of Aq. Hence, g does 
not change its sign. 
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We are now going to specify e. Note that 

{-ng){ko - 1) = -bko-i{e - gk„-i) + ak„-iigk„-i - gko-2) + Ck^-igko-i 
= i-^9)iko - 1) + bk„-iigko - £) 

= Aofffco-l + bko-l{9ko — ^) 
< '^oS'fco-l 

since £ > ^ gko- Note also that 

{-ng){ko) = -6fco(0 + afeo(£ - 9ko-i) + Cko^ = ^(afco + ho + Cko) - «feoffco-i- 

Next, since D{f) = (/, —0,f) for every / G and for each i, i^f{i) depends on 
three points i and i±l only, we obtain 

^(9) = X] f^i9i{-^9){'i) + IJ'ko-i9ko-i{-^9){ko - 1) + /^fco5fco(-^5)(^o) 

feo-1 

< Ao ^ IJ-igi + enkoHako + ho + Cfeo) - afeoffco-i]- 

i=0 

Because 

fco-l 

\\9f = Yl f^i9i + l^ko^'^ ^ 
for < Ao, it suffices that 

£[e(«feo + ho + Cfco) - ako9ko-i] < Aoe^- 

Equivalently, 

e(afco + ^fco + Cfcn - -^o) < ako9ko-i- 

This clearly holds for sufficiently small £ > 0. 
(c) If A = 0, then (2.5) becomes 

k 

l^kh{9k+i - 9k) = CiHigi, ke E, g^+i = if iV < 00. (2.10) 

i=0 

Clearly, if go = 0, then = by induction. Without loss of generality, assume 
that (7o = 1- By (2.10) and induction, it follows that gk+i — gk ^ ^ for all i E E. 
Actually, gk+i > 9k for all /c: i ^ k < N provided Cj > 0. □ 

In view of (2.5), the eigenfunction g may not be monotone if q ^ 0. 

For the remainder of this section, we assume that Cj = for i < N but cn > 

if N < 00. However, to simplify our notation, set Cj = but letb^ > if N < 00. 

In view of the definition of the state space E, the point N + 1 is regarded as a 

(2 2) 

Dirichlet boundary. From now on in the paper, when we talk about Xq , it is 
defined by (2.2) but in the present setting. 
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Proposition 2.2. 

(1) Let g be a non-zero eigenfunction of Aq > 0. Then g is either positive or 
negative. 

(2) Let g be a positive eigenfunction of X > 0. Then g is strictly decreasing. 
Furthermore, 

N ^ k 
k=n 1=0 

In particular, 

N N n 

^IJ'n9n'^[n,N] = ^Un^l^kgk = — — < oo, (2.12) 
n=0 71=0 k=0 

where ^'[£,m] = YlT-e'^'^' = {i^kbk)~^- Moreover, if (1.2) holds, then 
ffoo := liniAr^ooffAf = 0. 

(3) Let Ao = 0. Then N = oo and the eigenfunction g must be a constant 
function. 

Proof, (a) The first assertion follows from Proposition 2.1 (2). 

(b) Let A > 0. Since ^ > 0, by (2.5) with Cj = 0, it follows that gi is strictly 
decreasing in i. By (2.5) again, we have 

N N ^ k N 

9n - gN+1 = ^{9k - 9k+i) = — J-'^l^idi = A^MiSfji^i^ Vra,A^]. 

k=n k=n ^ i=0 i=0 

We obtain formula (2.11) and then (2.12). If goo > 0, then by condition (1.2), the 
left-hand side of (2.11) is bounded below by 

oo ^ k 

k=n t=Q 

which is a contradiction since the right-hand side of (2.11) is bounded from the 
above by 50/A < oo. Therefore, we must have goo = 0. 

With some additional work, condition (1.2) for = will be removed (see 
Proposition 2.5 below). 

(c) We now prove the last assertion of the proposition. When < 00, it is 
well known that Aq > 0. Now, let Aq = and then N = 00. By (2.6) with Cj = 0, 
we have ^ 

gi+i - 9i = -r-{9i ~ 9i-i), i > 0. 

Oi 

From this and induction, it follows that gn = go for all n ^ 1 since ao = 0. □ 

We remark that for finite state space with absorbing at A^+1 < 00, Proposition 
2.2 was actually proved in [4; proof d) of Theorem 3.4] with a change of the order 
of the state space. Next, when N = 00 and Aq > 0, in contrast with the ergodic 
case where g G L^(fi) (cf. [12; Proposition 3.5]), here one may have g ^ L'^{ji) 
and then g ^ L^{pi). However, g G L^{v) since gn is strictly decreasing and 
f„ < 00, which is a consequence of Theorem 3.1 below. 
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Corollary 2.3. Let Aq > 0. Then lirrij^oo Pt/(i) = for allt^O and f eJ^. 

Proof. It suffices to show that hrrij^oo Yl^=iPik{'t) = 0- We now prove a stronger 
conclusion: linii^oo Ptg{i) = for all t ^ 0, where (/ > with (/q = 1 is the 
eigenfunction of Aq. Since g is bounded, by using the well-known fact that 

e-^°*5i = Ptg{i), t ^ 0, 

the conclusion now follows from Propositions 2.2 and 2.5 (2) below. □ 

For a specialist who does not want to know many details, at the first reading, 
one may have a glance at the remainder of this section and the next section, 
especially Proposition 2.7, and then go to Section 4 directly. From here to the 
end of the next section, we are dealing with a case which is a dual of the one 
studied in Section 4. However, for the reader who is unfamiliar with this topic, it 
is better just to follow the context since we present everything in detail in these 
two sections. A large part of the details in Sections 4 and 6 are omitted since 
they are supposed to be known. 

To state the main results of this section, we need some notation. First, we 
define two operators as follows. 

They are called an operator of single sum (integral) or double sum. respectively. 
Here for the first operator, we use a convention: /jv+i = if AT < oo. The second 
operator can be alternatively expressed as 

1 m 

IIi{f) = YJ2t"^f''''l'^^^^^^ i^[^,rn]=J2^i, i^i = —r. (2.15) 
•'^ fees i=l * 

Next, define a difference operator R as follows. 

Ri{v) = ai[l-vl\)+hi{l-Vi), iG£;, 'y_i>0 is free, i;iv:= if iV<oo. (2.16) 

The domain of the operators //, I and R are defined, respectively, as follows. 

^77 = {/:/>0onE}, 

= {f : f > OIL E and is strictly decreasing}, 
'ri={v: for alH(0 ^i<N), ^^€(0,1) if T.ji^j <oo and G (0, 1] ii^ji/j=oo}. 

These sets are used for the lower estimates. For the upper estimates, we need 
some modifications of them as follows. 

^11 = {/:/> up to some m: l^in<N + l and then vanishes}, 

^/ = {/ : / is strictly decreasing on some interval [n, m] {0 ^ n < m < N + 1) , 

fi = fn for i ^ n, /m > 0, and = for f > m), 
■^1 = ai+i{ai+i + bi+i)-^ <v,<l - ai(t;-_\ - l)^-^ 

for i = 0, 1, . . . ,m — 1 and f j = for i ^ m} . 
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Here and in what follows, to use the above operators on these modified sets, 
we adopt the usual convention 1/0 = oo. Besides, the operator // should be 
generalized as follows: 

IIi{f) = j Y. ^E/^fc^' ^esupp(/). (2.17) 

Prom now on, we should remember that II,{f) is defined on supp(/) only. For- 
tunately, we need only to consider the following two cases: either supp (/) = 
{0, 1, . . . , m} for a finite m or supp (/) = E. 

To avoid the heavy notation, we now split our main result of this section into 
a theorem and a proposition below. 

Theorem 2.4. The following variational formulas hold for Aq defined by (2.2). 

(1) Difference form: 

inf sup Ri{v) = Ao = sup inf Ri{v). 

(2) Single summation form: 

inf sup/i(/)"^ = Ao = sup inf /j(/)"\ 
fe^i ieE fe^i 

(3) Double summation form: 

inf sup //j(/)-i = Ao= sup inf//i(/)"^ 
fes^ii iesupp(/) fe^n »eB 

Moreover, the supremum on the right-hand side of the above three formulas can 
be attained. 

The next result extends the domain of Aq or adds some additional sets of test 
functions for the operators I and //, respectively. Roughly speaking, a larger set 
of test functions provides more freedom in practice and a smaller one is helpful 
for producing a better estimate. 

Proposition 2.5. 

(1) We have 

Ao = inf{Z?(/) : ||/|| = 1, /jv+i = 0}, (2.18) 

where /oo := limiv->oo /jv in the case of N = oo. 

(2) When Aq > 0, the eigenfunction g satisfies gN+i = 0. 

(3) Moreover, we have 

Ao= inf sup/i(/)-i (2.19) 

= mf _ sup //i(/)-\ (2.20) 
fes^nvi^'jj jesupp(/) 

inf sup IIi{f)-^ = Ao = sup inf IIi{f)-\ (2.21) 
fe^i iesupp(/) fe^i 'S-^ 
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where 

f is strictly decreasing and positive up to some m:l^m<N-[-l 
and then vanishes^ C ^i, 
f>0 onE and fll{f) G 

Besides, the supremum supf^^^^ in (2.21) can also he attained. 

The condition "/at+i = 0" in (2.18) explains the meaning of "absorbing (Dirich- 
let) boundary at infinity" used in the title of this and the next sections. 

Among the different groups of variational forms, the difference form is the 
simplest one in the practical computations. For instance, when N = oo, by 
choosing t;j = c < 1, we obtain the following simple lower estimate: 

Ao ^ inf [6i(l-c)-ai(c-^-l)]. 

This is non-trivial and is indeed sharp for a linear model (Example 3.5, c = 1/2). 
The difference form of the variational formulas will be used in Section 5 to deduce 
a dual representation of Aq. In general, the estimates produced by the operator R 
can be improved by using the operator / and further improved by using //. The 
price is that more computation is required successively. The single summation 
form of the variational formulas enables us to deduce a criterion for Aq > 
(Theorem 3.1). Whereas the double summation form of the variational formulas 
enables us to deduce an approximating procedure to improve step by step the 
lower and upper estimates of Aq (Theorem 3.2). 

Next, we mention that when N = oo, for the upper estimates (the left-hand 
side of the formulas given in Theorem 2.4 or the formula given in (2.20)), the 
truncating procedure or the condition "///(/) G L^(;u)" cannot be removed. For 
instance, the formally dual formula info<,;^i supjg^i?j(v) of the lower estimate 
supo<i,^i infies Ri{v) {= sup^^-^^ inf^gs Ri{v)\ is not an upper bound of Aq, and 
is indeed trivial. To see this, simply take Vi = l{i< oo). Then Ri{v) = and so 

inf sup Ri{v) ^ sup Ri{v) = 0. 

More concretely, take 6j = 2 and Oj = 1. Then for = c < 1, we have 

inf sup Ri{v) ^ inf sup Ri{v) = 1 inf (1 — c) = 0, 
"eri c<i c<i 

but Ao = (\/2 — l)^ as will be seen in the next section (Example 3.4). Therefore, 
the quantity info<^^i sup^^^; Ri{v), as well as inf^>o sup^g^; Ri{v), has no use for 
an upper estimate of Aq. 

Proofs of Theorem 2.4 and Proposition 2.5. 

Part I. Recall that Aq denotes the one defined by the formula (#). In particular, 

(2 2) 

the notation Aq used from now on in this section is Ag ' . 



^n = {f-- 
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To prove the lower estimates, we adopt the following circle argument: 

Ao ^ aJ'-'') ^ sup inf = sup inf = sup mf L{f )-^ 
^ sup inf Ri{v) ^ Aq. (2.22) 

Clearly, Aq^'^^^ = Aq if iV < oo. However, the identity is not trivial in the case 

of N = oo. Besides, we will show that each suprcmum in (2.22) can be attained; 
and furthermore the eigenfunction g satisfies Qn+i = whenever Aq > 0. 

(a) Prove that Aq ^ A^^'^^^ ^ supjg,^^^ infjgB //i(/)-^ 
When N = oo, the first inequality is trivial since 

{11/11 = l,/e^}c{||/|| = l, /oo = 0}. 

The proof of the second inequality is parallel to the first part of the proof of 
[4; Theorem 2.1]. Let g satisfy gw+i = and \\g\\ = 1, and let (hi) be a positive 
sequence. Then by a good use of the Cauchy-Schwarz inequality, we obtain 

^ = ^IJ'i9l (since ll^ll = 1) 

i 

= XI IZ^^J ~ j (since ^jv+i = 0) 



Exchanging the order of the first two sums on the right-hand side, we get 

N 



j i^j k=i 

N 



^ D{g) supT^^/^i^ 



hk 



3eE hj ^ ^ ^lkhk 
=: D{g) swpHj. 

We mention that the right-hand side may be infinite but we do not care at the 
moment. Now, let / G satisfy c := supj^^ II j{f) < oo and take hj = 
Yli<jl^ifi- Then hj ^ cfj/vj < oo for all j. By the proportional property, we 
have 

I ^ h 1^1 

sup Hj ^ sup -- ^ -j- = sup ^ — r- mfi = sup IIj{f) < oo. 

jeE jeE Jj ^ l^kOk jeE Jj ^ f^kOk ^ jeE 
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Combining these facts together, we obtain Aq^'^^'' ^ infj^o IIj{f)~^ whenever 
supjg£^ IIj{f) < oo. The inequahty is trivial if supj^^^ ^^jif) — ^ so it holds 
for all / G ^jj. By making the supremum with respect to / G ^jj, we obtain the 
required assertion. 

(b) Prove that sup inf = sup inf = sup inf/i(/)~^. 

fe^n f€.^i f(z^j ieE 

Let / G C ^u. Without loss of generality, assume that supjg£;/j(/) < oo. 
By using the proportional property, we obtain 

1 ^ 1 

sup Iliif) = sup — V — — V Hkfk 

W , I N 



^ ^^P XI ~T' yii^kfk Y\{fj - fj+i) (since /jv+i ^ 0) 
ieEj^-Hbjttj I U 

^ sup - — \ — ( ^— V] iikfk] (note that fi > /j+i) 

ieE ji - ji+i xf'iOi rr. J 



= sup/i(/) < oo. (2.23) 

ieE 

Making the infimum with respect to / G , we get 

inf sup JJj(/) ^ inf sup/i(/). 

fe^i ieE fe^i ieE 

Since C the left-hand side is bounded below by infjgj^j.^ supjgg //i(/). 
We have thus proved that 

sup inf//i(/)-i> sup inf7/i(/)-i^ sup iui Ii{f)-\ 
fe^ii fe^i fe^i '^-^ 

There are two ways to prove the inverse inequality. The first one is longer but 
contains a useful technique. Let / G with c := sup^^^ //»(/) < oo. Set 

JV ^ AT 

5i = X m 1^''^'' ^ X X '^'^•^'^ ^ ^' * e gN+1 := if iV < oo. 

i=i k^j j=i k^j 

Then is strictly decreasing in i, gi < go ^ c/o < oo for all i. Hence, g G =^7. 
Noticing that 

N N 

gi - gi+i = ^^3^ f^kfk - X ^-J' X ^''f'' = ^» m ^'^J''^ 

j=i k^j j=i+l k^j k^i 

(here and in what follows, J2k=i nieans J2i^k<j+i ^^'^ = by the standard 
convention), we have 

ng{i) = bi{gi+i - gi) + ai{gi-i - gi) 

k^i k^i—1 

= — y2 ^^kfk H 7 — x ^^^-^^^ 

= -/i, 1 ^ i < N. 
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Actually, this holds also for z = and i = N it N < oo. Applying (2.7) to h = f, 
by (2.9), it follows that 

l^khigk - gk+i) = V i^jgjfj/gj > V i^jgj inf lli{f)~\ keE. 

j^k j^k 

That is. 

Making the supremum with respect to fe, we obtain 

inf 4(5)-i ^ mUIi{f)-\ 
keE leE 

and hence, 

sup inf Ik{g)-^ ^ inf IIi{f)-\ 

This lower bound becomes trivial if sup^^^ IIi{f) = oo, and hence, the inequality 
holds for all / G ^jj. Making the supremum with respect to / G ^n, we obtain 

sup inf Ik{g)~^ ^ sup inf IIi{f)~'^ . 

We have thus proved the required assertion. 
The second proof is to show that 

sup inf Ii{f)~^ ^ Ao 

and thus completes a smaller circle argument. To do so, without loss of generality, 
assume that Aq > 0. Let g > he the eigenfunction of Aq. Applying (2.9) to 
h = Xog, we obtain Ii{g) = Ag ^ for all i E E, and hence, inti^E Ii{g)~^ = Aq. 
Noticing that g G by Proposition 2.2, the assertion is now obvious. 

(c) Prove that supjgj?^^ infi^E IIi{f)~^ ^ sup„g^^ infjee i?j(t;). 
Note that by a change of the sequence {vi}fj^^: 

Ui = vqVi ■ ■ ■ Vi-i, i E E, v-i > is free, vn '■= it N < 00, 

the quantity Ri{v) becomes 

ttj ( 1 — '^h::! I -I- 1 _ !fl±i \ ieE, n_i > is free, ujv+i := if iV < 00. 
V Ui J \ Ui J 

To save our notation, we use Ri{u) to denote this quantity. Clearly, {ui} is 
positive and Vi ^ 1 for all i mean that {ui} is non-increasing. 
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Before moving further, we prove that if inf Ri{u) > for a positive sequence 
u = (uj), then Uj must be strictly decreasing in i. To do so, let 

fi = (oi + bi)ui - aiUi-i - biUi_^i. 

Then /, = UiRi{u) > for all i E E by assumption, and so / G ^jj. Noticing 
that 

IJ-kfk = Mfc+iafc+i('"fc — ""fc+i) — IJ-kakiuk-i — Uk), 

we obtain 

<^IJ'kfk = Hj+iaj+i{uj - Uj+i) = iijhj{uj - Uj+i). 

Hence, Ui is strictly decreasing in i (equivalently, Vi := Ui^i/ui < 1). This proves 
the required assertion. The reason of using rather than {v : Vi > 0,0 ^ i < N} 
should be clear now. 

We now return to our main assertion. For this, without loss of generality, 
assume that inf^g^; Ri{u) > for a given strictly decreasing u = (uj). Otherwise, 
the assertion is trivial. From the last formula, we obtain 

N N 
j=i k^j j=i 

Therefore, 

0<Ri{u) = ^^fi (j2 ^'kfk) = Ih{f)-\ i G E. 

* ^j=i k^j ' 

It follows that 

inf i?i(n) ^ inf ^ sup inf 77i(/)-\ 

The assertion now follows by making the supremum with respect to u. 
(d) Prove that sup^g-^^ infjgB Riiy) ^ Aq. 

Assume that Aq > for a moment (in particular, if AT < oo). Then by Propo- 
sition 2.2, the corresponding cigcnfunction g (with = 1) of Aq is positive and 
strictly decreasing. From the eigenequation 

-^lg{i) = Xogi, i G E, gN+i := if iV < oo, 

it follows that 

ai(l-^)+6.(l-^)=Ao, ieE. 
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Let Vi = gi^i/gi. Then Vi G (0, 1) for all i < N and so v = (vi) G 1^. Moreover, 
Ri{v) = Xq for all ieE. Therefore, we certainly have sup^g^^ infjee Ri{v) ^ Aq, 
as required. 

It remains to prove that sup^^-^^ infjgE Ri{v) ^ when N = oo. First, let 
Ylilo < oo. Choose a positive / such that 

oo oc 

X] l^kfkV^k < OO, ipk := ^ i^j. 

k=0 j=k 

Define 

oo 

hi = '^ Uj ^ fikfk, i ^ 0. 

j=i k^j 

Then 

oo oo 

= ^ fJ'kfk'PiVk ^ X] l^kfk'Pk < OO. 
k=0 fc=0 

Set Ui = hi+i/hi{i ^ 0). Then P G '^i since /ij is strictly decreasing. A 
simple computation shows that Ri{v) = IIi{f)~^ > for all i ^ 0. Hence, 
sup„g^^ infi^o Ri{i^) ^ 0. Next, let Yli'^i — oo and set Vi = 1. Then Ri{v) = 
and so the same conclusion holds. 

The proof of the last paragraph indicates the reason why in fi we define "f j G 
(0,1)" and ''vi G (0,1]" separately according to "X^j z^i < oo" or "J^j = oo". 
Although we have known from proof (c) that for inf^ Ri{v) > 0, it is necessary that 
V < 1 but this condition may not be sufficient for infji2j(i;) ^ 0. The extremal 
Vi = 1 is used only in the case of J2i z^i = oo in which we indeed have Aq = (cf. 
Theorem 3.1 below). 

We have thus completed the proof of circle (2.22). 

(e) We now prove that each supremum in (2.22) can be attained. The case 
that Aq = is easier since 

= Ao ^ inf IIi(/)-i ^ and = Aq ^ inf /i(/)-i ^ 

ieE ieE 

for every / in the corresponding domain, as an application of (2.22). Similarly, 

the conclusion holds for the operator R as seen from proof (d): noting that in 
the degenerated case that Vi = oo, we have Aq = and then = 1 by 
Proposition 2.2 (3). 

Next, we consider the case that Aq > with eigenfunction g: go = 1. Then 
for the operator R, the supremum is attained at Vi = gi^i/gi as seen from the 
first paragraph of proof (d). For the operator /, it is attained at / = as an 
application of Proposition 2.1 with c, = 0: Ii{g) = Aq ^. At the same time, in 
view of part (2) of Proposition 2.2, we have IIi{g) = whenever gN+i = 0. 

It remains to rule out the possibility that gN+i > 0. Otherwise, by part (2) of 
Proposition 2.2 again, we have N = oo and 

Mi ■-^''j^f^k e (0,oo). 

j^i k^j 
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Let g = g — goo- Then g G ^n- Noting that 

^ gj- goo 
Ao 

we obtain 



gooMi (by (2.11)), 



sup Hi (^) = sup 

i>0 i>Q 



1_ _ g^Mj 

Ao gi — 5oc 



1 . , Mi 
goo mf 

Ao 1^0 gi - g, 



oo 



By using the proportional property and (2.5), it follows that 

inf^^^inf ^15^ = 1. 
i>0 gi - g^ i^o gi - gi_^_i Aq 

Thus, we get 

suplli{g) ^ —{I -goo) < T-. 
Hence, inf IIi{g)~^ > Ao, which is a contradiction to proof (a): Aq ^ inf II{g)~^. 

i^O i(zE 

We have thus proved that ^oo = whenever Ao > 0. Note that this paragraph 
uses Proposition 2.2 and proof (a) only. 

Part H. Next, to prove the upper estimates, we adopt the following circle argu- 
ment: 

Ac ^ mf _ sup Iliif)-^ (2.24) 

f€.S^nLI^'n iGsupp (/) 

^ inf sup Iliif)-^ (2.25) 
iesupp (/) 

= inf sup Iliif)-^ = inf supli(/)-^ (2.26) 

fe.^i iesupp(/) fe^i ieE 

^ inf sup/i(/)-i (2.27) 

^ inf sup Ri{v) (2.28) 

veTi ieE 

^ Ao. (2.29) 

Since inequalities (2.25) and (2.27) are obvious, we need only to prove (2.24), 
(2.26), (2.28) and (2.29). 

(f) Prove that Aq ^ mi^^^^^^^,^ sup^esuppC/) ■ 

Wc remark that in the particular case that the eigenfunction / is in L'^^j,), then 
the function g := f 11(f) is nothing but just //Aq G L'^{fi). Hence, the infimum 
in (2.24) is attained at this / G and the equality sign in (2.24) holds. 
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We now consider the general case. Let / G ^n. Then there exists an m such 
that fi>OfoTi^m and /j = for z > m. Set g = Igupp (/)///(/). That is, 

1 0, i^m + 1. 



Clearly, g G and 



0, z ^ m + 1. 



gi - gi+i = 

We now have 

^is) = I^M9i+i-9if = Yi9i-9i+i)Yl^kfk = l^kff' Y i9i~9i+i)- 

i^m i^m k^i k^m k^i^m 

Since gm+i = 0, we get 

-^(f) = y2 I^kfk9k ^ y2 i^ksl max {fi/gi) = \\gf sup 

Dividing both sides by WgW^ G (0, oo), it follows that 

Xo^D{g)/\\gf^ sup IIi{fr\ f^^n- (2.30) 
iesupp (/) 

For / G the same conclusion clearly holds if < oo. When N = oo, since 
g G L'^i/J') by assumption, we have < 51 < oo. As a tail sequence of a convergent 
series (which sum equals go), we certainly have g^ X goo = as i t 00. Hence, the 
same proof replacing m with 00, plus the fact that Aq = Aq^'^^'* proved in Part I, 

shows that 

Ao = Af)^ sup IIi{f)-\ fe^'n- 

iesupp (/) 

Combining this with (2.30), we prove the required assertion. 

The proof indicates the reason why the truncating procedure is used for the 
upper estimates since in general the eigenfunction g may not belong to LF'{lJi) as 
shown by Proposition 2.2. 

(g) Prove that 

inf sup IIi{f)~^= inf sup IIi{f)~^ = inf sup/j(/)~^. 
fe^ii iesupp (/) fe^i iesupp (/) fe^i ieE 



Let / G Then there exist n < m such that = /ivrel{i^m}) fm > 0; and 
/ is strictly decreasing on [n,m]. Clearly, we have 

min//j(/) = min //»(/) and inf 7j(/) = min Ii{f) 

i^m n^i^m ieB n^i^m 
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since, by assumption, 1/0 = oo. By the proportional property, first we have 



min IIi{f) = min ^ ^fikfk/ ^{fj - fj+i) 



and then 



sup inf IIi{f)^ sup inf IIi{f) ^ sup inf 
*esupp(/) «esupp(/) f^^j^^^ 

since C ^77. 

As in proof (b), there are two ways to prove the inverse inequality. First, let 

/ G ^7/- As in proof (f), set g = Isupp (/)/^^(/)- Clearly, g e C &nd 
moreover, 



bi{gi+i - gi) + ai{gi-i - gi) = — ^V'/ifc/fcH ^ l^-kfk 

= ^ Mfc/fe H ^ Affc/fc 

= -fi, i ^ rn. 



When i = 0, the second term on the left-hand side disappears since oq = 0. It 
follows that 

fJ'ibiigi+i - 9i) + IJ^iaiidi-i - 9i) = -fJ'ifi, i ^ rn, 
and furthermore, 

l^khigk - 9k+i) = yZH9jfj/9j ^ y^Hdj Ih{f)~^, k ^m. 

j^k j^k 

That is, 

min Iliif) < , r H9j =Ik{g), m. 

O^j^m IJ'kOk{9k - 9k+l) 

3 

Making the infimum with respect to k, we obtain 

max Ik{g)~^ ^ max IIi{f)~^. 
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One may rewrite maxo^k^m as sup/^g^; on the left-hand side since Ik{g) = oo for 
all A; ^ m + 1. Since g G C =^7, we now have 

inf sup 4(5)"^ ^ inf sup 4(5)"^ ^ sup Ui{f)~'^. 
ge^i keE ge^'j keE jesupp(/) 

Next, making the infimum with respect to / G ^n, we obtain 

inf sup7fe(5)~^ ^ inf sup 4(5)"^ ^ in| sup IIi{f)~^. 
ge^i keE ge^'j keE fe^n iesupp(/) 

The second proof for the inverse inequality is to show that 

inf sup/i(/)"^ ^ Aq. 
fe^'i ieE 

For this, recall the definition 

Ao = inf{D(/) : ||/|| = 1, = for all i > some m:l^m<N + l}. 
Because of 

(ll/ll = = for alH > m : 1 ^ m < iV + 1} 

C (ll/ll = 1, = for alH > m + 1 : 1 ^ m < iV + 1}, 

it is clear that 

A^,""^ := inf{L>(/) : ||/|| = 1, = for alH > m : 1 ^ m < + 1} | Aq 

as m t Note that Aq™'' is just the first eigenvalue of the Dirichlet form 
{D, S){D)) restricted to {0, 1, . . . , m} with Dirichlet (absorbing) boundary at m-|- 
1. Now, let g = g^'^^ be the eigenfunction of Aq™"* > with qq = 1. Extend g 
to the whole space by setting gi = for all i > m. By using Proposition 2.2, it 
follows that g G ^'j with supp (g') = {0,1,... ,m}. Furthermore, by (2.9) with 
h = Xog, we have Ii{g)~^ = Ag™^ > for all i ^ m, and hence, 

supli(c/)"^ = sup Ii{g)~^ = A^™^. 

ieE i^m 

Thus, 

A^""' = suplj(5')"^ ^ _ inf sup/i(/)"^ ^ inf sup/i(/)"^ 

ieE fe^'p supp if )={0,l,...,m} ieE fe^'r ieE 

The assertion now follows by letting m ^ N. 

(h) Prove that inf^^j^^^ sup^g^^pp (y.) Ihif)''^ ^ inf^eri s^Pig^ i?i(u). 



SPEED OF STABILITY FOR BIRTH-DEATH PROCESSES 23 

Let u with supp (u) = {0,1,... ,m} be given such that Vi := Ui+i/ui G 'fi. 
Then, the constraint 

Vi <1- ai{v~\ - 1)6^^ ^m, Vm = 0, 

is equivalent to minQ^i^^ Ri(v) > 0, and the constraint 

> aj+i(aj+i + 

comes from the requirement that Vi > for all i < m. Since the case of i = m 
in the first constraint is contained in the second one, we obtain the constraint 
described in In particular, we have 

ai(ai + bi)~^ <vo<l — ao[vZl — l) = 1 

and so vq G (0, 1). By induction, we have Vi G (0, 1) for all i <m. The existence 
of such a n is guaranteed since m < oo, as will be shown in proof (i) below. Now, 
let 

{ai + hijUi - aiUi-i - hiUi+i, i^m, 
0, i > m. 



Then by assumption, fi/ui = Ri{u) > for i ^ m. Hence, / G ^ij. Next, we 
have 

< ^ fXkfk = Hbj{uj - Uj+i), j ^ m. 

Hence, 

m 

i^j XI A^fc/fc = Ui- Um+i = Ui>0, i ^m. 

j=i k^j 

Therefore, we obtain 
and then 

sup i?j (ti) = max (u) = sup IIi{f)~^ ^ inf sup IIi{f)~^. 
i^E i^'m jesupp(/) fe^ii iesupp(/) 

To be consistent with the convention of Ri{v), here we adopt the convention: 
Ri(u) = — oo for all i > m. The assertion now follows by making the infimum 
with respect to u. 

(i) Prove that inf^^^^ supj^^; Ri{v) ^ Aq. 

As in the last part of proof (g), denote by g (with go = 1) the eigenfunction 
of Aq"*"* > 0. Then supp (51) = {0,1,... ,m}, and g is strictly decreasing on 
{0, 1, . . . , m} by part (2) of Proposition 2.2. The definition of g gives us 

bi{9i - 9i+i) - aiigi-i - gi) = xi"'^gi, i ^ m, g^+i = 0. 
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That is, 

aJl — — 1 = A5 % z ^ m. 

Let Vi = Qi^i/gi for i ^ m and i;^ = for i > m. Then Vi G (0,1) for i G 
{0, 1, . . . ,m — 1}, and Ri{v) = Ag"^ for all i ^ m. It is now easy to see that 
V e We have thus constructed a u{= g) required in proof (h). Clearly 
Ri{v) = —00 for all i > m. Therefore, 

Ao"''' = max Ri{v) 

^ _ inf max Ri{v) 

vefi: supp(v)={0,l,... ,m-l} 0^i<m 

^ _ inf sup Ri{u) 

vE'f'i: supp (u)={0,l,... ,n} for some n^O iEE 

= inf supii!j(i;). 

Letting m N, wc obtain the required assertion. 

We have thus completed the circle argument of (2.24)-(2.29) and then the 
proofs of Theorem 2.4 and Proposition 2.5 are finished. □ 

Before moving further, we mention a technical point in the proof above. Instead 
of the approximation with finite state space used in Part II of the above proof, 

it seems more natural to use the truncating procedure for the eigenfunction g. 
However, the next result shows that this procedure is not practical in general. 

Remark 2.6. Let g be the eigenfunction o/Aq > and define g^"^"^ = gl^m- 
Then 



min IliU'''^) = — 

upDfo("»)') An 



9m 



i e supp Aq 

In particular, the sequence { min^gg^pp (^(m)) i7i(5*^"*-') }^^^ may not converge to 

Xq"^ as qo. 

Proof. Note that 



^ m 

min Iliiq^"^^ = min — 7^ i^iT^ UkQk 



= mm ]_^,^l^MiL_JLl±ll (by (2.5)) 



O^i^m Oj ' An 
j=i 

min -^{gi- gm+i) 

O^i^m Aogi 

1 



^ _ 9m+l 



9m 

This proves the main assertion. For Example 3.4 in the next section, we have 

lin, fl_^')=l_ <1, 

\ 9m J \ 
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and so 

lim mill IIi(g^"'A < Xo-\ □ 

Tn->cxD j g supp (gt"")) 

To conclude this section and also for later use, we introduce a variational 
formula of Aq in a different difference form. 

Proposition 2.7. On the set y := {v : Vi > 0, ^ i < A'"}, redefine 

Ri{v) = Oj+i + bi- ai/vi-i - bi+iVi, i e E, v-i > is free, 

where a^+i = b^+i = and vn is free if N < oo. Then 

(1) we have 

sup inf Ri{v) ^ Ao. (2.31) 

The equality sign holds once YliLo fJ-i = If^ this case, we indeed have 
Ao = sup inf Ri{v) = sup inf Ri{v), 

y^-y ieE ^^Y^ ieE 

where 

% = {v: Vi-i > ai/hi, ^ i < AT + 1}. 

(2) In general, we have 

Ao = sup inf Ri{v), (2.32) 

and the supremum in (2.32) can he attained. 

Proof, (a) First, we prove that sup^g^^ infigE Ri{v) ^ 0. Given a positive, non- 
increasing /, /jv+i = if iV < oo, define 

Ui = {liihy'^^Hjfj e{{),oo), i<N + l. 

Then 

biUi — QiUi-i = fi > 0, i e E, u-i > is free. 
This implies that {vi := ■Uj+i/uj : i < N) e "V*. As before, we also use 

i?j(u) := aj+i + 6j - , i^E 

Ui Ui 

instead of Ri (v) . Clearly, 

R.(u) = " ^'^^ ^0, ieE. 
Hence inf^^^; Ri{u) ^ and the required assertion is now obvious. 
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(b) By (a), without loss of generality, assume that Aq > 0. Then by Proposi- 
tion 2.2, the corresponding cigcnfunction g of Ao is positive and strictly decreasing. 
With Ui := gi — ^j+i > (i G E), the eigenequation 

-^g{i) = biUi - aiUi_i = Xogi, i G E, c/^+i := if A/" < oo, 

gives us {vi := Ui-\-i/ui : i < N) e y*. Next, by making a difference of —Q,g{i) 
and —^g{i + 1) and noting that flg{N + 1) is setting to be zero if iV < oo, we 
obtain 

(oi+i + bi)ui - bi+i-Ui+i - aiUi-i = X^Ui, i G E. 

Thus, we have Ri{v) = Aq for all i £ E. Therefore, (2.31) holds. 

(c) To prove the equality sign in (2.31) whenever J2iLo l^i = cso, in view of Part 
I of the proofs of Theorem 2.4 and Proposition 2.5 and (b), it suffices to show 
that 

sup inf Ij(/)~^ ^ sup inf Ri{v). 

In view of (a), without loss of generality, assume that inf^g^; i?i(n) > for a 
given u > 0. Define /j = hiUi — UiUi-i for i E E, /jv+i = oi N < oo. Then it is 
clear that 

{fi - fi+i)/ui = Ri{u) > 0, ieE. (2.33) 

Hence, / is strictly decreasing. 

We now prove that / G whenever /Xj = oo. First, we have 

'^l^kfk = '^l^k{bkUk - akUk-i) = '^{i^kbkUk - i^k-ih-iUk-i) = l^Aui > 0, 

k^i k^i k^i 

i G E. (2.34) 
In particular, /o > 0. If fk^ ^ for some ko ^ 1, then /fcp+i < and 
^ IJ'ifi < /feo+i l^i^ as n ^ OO 

since X^^o/^* ~ This implies that 

fXifi — >■ — oo as n — >■ oo. 

Now, by (2.33), we would get 

< UnKun = XI = m l^^-f^ + XI ~^ as n ^ oo, 

which is impossible. Therefore, / > and then / G 
Combining (2.33) with (2.34), we obtain that 
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Hence, we have first 

inf Ri{u) = inf liif)-^ ^ sup inf 
ieE ieE fe^i 



and then 



sup inf Ri{u) ^ sup inf liif)' 



as required. We have thus proved the equahty in (2.31) under '^j^/J'i = oo. 

Actually, we have proved in the last paragraph that [fi =)hiUi — aiUi-i > 
for all i e E and so {vi := Ui^i/ui) G whenever inf^g^; Ri{u) > 0. This means 
that the set 1^ \ % is useless since for each v £ y\y^,, we have infj^^; Ri{u) ^ 0. 
Now, because of C 1^ and (a), using the equality in (2.31), we obtain the last 
assertion of part (1). 

(d) To prove part (2) of the proposition, note that the inequality is proved 
in (b). For the inverse inequality, recalling that the main body in proof (c) is to 
show that the function fi (i G E) defined there is positive, this is now automatic 
due to the definition of y*. The equality sign in (2.32) has already checked in 
proofs (a) and (b) in the cases Aq = and Aq > 0, respectively. □ 

Remark 2.8. For the equality in (2.31), the condition J2il^i — °° cannot be 

removed. For instance, consider the ergodic case for which - /Xj < oo but Aq = 
by Theorem 3.1 below and so (2.31) is trivial. However, as proved in [3; The- 
orem 1.1] (cf. Theorem 6.1 below), the left-hand side of (2.31) coincides with 
another eigenvalue (called Ai) which can be positive. In this case, the equality in 
(2.31) fails. This also explains the reason for the use of 'K. 

Remark 2.9. The test sequences with the same notation (vi) used in Theo- 
rem 2.4 and Proposition 2.7 are usually different. Corresponding to the eigen- 
function (^j) of Aq, the sequence constructed in proof (d) of Theorem 2.4 and 
Proposition 2.5 is Vi = gi+i/gi, but the one constructed in proof (b) of Proposi- 
tion 2.7 is 

^ _ ffi+l ~ 9i+2 _ 1 — ffi+2/ffi+l 

9i-9i+i 

Thus, the mapping from the first sequence to the second one is as follows: 

{viWi<N ^ (^^P^) , (2-35) 

where on the right-hand side, vn is set to be zero if AT < 00. 

3. Absorbing (Dirichlet) boundary at infinity: 
criterion, approximating procedure and examples 

This section is a continuation of the last one. As applications of the variational 
formulas given in the last section, a criterion for the positivity of Ao and an 
approximating procedure for Aq are presented. The section is ended by a class of 
examples and then the study on the first case of our classification is completed. 
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Theorem 3.1 (Criterion and basic estimates). The decay rate Aq > ijf 
S < oo, where 



n ^1 

6 = sup /x[0, n] u[n, N] = sup \ fij 2, t ■ (^-l) 

j=0 k=n 

More precisely, we have {4:S)~^ ^ Aq ^ . In particular, when N = oo, we have 
Ao = i/ the process is recurrent (i.e., oo) = oo) and Xq > if the process is 
explosive (i.e., condition (1.2) does not hold). 

Proof, (a) Let (pn = Y^^=n^o ~- ^['>T',N],Uj = (hjiij)~^. To prove the lower 
estimate, without loss of generality, assume that (po < oo. Otherwise, 5 = oo 
and so the estimate is trivial. Next, let M„ = ju[0,n] := 'Yllt=Ql^k- By using the 
summation by parts formula 

n n — 1 n 

^XkVk = XnVn - "^Xkiyk+i - Vk), := ^Xj, (3.2) 

fc=o fe=o i=o 

in viewing the definition of S and using the decreasing property of cp, we get 

n n— 1 

j=0 k=0 
^ S 



/C^n ^ ^k 



Noting that 
we obtain 

Therefore, 



^ ^ 25 



1 2(5 2(5 

-fn(v^) ^ — , / , N • —;= = —;= (v^ + y/'^n+i ) ^ 4(5. 

By part (2) of Theorem 2.4, we have Aq ^ (45)~^- 

(b) Next, fix arbitrarily n <m and let fi = v[iy n,m] Then / G =^7. 

To compute Ii{f), note that when i < n or i > m, wc have fi — fi^i = but 
^j^iHfj ^ Afo/o > 0; and when n ^ i ^ m, we have fi - /j+i = z^j = {biiJ,i)~'^. 
Hence, we have 



liif) 



00 (by convention, 1/0 = 00), otherwise. 



SPEED OF STABILITY FOR BIRTH-DEATH PROCESSES 



29 



Clearly, achieves its minimum at i = n, 

inf Ii{f) = /x[0,n] z^[n,m]. 

Since n,m{n < m) are arbitrary, by letting m ^ N and making the supremum 
in n, it follows that 

sup inf Ii{f) ^ sup /i[0, n] i/[n, AT] = d. 

By using part (2) of Theorem 2.4 again, we obtain Aq ^ S~^. Note that in this 
proof, we do not preassume that S < oo. 

(c) The particular assertion for the recurrent case is obvious. The explosive 
case is also easy since 

oo n 

OO > fii ^[i, oo) > jii ^[i, oo) > /x[0, n] v[n, oo) 
i=0 j=o 
for all n, and so (5 < oo. □ 

The next result is parallel to [7; Theorem 2.2], and is a typical application of 

parts (2) and (3) of Theorem 2.4. It provides us a way to improve step by step the 
estimates of Aq- In view of Theorem 3.1, the result is meaningful only if (5 < oo. 

Theorem 3.2 (Approximating procedure). Write Vj = {iijhj)~^ and tpi = 

(1) Whenifo < oo, define fi = /„= fn-JIifn-i) and6n = sup^^E IIi{fn). 
Otherwise, define 5„ = oo . Then 5n is decreasing inn ( denote its limit by 
5oo) and 

Ao^C=s---^'5r'^(4<^)"\ 

where 6 is defined in Theorem 3.1. 

(2) For fixed i,m G E, i < m and m ^ 1, define 

/(^''-) = l^„/i!:7)//(/(!:7)), n^2, 

where l^m 'is the indicator of the set {0, 1, . . . ,m}, and then define 
S'^= sup rnin//,(/(^'-)). 

e,m:e<m ^^f" 

Then S!^ is increasing inn ( denote its limit by S'^ ) and 

Next, define 

II f(^."i)||2 

Then ^ Aq, ^n+i ^ 5'^ for all 1 and 5i = S'^. 

As the first step of the above approximation, we obtain the following improve- 
ment of Theorem 3.1. 
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Corollary 3.3 (Improved estimates). We have 

where 

Si = sup ^ 



(3.3) 



keE 



sup 



3/2 



k=0 V i+l^k<N+l 



eeE '^t 



keE 



teE 



k=l+l 



G [<5, 25]. 



(3.4) 
(3.5) 



Proofs of Theorem 3.2 and Corollary 3.3. (a) First, we prove part (1) of Theorem 
3.2. Noting that if (/?o = oo, then 5 = co and Sn = oo for aU n ^ 1, the assertion 
becomes trivial in view of Theorem 3.1. Thus, we can assume that <^o < oo- 
By (2.23), we have 

Si = sup//i(/i) ^ sup/i(/i). 

ieE ieE 

Proof (a) of Theorem 3.1 shows that the last one is bounded from above by AS. 
This gives us the lower bound of S^^ as required. 

We now prove the monotonicity of {Sn}. By induction, assume that /„ < oo 
and Sn <oo. Then fn+i < oo- Note that 

^Hjfn+lU) =YlJ'jfnU)fn+lij)/fn{j) 
j^k j^k 

< SUpJ/i(/„) V//j/„(j) 
j^k 

= Sn^fijfnU)- 

Multiplying both sides by Vk and making a summation of k from i to A^, by (2.14), 
it follows that 

fn+2{i) ^ Snfn+lii)- 

Because Sn < oo and fn+i{i) < oo, we obtain fn+2 < oo and IIi{fn+i) ^ Sn < oo. 
Now, making the suprcmum over i, we obtain Sn+i ^ Sn. < oo. 
We have thus proved part (1) of Theorem 3.2. 

(b) To prove the monotonicity of S'n given in part (2) of Theorem 3.2, we use 
the proportional property twice: 



k^j j=i k^j 



^min^/.,/^-)(fc)/5;M.er(fc) 



k^i 



k^i 



^rnin/C^'-)(z)/eT\0- 
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This implies that (5^_,_i ^ S!^. 

By part (2) of Theorem 2.4, we also have 5'^ ^ Aq ^ for all n ^ 1. The assertion 
that 5n ^ Aq^ is obvious. Next, let / = fn'"^\ Then g := Igupp (/) = 

fn+i^- ^ consequence of (2.30), we obtain 6n+i ^ (5^. 

We have thus proved part (2) of Theorem 3.2 except the last assertion that 
5i=6[. 

(c) We now prove (3.4) and S[ ^ S. By (2.15), we have 

fn+i{i) = l^kfn{k)iy[iV k,N] 

keE 

= X] f^kfnik) ViVk 
keE 

i 

= ifi'^IJ.kfnik) + X flk^kfnik). (3.6) 
k=0 i+l^fe<Ar+l 

In particular, with /i = y/Ip, we get 

i 

f2{i) = (Pi^lJ'kV^+ XI l^k'/J'^- (3.7) 

fe=0 i+l^k<N+l 

From this, we obtain (3.4). 

To prove S'^ ^ 8, we need some preparation. As an analog of (3.6), we have 

/(^_;T^(i) = t{,^^} X ^[i V kM- (3-8) 

k^m 

In particular, 

f!f''^\i) = ^ V£,m]z/[i VA;,m]. (3.9) 

Since the right-hand side is decreasing in i for i ^ i < m, fi^'"^\i) = fi^'"^\'^) for 
all i ^ £, and fi^'"^\i) = for i > m, it follows that 

min//,(/(^'-))= min I/,(/f'"^') 

^ m 

j=% k^j 
m I m 

= X X V ^, m] / X 

^ min yUfe z/[A; V m] = inf i'j(/}^''"^) 



k^i 
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Here in the last step, we have used the proportional property. Since the sum on 
the right-hand side is increasing in it is clear that 



min > /ifc ^'[A; V m] = m] /Xfc = /x[0, m]. 

We have thus proved that 

d[ = sup min/Ij (/-l^'™'-') ^ sup //[O,^] m] ^ sup /Li[0,£] ipi = 5. 

A different proof of this is given in proof (d) below, 
(d) We now compute S[. Note that by (3.9), we have 

m 
k=0 

Since f2 '"^\i) is decreasing in i and f[^'"^\i) is a constant on {0,1,... ,£}, it 



when £ ^ i ^ m, we have 



is clear that mm,^^ f^''^> {i)/ f^'"'"'' (z) = min,^,^„ Besides, 



k=0 e+l^k^i i+l^k^m 

It follows that 

Jl {'') ^ ^ L ;.^o k=e+l ' ^ k=i+l 

We show that the sum on the right-hand side is increasing in i. That is, 

m 

^-^ iyh,m\ ^-^ 

^+i<fe<j ^ ' ^ fc=j+i 

i+i ^ 

fe=€+l L I I J i+2^k^m 

Collecting the terms, this is equivalent to 

1 / 1 1 \ 

-fii+ii'[i+l,m\'^ ^ iJ,i+iu[i+l,m\+{ V /ifci/[A;, m]^. 

i'[i,m\ \iy[i + l,m\ i^i^T^lJ ^^—^^ 

Now, the conclusion becomes obvious because by the decreasing property of m] 
in i, the first term is controlled by the second, and the last one is nonnegative. 
We have thus obtained that 



= ^[^'-] E/^'^ + E ^^>^-[kM'■ (3.10) 

/i fc=o ^ ' k=e+i 
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As will be seen soon that the right-hand side is increasing in m (> hence, we 
obtain 



S[ = sup min — r — = sup 



• fc=0 e+l^k<N+l 



(3.11) 



Prom this, it follows once again that 6'i ^ 6. We now turn to prove the monotone 
property: 

2 m+l 

fi[0,e]u[£,m + l] + — — V fiiu[i,m + l]^ 

iy[£,m + l\ 



Equivalently, 

This becomes obvious since the term in the last bracket is positive: 



(e) To show that (5^ ^ 2(5, assume ^ < oo. By using the summation by parts 
formula (3.2) with Xk = fik, Xk = J2j=o Vk = ^Ivi^ we get 

M M-1 



k=0 k=0 

M-1 

= V^mXm + X Xfe [ipl - ipl_^_-i] 

k=i 
M-1 

= Vm^m + X XkVkiVk + Vk+l) 
k=i 
M-1 

< VmXm + 2 X XkVkVk 

k=i 
M-1 

^ 5ipM + 25 X i/fc (since X^ifk ^S), i < M < N + 1. 



k=i 

If iV = GO, letting M ^ iV, it follows that 

N 



fc=0 
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The same conclusion holds in the case that N < oo since 

iV-l N 

S(fN + 2(5 ^ Uk < 25 ^ Uk = 26ipi. 

k=i k=i 

Hence, 



1 ^ 

5[ = sup — V Mfe¥'fcvi ^ 25. 
ieE (p-- ^ 



fc=0 

(£,m) / .X 



(f) Now, it remains to compute 5\. Since /{ ' (i) = V^, m]l{j^^}, we 
have 



l|/i II =^^^^iy{i^ ^,mfts^i^^} = [i^,l\v%mf ^ ^ /x^I/[^,m]^ 

and 



m — 1 

E /ij6i(z/[i + l,m] - 2/[i,m]) + /Xm V^^^ 



m—1 



v\i, m]. 



Thus, 



/■(^,m)||2 ^ m 



Hence, we have returned to (3.10). Since the right-hand side is increasing in m 
as we have seen in the proof of (3.11), we obtain 

||^(£,m)||2 

5\ = sup — ^ = 5\. □ 

.<mZ?(/(^-)) 

To conclude this section, we present some examples to illustrate the power of 
our results. The first one is standard having constant rates. 

Example 3.4. Let bi = b > 0{i ^ 0), Ui = a > 0{i ^ 1), b> a. Then 

(1) Aq = {^/a — VbY with eigenfunction g: 

9n=(^^^ (^n + l-n^lj, n^O, g ^ L\^i) U L\fi) . 

(2) S = b{b - a)~^, 5[ = (a + b){b - o)~^ = 5i > 5, and 5i = Aq ^ which is 
exact. Note that 5i/d[ < 2 whenever a ^ b and limh^aSi/6[ = 2. When 
a = b, we have Aq = 5^^ = 5'i ^ = 0. 
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Proof, (a) First, we have /x„ = {b/a)'^, n ^ 0. Hence, 

^/X„=00, ^^„5r„ ^ ^ ^ 1 = OO. 



Next, since 



we have 



and then 



1 (a 



Hihi b\bj ' 



b — a\b 



I 

a ' 



CXD OO ^ OO 

'uiifi = OO. 



Hence, (1.2) holds. It is easy to check that (2.12) holds: 

OO ^ OO / ^\n/2 / 

J]/x„5„i.[n,oo) = — — ^ ( ^) (n + l-nJ^ 

n=0 n=0 V / \ » 







(b) To study Aq, according to (a), the Dirichlet form is regular and so the 
condition "/ G in the definition of Aq can be ignored. Thus, 

Ao = inf Dij) = 6 inf y //»(/,+! - jif. 
11/11=1 11/11=1 

It suffices to consider the case that 6=1. Write 7 = 6/a > 1. Then we have 
and the required quantities are reduced to 

, _ (V7-1)' . _ 7' . _ 7 e, _ 7(7 + 1) 

'^"(7^' ^^"(^^' ^^"W^- 

Now, to prove part (1) of Example 3.4, write ^ = (^/7 — 1)^7"^ for distin- 
guishing with Aq. Since (5,^) satisfies the eigenequation, applying anyone of the 
variational formulas for the lower estimate given in Theorem 2.4 with /j = gi or 



5i vfeV 1 + 1(1-^/^) J V l + ^(l-7-'/')y' 



it follows that Aq ^ We have seen that the equality sign holds once g G L^{^). 
Unfortunately, we are now out of this case. Therefore, we need to show that 
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Ao ^ To do so, one may use the truncated function of g: gf^^ = gi'^{i^m}- 
Then by the Stolz theorem, we have 



An ^ hm 



115 



(m) 112 



hm 

m— >oo 



bm + 



Mm— l^m— 1 
Mm 



1 - 



25 



m— 1 



9r, 



(3.12) 



The last Umit equals ^. Alternatively, noting that the leading order of g ^ 
is 7"^^/^, one may adopt the test function fi = z"*/^ for z > 7. Then / € L^ilJ-)- 
The required assertion follows by computing D(/)/||/|p and then letting 2; i 7. 
This proof benefits very much from the explicitly known expression of Aq. 
(c) The computation of S is easy: 



6 = sup ipn^ 1^3 = _ SUP7" 



»+i(y»+i _ 1) 



J=0 



n>0 



7 



(7-1)2- 

(d) To compute 5i, by (3.7), we have 

i oo 

k=0 k=i+l 

_ 1 [- -i+l V-,fc/2+l/2 , fc/2+3/2\ 

1)3/2 IT Z^T + 7 j> 



k^i+1 



7 



-i/2+3/2 



Therefore, we obtain 

f2{i) 7 / rr-. , ^\ 7 ^ J_ 

Ao' 



(5i = sup 



^Wi(0 (7-l)(V7-l) 

Noting that even if neither /i nor /2 is the eigenfunction, we still obtain the sharp 
estimate. 

(e) To compute S[, by (3.5), we have 



d[ = sup 

eeE 

1 



^ 1 

</5£VMfc + — V fj-k^pl 

^ — ' (Z'/' ^ — ' 



■ sup 



7 



-€+1 



7 - 1 £eE 
1 



7 



fe=0 



k'^e+i 



( T^^ sup [7^-7 ^+'+7] 

7(7 + 1) 
(7-1)2- 



□ 



The next example is a typical linear model for which, interestingly, we have a 
very simple and common eigenfunction. Moreover, the eigenvalue Aq is determined 

by the constant term 27 in the rates, but not the difference of the coefficients of 
the leading term i, as in the ergodic case (cf. Example 6.8 below). 
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Example 3.5. Let bi = 2{i + 7) > 0), 7 > 0, = z (z ^ 1). Then 

(1) Ao = 7, 9n = 2"" for n^O, and g e L'^ifi). 

(2) When 7 = 1, we have (5 = log2 « 0.69, 5[ ?a 0.84, and di ?a 1.09. Then 
6i/6[ « 1.3 < 2. 

Proof. The uniqueness condition (1.2) is trivial since the birth rates are hnear: 

00 ^ k oc - ^ ^ k — 1 - 00 



(a) Because 



2"7(l + 7)---(n-l + 7) . ^ 

Mo = 1, Mn = j , n ^ 1, 

n! 



it follows that fin > 72" /n and so X^„Atn = 00. Next, since 

n! 

we have X]n(Mn^n)~^ < Furthermore, we have 

2-"7(l + 7) • • • (n - 1 + 7) 



71' 

n=0 n=0 



The ratio test tells us 51 G L^(|u). Since Aq is explicit and g G L'^{p,), it is simple 
to check that is the cigcnfunction of Ao- Hence, the proof of part (1) is done. 
For this example, the sequence (vi) takes a simple form: Vi = 1/2. 
(b) When 7 = 1, we have Aq = 1, 

k^i+l 

In particular, tpo = log 2, ipi = log 2 — 1/2. Numerical computations show that the 
supremum in the definition of S, S[ and Si are attained at 0, and 1, respectively, 
and moreover, 

6 = (fono = (/^o = log 2 Ri 0.69, 



S[ = fof^o + —72 t^k^l = log 2 + V 2 Vfe ~ 0.84, 
Vo — log 2 ^ — ' 



k>2 



2 log 2 - 1 + 1 V(2 log 2)(2 log 2 - 1) + J ^J^ _ ^ ^ 2 V^^' 



21og2 _ 

^ k^2 

PS 1.09. 

We have thus proved part (2) of the conclusion. □ 

The next example is often used in the study of convergence rates. For which, 
the first eigenfunction is unknown but Aq can still be computed. 
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Example 3.6. Let bi = {i + if and ai = . Then 5 = tt^/G ^ 1.64, S'l ^ 2.19, 
and (5i = 4 which is sharp (Xq = 1/4^. Besides, Si/6[ ^ 1.83 < 2. 

Proof, (a) Since Hi = 1, Ui = {i + 1)~^, we have fi[0,i] = i + 1 and (fi = 
■ ^ 6'i, the supremum is attained at 0, therefore, 



^0 



V— - — 



and 



^ oo 



« 2.19. 



fc=0 

(b) For Si, the supremum is attained at oo and is equal to 4. By Corollary 3.3, 
this means that Aq ^ 1/4. This can be also deduced by part (1) of Theorem 2.4 
with Vi = 1 — (2i + 4)"-*^ for which the minimum of Ri{v) is attained at i = and 
i = oo. It is even more simpler to use Vi = 1 — {2i + 3)~^. Next, it is known 
that Ao ^ 1/4 (cf. Example 5.5 below), hence, the estimate is sharp. A direct 
proof for the upper estimate goes as follows. Since the lower estimate is sharp, it 
indicates to use the test function 



E 



U + 1)^ 



1/2 



Vi + T' 



However, the last function is not in L'^{ij,), and so one needs an approximating 
procedure. Now, a carefully designed test function is the following: 



(a) _ 



y/{i + l)a^+i ' 



a > 1. 



Then 



{i + 1)q;»+i 



1=0 

oo 



oo I 

^ log[a(a-l)-i]< 



oo, 



i=0 



oo .9 



i=l 
1 



V(^ + 2)a^+2 ^{i + l)ai+i 



LV(^ + l)a Vi. 



[ii + l)a 



^ (z + l)a^+i y{i + l)a + ]2 

1 °° 1 

i=l ^ ' 

= i(2 + log[a(a-l)-i]). 
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The required assertion now follows from 



2 + log[a(a-l)-i] ^ 1 ^ , ^ ^ 

Ao ^ —, — p-^ -r—r. > - as a 4, 1. □ 

41og[Q;(a — 1) 4 



The last example below does not satisfy the non-explosive condition (1.2). 

Example 3.7. Letbi = (i + l)'^ and ai = i{i — l/2){i'^ +3i + 3). Then ^ ■ /li < oo, 
< oo, Ao = 1/2, 5 1.83, 5[ Ri 1.9, and 5i 2. Moreover, Si/S[ Ri 1.05 < 

2. 

Proof. A simple computation shows that 

^ i!^ ^ ^ nt=i(fc-l/2)(P + 3fc + 3) 

nL=i(fc-V2)(fe2+3A; + 3)' (i + l)(i + l)!3 

From this, it follows that Ylit^i < °° ^i^i < as an application of the 
typical Rummer's test: for a positive sequence {xn}, Yin converges or diverges 
according to k > 1 or k < 1, respectively, where 

K= lim n( -^-l). (3.13) 

n^oo \Xn+l J 

For each of S, S[ and Si, the supremum is attained at 0. 
To see that Aq = 1/2, first we check that Ri{v) = 1/2 for 

'^ = '"2(iTiy 

This gives us Aq ^ 1/2 by part (1) of Theorem 2.4. Since the corresponding 
eigenfunction g, 

Vt^ (2i-l)! 

k=0 ^ '' 

decreases strictly to and Y^if^i < '^^ have g G L'^i/J')- Now, because —0,g = 
Ao5, 5oo = 0, and D{f) = -{g, Ug), it follows that Aq = 1/2 by (2.18). □ 

4. Absorbing (Dirichlet) boundary at origin 

AND REFLECTING (NeUMANN) BOUNDARY AT INFINITY 

This section deals with the second case of the boundary conditions. The process 
has state space E = {i : 1 ^ i < N + 1} {N ^ oo), birth rates 6j > but 6jv = 
if < oo, and death rates > 0. The rate oi > is regarded as a killing from 
1. Define 

Ao = inf{Z?(/)//x(/2) : / / 0, Z)(/) < oo}, (4.1) 



40 



MU-FA CHEN 



where n{f) = Y,keEl^kfk, and 

D{f ) = ^Hkak{fk- fk-iT, /o:=0, 

keE 

6i • • • bk-^ 

/xi = l, lik = - —, 2^k<N + l. 

02 • • • Ofe 

The constant Aq'^'^^ describes the optimal constant C = Aq ^ in the following 
weighted Hardy inequality: 

^ CD{f), fo = 

(cf. [9]). In other words, we are studying the discrete version of the weighted 
Hardy inequality in this section. To save the notation, in this and the subsequent 
sections, we use the same notation Aq, /, //, R and so on as in Section 2. Each 
of them plays a similar role but may have different meaning in different sections. 

To study Ao, as in Section 2, we need some parallel notation originally intro- 
duced in [3, 7]: 

1 ^ 1 ' 1 ^ 

Here, for the first operator, we adopt the convention: /o = 0. The second one can 
be re-written as 

Next, define 



' k=l j=i ^^^^ 



Ri{v) = ai[l- + bi(l - Vi), i e E, vq := 00 

{vn is free if AT < 00 since 1^ = 0) and 

=^/7 = {/:/. >0 for alHG^}, 

= {/:/> and is strictly increasing on 
ri = {v:vi>l for all i G E}. 

The modifications of and are as follows: 

= {/ : there exists m e E such that fi = fiAm > for i e E}, 
= {/ : there exists m £ E such that fi = fi/^m > ^oi i £ E and / is 
strictly increasing in {1, . . . , m}}. 

Here, we use again the convention: 1/0 = 00. Note that for the localization, / is 
stopped at m rather than vanishing after m used in Sections 2 and 3. This is due 
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to the fact that the Neumann boundary is imposed at m but not the Dirichlet 
one. Besides, for the operator II here, the restriction on supp (/) used in Section 
2 is no longer needed. Finally, define a local operator R (depending on m) acting 
on 

^1 = ^meE{v ■.l<Vi<l + ai{l - vr\)br^ for i = 1,2, . . . ,m - 1 
and Vi = l for i ^ m} 

by replacing a^, with := l^mC^m/ Yl^=ml^k in Ri{v) for the same m as in Yx- 
Again, the change of is due to the Neumann boundary at m. Note that if 
Vi = \ for all i ^ m, then Ri{v) = Ri{v) = for all i > m. 

Before stating our main results in this section, we mention an exceptional case 
that Hi = oo. On the one hand, by choosing /o = and /j = 1 for i ^ 1, it 
follows that 

D{f) = Hiai < oo, = ^ /ij = oo 

and so Aq = 0. On the other hand, if ^^^^ fJ-i < oo, then for every / with 
= DO, by setting Z^'") = f./,m G ^^(yw), we get 

m 

oo>Z)(/(-)) =^^,ai(/i-/,_i)2 t D{f) asm^oo, 

i=l 

m 

OO > ^)^Y1 ^ifi ^ OO = ^i{f ) as m ^ OO. 

In words, for each non-square-integrable function /, both iJi{P) and D{f ) can be 
approximated by a sequence of square-integrable ones. Hence, we can rewrite Aq 
as follows: 

Ao = inf{i^(/) : = 1}. (4.2) 

In this CctSG, ctS will be seen soon but not obvious, we also have 

Ao = inf {£>(/): /x(/^) = l,/i = /iAm for some meE and all ieE], (4.3) 

Besides, we mention that the Dirichlet eigenvalue Aq is independent of 6o ^ (cf. 
[4; Theorem 3.4] or [12; Theorem 3.7]). 

For a large part of the paper, we do not use the uniqueness condition (1.2) 
(note that a change of a finite number of the rates and bi does not interfere 
in the uniqueness). Under (1.2), the process is ergodic iff ^^fJii < oo (see [10; 
Theorem 4.45 (2)], for instance). If (1.2) fails but N = oo, then the decay rate 
for the minimal process is delayed to Section 7. In (2.2), the condition "/ G ,J^" 
means that we deal with the minimal process. This condition is removed in (4.2). 
It means that we are in this section dealing with the maximal process in the sense 
that the domain &"^^'^{D) of D ignored in (4.2) is taken to be the largest one: 
{/ G L'^iP') ■ D{f) < oo} (that is the maximal process described at the beginning 
of Section 6 but killed at 1). When N = oo, even though there is now a killing 
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at 1 (i.e., oi > 0), the regularity for (or the uniqueness of) the Dirichlet form is 
still equivalent to (1.3): 

f;f-i-+A.fc') =oo (1.3)' 

since a modification of a finite number of rates does not change the regularity (cf. 
Theorem 9.22 for further information). In this section and Section 6, starting from 
any point in E, even though the process can visit every larger state, it will come 
back in a finite time. In this sense, the point infinity is regarded as a reflecting 
boundary. 

It is the position to finish the comparison of (4.1) and (4.2). We have seen 
that A^^-') = aJ"-') once X]iA*i < °o. We now claim that they can be different 
otherwise. To see this, note that on the one hand, A^^-^^ = if Ei 

fii = DO, as 

proved above. On the other hand, once (1.3)' holds (in particular, if Ylil^i = oo, 
then) by Proposition 1.3, Aq^'^^ coincides with 

inf{Z)(/):/G jr, //(/')= 1}, 

which is the one used in (7.1) below and can often be non-zero. Thus, in general, 
A^"-') ^ Xi^-'^ and they can be different. As will be seen in Theorem 7.1 (2), in 
the special case that both of the series in (1.3)' are divergent, we have Aq'^'^^ = 

a(-) = 0. 

Theorem 4.1. Assume that YliLil^i < oo. Then the following variational for- 
mulas hold for Ao defined by one of (4-1) — (4-3)- 

(1) Difference form: 

inf sup Ri{v) = Ao = sup inf Ri{v). 

(2) Single summation form: 

inf sup/i(/)"^ = Ao = sup inf Ii{f)~^, 

(3) Double summation form: 

Ao = sup inf Iliif)''^ = sup inf IIi{f)~^, 

Ao = inf sup//i(/)"^ = inf sup//i(/)"^ = mf _ sup//i(/)"\ 
fe^i ieE fe^ii ieE f^^jju^'^^ ieE 

where = {f : fi > for all i e E and fll{f) G ^^(^u)}. 

The next result was proved in [6] except the exceptional case that ^ ■ /i^ = oo 
in which case Aq = (and 5 = oo) and so the assertion is trivial. See also 
Corollary 5.2 below. Note that {vj) below is different from (2.15). 
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Theorem 4.2 (Criterion and basic estimates). The rate Aq defined by (4-1) 
(or equivalently by (4-^) provided ^^^e l^i ^) positive iff S < oo, where 



S = sup n] uln, N] = sup fij. 

^ 1=1 J=n 



neE 



(4.4) 



More precisely, we have (45) ^ ^ Aq ^ 5 ^. /n particular, we have Xq = if 
^i^E l^i — ^ '^^^ Aq > i/ either N < oo or (1.3)' fails. 



Theorem 4.3 (Approximating procedure). Assume that X]il=i Mi < oo and 

6 < CO. Write (po = 0, (pi = i^[l,i] := Yl]=i(.l^j^j)~^ j i ^ E. 

(1) Define fi = /„ = fn-\II{fn-\) and 5^ = sup^ge //i(/„). Then 5^ is 
decreasing in n and 

Ao ^ ••• ^ (4<5)-^ 

(2) For fixed m £ E, define 

/("^) = <^(-Am), 

/^)=[il-Wil-l)](-Am), n^2 

and i/ien define 5'^ = sup^gg inf jg^; //^(/i,™''). T/ien 5^ is increasing in 
n and 

^ s[-^ ^ • • • ^ C ^ Ao. 

Next, define 

(m) 2\ 



(5„ = sup 



neE. 



Then S^^ Aq, ^n+i ^ /'''^ all n ^ 1 and Si = S[. 

As the first step given in Theorem 4.3, we obtain the fohowing improvement 
of Theorem 4.2. 

Corollary 4.4 (Improved estimates). For the rate Aq defined by (4-1) {or 
equivalently by (4-^) provided Y^^^^l^i ^ °°)' have 

where 

5l = sup Hk^PiAkV^ 



ieE J'-Pi 



N 



sup f^k^i^^ + v^Piy^ukVvk 

^ iV r ^ m-1 
= sup ^ l^k'PlAm= sup ^ Hkfl + fPml^lm, N] 



■ k=l 



G [5, 2<5]. 



(4.5) 
(4.6) 
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Proof of Theorem 4.1. Note that 

1 ^ 1 ^ 

Hence, /»(/) coincides with used in [3, 4, 6, 7, 12], whenever 6o > 0. The 

same change is made for the operator //(/) in this section. 

Throughout this proof, wc use Aq = Aq"^'^"* to denote the one given in (4.3). 
Similar to the proofs of Theorem 2.4 and Proposition 2.5, we adopt the following 
circle arguments: 

Ao ^ A^^-') (4.7) 
> sup inf //,(/)-!= sup inf//,(/)-i= sup inf /,(/)-! (4.8) 

^ sup inf Ri{v) (4.9) 

^ Ao (4.10) 

and 

Ao ^ mf ^ sup//i(/)-i (4.11) 

^ inf sup//i(/)-i= inf sup//i(/)-i= inf sup/i(/)-i (4.12) 

^ inf sup Ri{v) (4.13) 

veTi ieE 

< Ao (4.14) 

Assertion (4.7) is obvious. The following assertions are proved in [4; Theorem 3.3], 
or [12; §3.8] and [7; §2] (see also the remark given in the next paragraph): 

sup inf/i(/)-i= sup miIIi{f)-^= sup inf /Ii(/)-i ^ Af,*''^. (4.15) 

inf sup/i(/)-^ = inf sup//i(/)-i = in| sup 11 i{f)-\ (4.16) 
fe^i ieE fe^i ieE fe^n ieE 

In particular, we have known (4.8) and (4.12) since the inequality in (4.12) is 
trivial. It remains to prove (4.9)-(4.11), (4.13) and (4.14). 

In [7; §2] and [12; §3.8], only the ergodic case under condition (1.2) is consid- 
ered. But for (4.15) and (4.16), one does not need (1.2). Actually, one can now 
follow the proofs of Theorem 2.4 and Proposition 2.5 with a little change. For 
instance, to prove the last inequality in (4.15), following proof (a) of Theorem 2.4 
and Proposition 2.5, let g satisfy \\g\\ = 1 and go = 0. Then 



^ = ^l^i9l (since ll^ll = 1) 
ieE 

= ^t^i[^(9j-9j-in (since go = 0) 

ieE ^j=l ^ 



ieE j=i ^ k=i '^'^ 
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Exchanging the order of the first two sums on the right-hand side, we get 




=: D{g) sup Hj. 
jeE 



The next step is to choose hj = l^'ifi for a given / G J^u with sup^g^; Hjif) < 
oo. From these, it should be clear what change is required in order to prove (4.15) 
and (4.16). 

We now begin to work on the additional part of the proof. 

(a) Prove that supj^^^^ infig^ //i(/)"^ ^ sup^g^^ infig^ Ri{v). 

As in proof (c) of Theorem 2.4 and Proposition 2.5, we use Ri{u), 

R,iu) = ajl- ^) +bjl- ^] , ieE, 

\ Ui J \ Ui J 

(uat+i is free if AT < oo since fejv = 0), instead of Ri{v), where Ui > for i e E 
and uo = 0. Then Vi > 1 {i G E) means that Uj+i > Uj > 0, and Vi = 1 for i ^ m 
means that Ui = Ui^m > 0. 

Without loss of generality, assume that infi^E Riiu) > for a given strictly 
increasing u with uq = 0. Define fi = (oj + bi)ui — Oj-Uj-i — ftj-Ui+i [ = UiRi{u)] 
for i e -E and /o = 0. Then by assumption, 

fi/ui = Ri{u) > 0, ieE. 

Hence, / G Next, since 

< iJ-kfk = fJ-kakiuk - Uk-i) - fJ-k+iak+iiuk+i - Uk) 
and the strictly increasing property of Uj in i, it follows that 

N 

< ^^kfk ^ tJ-jdjiuj - Uj-i), 

k=j 

and so 

i i N 

Ui = '^{Uj - Uj-i) ^ ^ i^j ^ fJ-kfk > 0. 
j=l j=l k=j 

We obtain 

Ri{u) = fi/ui ^ ni{f)-\ ieE. 
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Therefore, we have first 

inf Riiu) ^ inf Iliif)-^ ^ sup inf IIi{f)-\ 
ieE ieE fe^ii 

and then 

sup inf Ri{v) ^ sup inf IIi{f)~^, 

as required. 

(b) Prove that sup^^^^ inf^gE Ri{v) ^ Aq- 

First, we show that sup„g^^ infjee Ri{v) ^ 0. For a given positive / G L^ilJ), 
let u = fll{f). Then Ui+i/ui > 1 and Ri{u) = fi/ui > for all i G E. With 
{vi = Ui^i/ui) G 1^1, this implies inf Ri{v) ^ and then the required assertion 
follows. 

Alternatively, since oi > 0, the eigenfunction is still strictly increasing when 
Ao = by part (3) of Proposition 2.1. Hence the proof in the case of Aq = 
can be combined into the next paragraph, and then the last paragraph can be 
omitted. 

By assumption, we have "YlieE l^i ^ When Aq > 0, it was proved in proof 
(d) of [12; Theorem 3.7] that the eigenfunction of Aq"*'^-* is strictly increasing. Even 

though Ao could formally be bigger than Ag^'^-* , the same proof still works for the 

eigenfunction g of Ao since the modified function g used there satisfies gi = gi/\n 
for some n. Having this at hand, the proof is just a use of the eigenequation: 

-ng{i) := -bi{gi+i - gi) + ai{gi - gi-i) = Xogi, ie E, go:=0 

{9n+i f^^® if AT < oo since = 0). With Vi := gij^i/gi > 1 for z < iV, this 
gives us G and Ri{v) = Ao, and so the assertion follows. 
We have thus completed the circle argument of (4.7) — (4.10). 

(c) Prove that Aq ^ inf^^^^^^^^^,^ sup^g^; //i(/)"^ 

In the original proof of [7; Theorem 2.1], when N = oo, from the estimate 

D{g)^fi{g'') suplliif)-' 
ieE 

for / G and g := [///(/)](• A m) to conclude that Ao ^ D{g)/ fi{g'^), one 
requires an additional condition g G L'^{fi), provided m = oo is allowed. This is 
the reason why the set in part (3) of Theorem 4.1 is added. Anyhow, with 
the modified conditions, the same proof gives us the required assertion (cf. proof 
(f) of Theorem 2.4 and Proposition 2.5). 

(d) Prove that inf^^j^^^ supjgg //i(/)-^ ^ inf^g^i ^^PieE Ri{v)- 

Given u with uq = Q and Ui = UjAm for all i G .E so that {vi := Uj+i/wj) G Yi, 

let 

^ _ 1 (oi + bi)ui - ttiUi-i - biUi+i, 

I am{um - Um-i), i > m. 
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It is simple to check that /o = 0, 

fi/ui = Ri(u) > for i G {1, . . . ,m} and fi = fm for i > m, 
and so / G Moreover, since 

m—l 

fJ-kfk = HCljiUj - Uj-l) - (J-mamiUm - Um-l) 

N 

= tijaj{Uj - Uj-i) - /m X] 

k=m 

we get 

N 

< ^fikfk = fijajiuj - Uj-i). 

k=j 

It follows that 

i i N 

< Mi = ^{uj - Uj-i) = XI XI « G {1, ■ ■ ■ , m}, 

j=i j=i k=j 

and then Ri{u) = fi/ui = IIi{f)~^ for i G {1, 2, . . . , m}. Therefore, we have 

max Ri{u) = max IIi[f)~^^ _ inf max IIi(f)~^^ inf sup/Ij(/)~\ 

and then ^ 

inf sup ^ inf sup IIi{f)~^. 

(e) Prove that inf^^^^ sup^g^; Ri{v) ^ Aq. 
Recall the definition of Aq: 

Ao = inf {£)(/) : = 1, = f^^^ for some m G -B and all ieE}. 

Clearly, we have 

A^™) := inf {£>(/) : fi{f) = 1, fi = fi^m for alH G £;} 4 Aq as m t iV. 
We now explain the meaning of Aq'"'' as follows. Let 

N 

jli = IJ.i, 1 ^ i < m, /x„ = X f^i^ 

i=m 

m 

D{f) = Y,ftMfi-fi-if- (4.17) 
j=i 
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Then fiiCii = /Xjaj tor i = 1, . . . ,m, D{f) = D{f) and = for every 

/ with / = /.Am- Thus, Aq™^ is just the first eigenvalue of the local Dirichlet 
form (15,^(1?)) having the state space {!,... ,m}, with Dirichlet (absorbing) 
boundary at and Neumann (reflecting) boundary at m. Let g {go = 0) be the 
eigenfunction of the local first eigenvalue Ag'"''. Extend g to the whole space by 
setting gi = gi/^m- Next, set Uj = gi for i < N. Then 

0, i > m. 

Furthermore, for Vi := Ui+i/ui, we have vq = oo, > 1 on {1, . . . , m — 1}, and 
= 1 for i ^ m. Thus, by (4.18), it is easy to check that v = (fj) G Y^. Therefore, 

Aq™"^ = max Ri{v) 

^ _ inf max Ri{v) 

vefi: Vi=l for j^m 

^ _ inf sup ii!i(v) 

v&'f\: Vi=l for j ^ some n>l i&E 

= inf sup Ri {v) . 

The assertion now follows by letting m ^ N. □ 

Proof of Theorem 4.3. 

(a) We remark that the sequence {/n™''}^^^; is clearly contained in ^j. But 
the modified sequence used in [7; Theorem 2.2], 

/^) = ^(-Am), /(-) = /a(-Am)7/(/a(-Am)), n ^ 2, 

is usually not contained in However, 

sup inf //,(/(-)) 

= sup min 

= sup min //,(/('")(• Am)) 
= sup inf //i (/(-)(• Am)). 

Here in the last step, we have used the convention 1/0 = oo. Hence, these two 
sequences produce the same {S!^}. 

(b) The approximating procedure given in Theorem 4.3 is mainly a copy of [7; 
Theorem 2.2] (cf. the proof of Theorem 3.2). For later use, here we review the 
proof of part (1). From [6; proof of Theorem 3.5], we have known that 
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Hence (Alternatively, by the proportional property), 

i -. i 

/2« = Yl ^ 4<5^(/i(i) - - 1)) = 4(5/i(z). 

j=l 3 k^j j = l 

This gives us the assertion Si = sup^^j^ Hiifi) ^ 4(5. 

To prove the monotonicity of {5n\ and {/n} C L^ifJ-), we adopt induction. As 
we have just seen, 

5i=sup^^45. 

i>i /iW 

This means that /i G L^ifj) (or equivalently, J2 < 00) and 5i < 00 since 5 < 00 by 
assumption. Assume that /„ G L^ifJ-) (or equivalently, fn+i < 00) and 5„ < 00. 
Then 

^l^kfn+lik) =Yl^kfn{k)[fn+l{k)/ fn{k)] ^ K^l^kfnik). 
k^j k^j k^j 

Multiplying both sides by uj and making summation from 1 to i, it follows that 

fn+2{i) ^ Snfn+l{i), i ^ 1- 

Since fn+i < oo and 5„ < oo by assumption, we have fn+2 < oo, and 

fn+2{i) 



/n+l(i) 



Iliifn+l) < OO, i^l. 



This proves not only fn+i & L^if^) but also 5n+i ^ Sn < oo. 

The assertion that ^ Aq is obvious by (4.2). Similar to proof (b) of 
Theorem 3.2, the assertion Sn+i ^ (5^ is a consequence of the last part of the 
proof of [7; Theorem 2.1]. □ 

Proof of Corollary 4.4. 

(a) The degenerated case that X] ■ /U^ = oo is trivial since Aq"^'^^ = and 5 = 5i = 
S'l = oo. The main assertion of Corollary 4.4 is a consequence of Theorem 4.3. 
Here, we consider (4.6) only since the proof of (4.5) is easier. Note that 

-H'i(/^"'^) = ^ ^ fJ'k^kAm- 

j=l r J J 

The right-hand side is clearly increasing in i for i ^ m and is decreasing (not 

hard to check) in i when i ^ m. Hence, Ilj {f["^^) achieves its minimum at i = m. 
Then, by exchanging the order of the sums, it follows that the minimum is equal 
to 



1 ^ 

— y2 IJ'kVkAr 



k=l 
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This observation is due to Sirl, Zhang and PoUett (2007). We have thus proved 
the first equality in (4.6). 

Next, following the proof of [6; Theorem 3.5], we have 



and 



N 



Ml 



(/r') = E 



k=l 



Combining these facts together, it follows that 5i = 5^. 

(b) Finally, we prove the estimates in (4.6). The lower estimate of 5'i is rather 
easy since 



N 



1 

— Z^l^kVkAm ^ 



A' 



N 



'~Pm 



(frn — 

k = l k=m k=m 

For the upper estimate, use the summation by parts formula: 



N 

E 



f^kVkAr 



N 



1 \ ^ \ ^ ^k + V'fc-l \ ^ 

j=k k=l J=k 



N 



k=l 

It follows that 

AT 

L ■ 



k=l 



1 A 2 2 ^ 1 r ^ 1 25 ^ 1 

> fJ'k'PkAm < yf^j ^ > 



25. 



k=l '^'"fe=l'^"- j=k fc=l 

The estimate now follows by making the supremum with respect to m G .E. □ 

5. Dual approach 

This section is devoted to the duality of the processes studied in the previous 
sections, as well as a duality to be used in the next two sections. Again, the 
section is ended by a class of examples. 

Suppose that we are given a birth-death process with state space E = {i : ^ 
i < N + 1}{N ^ oo), birth rates 6^ > (6o > 0, especially) but ^at ^ if iV < oo, 
and death rates > but oq = 0. The case that > is used in this section 
while the case of 6iv = is for use in Section 7. Define a dual chain with state 
space E = {i : 1 ^ i < N' + 1} and with rates as follows: 



6o = 0, bi = ai, ai = bi-i, ieE, 
where ajv+i = bN+i = if < oo by convention and 



(5.1) 



N' 



N, N < oo and b^ = 0, 

+ 1, A^ < oo and > 0, 
oo. A" = oo. 
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The dual process with rates (ai,6j) has an absorbing at 0. When N < oo, for 
the dual process, the state + 1 is absorbing if 6iv = (then oat+i = but 
Sn > 0); otherwise, it is a reflecting boundary since cin+i = bf^ > 0. In a word, 
the absorbing boundary is dual to the reflecting one and vice versa. This dual 
technique goes back to Karlin and McGregor (1957b, §6). Next, define 

Ai = l, An = " " " , 2^n<N' + l. (5.2) 

When N < oo and > 0, then ajv+i > 0, and so /i„ can be defined up to 
n = N + 1. Otherwise, it can be defined up to n = iV only. It is now easy to 
check (noticing the difference of {uj) and (Oj)) that 

/i„ = = feol^n-l, K ■= = T-fJ-n-i, l^n<N' + l. (5.3) 

Actually, the rates (a^, 6^) in (5.1) are determined by the transform given in (5.3): 
An = bo^n-i and On = iin-i/bo- From this, it follows that 

A^n = 6oZ^n+l = Ol^'n+l, ^^n = r~An+l = :^An+l, ^ n < TV', 



(J'N = aif^fiNbN) ^ 'd N < oo and 6jv = 0, (5-4) 



and so 



N' ^ N' ^ N'-l 

V = H '^n = 7- V ^n, 

f^nan bo 

N' N' N'-l ^ 

^ An = f'O ^ Z^n-1 =boYl -Tir- 
n=l n=l n=0 

Note that by (5.1), 

Ui f . bi ^ 
ttj+i + bi Oi+iVi = 6j+i + Oj+i aj+2 

Vi-l Vi-i 

By a change of the variables (vi) G y-. 

bi+i 



(5.5) 



(5.6) 



O-i+2 

or 

Vi = Vi-l = — Vi-l, (5.7) 
bi ai 

we get 

Oj+i + 6j — bi+iVi = fij+i ( 1 - ^ ) + 6i+i(l - Uj+i). 
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Since 60 > 0, f_i > but oq = 0, from (5.7), it is clear that we should set vo = 00. 
Next, by (5.7) again, 

Vi-l > — ■<=^ Vi > 1. 

bi 

It remains to examine the boundary condition on the right-hand side when N < 
00. 

(1) First, let Bn = 0. Then v = {vi > : ^ i < N — 1) , v-i and vn-i are 
free. The dual state space is E = {1,2, . . . ,N}. The dual test sequence 
is V = [vi > : 1 ^ i < N) , vm = 0. 

(2) Next, let b^- > 0. Then v = {vi > : ^ i < N), V-i and arc free. 
The dual state space is E = {1,2,... ,A'' + 1} with reflecting at + 1. 
Hence, v = [vi > : 1 ^ i < N + l) , vn+i = 0. 

We have thus proved the following result. 

Proposition 5.1. For the dual processes defined above, the following identities 
hold: 



sup inf 



ai+1 + bi 



Vi-l 



bi+iv, 



sup inf 

5 l^i<N'+l 



dii 1 



Vi-l 



+ bi{l-Vi) 



(5.8) 



where v = {vi > : ^ i < N' — 1) with free V-i, and v = (vi > : 1 ^ i < N') 
with vq = 00, vn'-i is free and vn' = if N < 00; 



sup inf 



sup inf 



o-i+i + bi- 
aj 1 



a,; 



Vi-l 

1 



Vi-l 



bi+iVi 
+ bi{l - Vi 



(5.9) 



in the case that 6jv > i/iV < 00, where 'fl is given in Proposition 2.7, and Yi 
is defined in Theorem 4-1 replacing N by N + 1 when N < 00. 
In these formulas, ajv+i = 6jv+i = if N < 00 by convention. 

Corollary 5.2. Given rates {ai,bi) as in Section 2 (then bjss > Q if N < 00), let 



An 



A 



(2.2) 



and define 6 by (3.1). Next, define the dual rates [ai,bi) as above. 
Correspondingly, we have Aq o-iT'd S defined by (4-1) o,nd (4-4) replacing N by 
N + 1 if N < 00, respectively, in terms of the dual rates. Then we have Aq = Aq 
and 6 = 6. 

Proof. Having relationship (5.9) at hand, the assertion that Aq = Aq follows by 
a combination part (2) of Proposition 2.7 and part (1) of Theorem 4.1, provided 
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Next, by (4.4), (5.3), and (3.1), we have 



sup fij 

i N+1 

sup V Uj V Afc 



i iV+1 



sup 



bovk- 



N 



sup 



0^i<N+l'—Q 



k=i 



fJ-kbk 



= s. 



This proves that 6 = S. In particular, if X^j/ii( = X^j z^i) = oo, then by Theo- 
rem 3.1 and Corollary 4.4, we get Aq = Aq = 0. We have thus completed the proof 
of Ao = Aq. □ 

As will be seen in Theorem 7.1 (2), in the degenerated case that /^i = oo 
and '^i{fJ'ibi)~^ = oo, the dual of the process studied in Section 2 also goes to 
the one studied in Section 7. 

Before moving further, let us discuss the duality used here. Very recently, Chi 
Zhang provides us a nice explanation which leads to a deeper understanding of 
the duality (5.1). Consider a simple example as follows: 

-ai - bi 



Q 



f-bo 
ai 


V 



02 






bi 

— 02 — 62 

as 






62 

-as -bsJ 



ai,bi > 0. 



Introduce an invertible matrix: 



M-- 



fJ'obo 


-f^obo 











Hibi 


-Hibi 











y^2^2 


-IJ'2b2 











IJ-sbs 



M- 



/ ^ 



V 



1 




1 
1 

M262 




Then 



MQM- 



f 

V 



ai - 


bo 


Ol 











bi 




-02 - 


bi 


02 












b2 




-03 - 


b2 


os 












b3 




-6s/ 


ai - 


bi 


bi 


h 







\ 


0.2 




-0.2 - 


h 












0-3 




-os - 


bs 


h 












04 




-04/ 



1 

Ms 63 

1 

Ms 63 
1 

Msbs 



Q. 
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Hence, the dual matrix Q is just the classical similar transformation of Q and 
so they have the same spectrum. In particular, the eigenequation Qg = —Xog 
{g 7^ 0) is transferred into 

Q{Mg) = {MQM-^){Mg) = X^Mg = \o{Mg). 

Hence, the eigenfunction g of Aq is transformed io g = Mg of Aq = Aq. Corre- 
spondingly, the test function / is transformed to / = Mf. Prom this, it should 
be clear that all the operators R and R, I and /, // and // are closely related to 
each other and then so are the variational formulas. 

Having these facts at hand, one can simplify a part of the previous proofs. 
However, we prefer to keep all the details here since they are needed when we go 
to the more general situation, so called the Poincare-type inequalities (Section 8), 
or can be used as a reference for studying the continuous case. For the Poincare- 
type inequalities, the current duality seems not available. 

By Corollary 5.2, we have two ways to estimate Aq = Aq: using either the rates 
{ai,bi) or (ai,6i). The corresponding formulas for S[, S[, 5i and are collected 
in Tables 5.1 and 5.2. 

Table 5.1: Expressions ol 6 = S, 5'i, 6[, 6i and 6i in terms of the rates (6j,aj): 



N 



S = 5= sup ij.[0,i]u[i, N] = sup / ^ f^j/ , ^k-, 
0^i<N+l o<i<Af+ip^ ^ 

1 ^ 

5[= sup -^—Ttiki^[kyi,Nf 



(5.10) 



sup 

0<i<iV+l . 



N 



/x[0,z]z/[z,Ar] + 



u[i,N\ 



k=i+l 



(5.11) 



N 



^1= sup --^Vi/fc/i[0,feA z]2 



fc=o 



sup 

0^i<JV+l . 



i-1 



//[0,i]z/[i,iV] + -i-^^z/fcM[0,fe]' 



' fc=0 



(5.12) 



Si= sup V fiki^[i V k, N] y/iy[k, N] 

0^i<Ar+l y^ull, N] ^ 

- i ^ 



sup 

0^i<Ar+l 



N 



k=Q 
N 



(5.13) 



Si = sup — ■j= 

o^i<Ar+i V^[0,i] 



'Yukn[0,kAi]\/iJ.[0,k] 



sup 



1 



N 



^kl^[0, k]^/^ + VMM "kVt^^k] 



k=0 



k=i 



(5.14) 



SPEED OF STABILITY FOR BIRTH-DEATH PROCESSES 55 
Table 5.2: Expressions of S = 6, 6'i, S'^, Si and in terms of the rates (6j,aj): 



N 



6 = S= sup i'[l,i]il[i,N] = sup V'^^feV'/ij, 

Ki<iV+l l^i<iV+l^ ^ 

1 ^ 

(5i = sup V i>kfi[k V i, A^]^ 



sup 

l^j<iV+l 



AT 



J k=i+l 



1 ^ 

^1 = sup ^qrj— rV] AfcZ>[l,fcA 



sup 

l<i<JV+l . 



i-1 



1 ^ 

(5i = sup -7== V i>fe A[fe V i, A^] VMMV] 



sup 

l^i<iV+l 



N 



k=l 

N 



sup 

l<i<JV+l 



i-1 



fe=i 



(5.15) 



(5.16) 



(5.17) 



(5.18) 



(5.19) 



The next four examples are dual of Examples 3.4-3.7, respectively. 

Example 5.3. For Example 3.4, we have cii = b{i 1), 6j = a(o > 0), 6 ^ a. 

Then 6 = 5 = b{a- b)'^ , 5[ = S[ = (a + b){a - b)'^ , and 5i = 5i = Aq ^ = 
[^/a — Vb) . In particular, if we take a, = 4 and 6j = 1 (i ^ 1), then Aq = 1, 

5[ = 5/9 = 0.5, S2 = 0.64, 63 ^ 0.71, ^4 w 0.755, 6'^ « 0.79; 
^1 = 0.5, 62 ~ 0.71, S3 «i 0.79, S4 « 0.835 ^5 0.8647. 

Thus, S'^ and Sn are increasing and close to Xq^ as n^. 



Proof. To compute S[ and Si, we use Table 5.1. For simplicity, write 7 = 6/a > 1. 
Then 

7 — 1 
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(a) Note that 



i—l oo 

' -I fc=0 fe=i 



1 

7-1 v-.-fc/T 



,7 



fc=0 



k=i 



7-1 



+ 



7 



1 



6(7-1) 
1 



fc=0 k>i 



7 

' fc=0 

7(1+7) 2(i + l)7' 



7-1 



7 



i+l _ 1 



Since the second term in the last [■ ■ ■ ] is negative and 7 > 1, the right-hand side 
attains its supremum at i = 00. By (5.12), we have thus obtained 



7(1+7) a + b 



6(7-1)2 {a -by 



(b) Next, note that 



i-l 



k=i 



1 



i-l 



V 7* 1 fc=o fc^i 

^ V 'y('=+=^)/2 , ^ i+l _ I ^-fe/2+1/2 

^/^i+l - 1 ^ ^ 

V / fc=o 



k>i 



1 ^3/2^y/2 _ ^-i/2+l/2^^i+l _ I 



^7^+1 - 1 v^- 1 
7 _^ 7^/7 



1-1/V7 



V7-1 V7-1 



6(7-1) 
1 

6(7"^ 
1 

1 

6(7"^ 

= 7 
6(V7-1)' 
1 

" i^-Vby' 

By (5.14), this means that Si ^ Xq^ and so the equahty sign must hold because 
S^^ is a lower estimate: Aq ^ 

(c) We now compute the approximating sequences {S'^} and {^n} for the upper 
estimate, using the dual rate (aj,6i). In the particular case, we have 
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The approximating sequences can be computed successively by using the following 
formulas: 

/^^« = ^, iG{l,2,...,m}, 

^ i— 1 m — 1 

/r^o = ^ E(i - + (4^ - 1) E 

fc=l k=i 

+ ^(4^_l)4i— /(-)(^)|, ^G{1,2,... ,m}, n^2. 

Then S'^ = sup^^^ mini^^^^ For the first five of {6'^}, the 
minimum are all attained at m and so the computations become easier. 
To compute Sn, simply use the formula 

- V^"i Ai-i f('^'>( 4)2 , -3-1 Ai-mA^)( \2 

^™,4^-(/^)(^)-/^)(^-l))^ 

Example 5.4. For Example 3.5 with 7 = 1 (6j = z, ai = 2ij, we have 6[ « 
0.75 < 5i 0.84 ond ^1 1.12 > (5i Ri 1.09. Besides, Si/6[ ^ 1.5. 

Example 5.5. For Example 3.6, we have ai = hi = i^ {i ^ 1), 60 = 0, = 2 < 
d[ Ri 2.19 and Si = 5i = 4 which is sharp. Besides, Si/S[ =2. 

Proof. By Example 3.6 and Corollary 5.2, it follows that Aq = Aq = 1/4. Here, 
we present an easier proof for the upper estimate. Note that when = 6j for 
i ^ 2, we have 

/xi = 1, fii = ^ = T^, i ^ 2; fiibi = bi, i^l. (5.20) 

a2 • • • ai 

In the present case, we have /ij = i~^ {i ^ 1) and /tjOj = 1. Let f^"^^ = \/i t\m. 
Then 

z=l 

i=i i=i (V« + V«- 1) ^ i=i ^ 



Hence, 



i)(f(™)) 1 

Ao ^ lim = -. □ 
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Example 5.6. For Example 3.1, we have Oj = {i ^ 1), hi = 1/2) (i^ +3z + 
3), Ao = Ao = 1/2, 5[ ^ 1.83 <5[k 1.9 and 5-i_^5i^ 2. Besides, 5i/K = 1.09. 

Proof. First, we have 

. m=i(fc-V2)(fe^+3fc + 3) (^-1)!^ 

' ' ni=\(fc-l/2)(A:2 + 3fc + 3)' " 

By (5.5) and Example 3.7, we have X^j/ti < oo and X^j < °o, and so the 
minimal dual process is explosive (but here we are dealing with the maximal 
one). The sharp lower bound can be deduced from part (1) of Theorem 4.1 with 
the dual test sequence 

Vi = I -\ — r, ^ 1. 

+ 3z + 3) ' 

From this, it follows that the corresponding eigenfunction 

i-l 

9i = Y[vk, 2, gi = 1, 

fe=i 

increases strictly to a finite limit since Xi>i^~^(^^ + 3i + 3)~^ < oo. The se- 
quence (vi) comes from the one computed in Example 3.7 plus a use of (2.35) and 
(5.7). □ 

The precise value of Aq for the next example is unknown. Its eigenfunction 
is non-polynomial. It is interesting to compare this example with the ergodic 
one given in §6 for which Ai = 2, as well as the one with rates ai = i + 1 and 
bi = {i^ 1) given in §7 for which Aq = 2. 

Example 5.7. Let bo = 0, hi = i + 2 {i ^ 1) and hi = i^. It is the dual of the 

process studied in §S with rates = i + 2 (f ^ 1) and bi = (i + 1)^ (i ^ 0). Then 
Ao G (0.395,0.399), S[ « 2.37 < 5^ « 2.48 and di « 2.63 > di ^ 2.61. Besides, 
dJ5[ « 1.1. 

It is interesting that for all of Examples 5.3-5.7, we have d[ ^ 6[ and ^ 5i 
which then means that Corollary 3.3 is more effective than Corollary 4.4. The 
effectiveness of the bounds 5i and S[ given in Corollary 4.4 was also checked by 
Sirl, Zhang and PoUett (2007) for some models from practice. 

Remark 5.8. It is now a suitable position to mention a method for the numerical 
computation of Aq defined in §4. The idea is meaningful in the other cases. From 
proof (b) of Theorem 4.1, it follows that there is a sequence {vi : Vi > 1,1 ^ i < N) 
such that 

Ri{v) = ai{l - v~}-^) + bi{l - Vi) = Ao, vo = oo, vn = ii N < oo. 
Hence, we have 

vi - 1 = {a, - Xo)bi\. 

Vi-1= [ai{l - i;r_\) - Xo]br\ 2 ^ i < N. 
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In other words, replacing — 1 by Uj, when z = Xq, the equation 

. . , (5.22) 

Ui = [aiUi_i{l + Ui_i) - z\bi , 2 < z < iV, 

has a positive solution (ui = Ui{z))i^i<N ■ Thus, one may use the maximal z 
so that (5.22) has a positive solution as an approximation of Aq (based on part 
(1) of Theorem 4.1). In this way, we obtain the approximation of Aq given in 
Example 5.7. 

6. Reflecting (Neumann) boundaries 

AT ORIGIN AND INFINITY (eRGODIC CASE) 

We now turn to studying the first non-trivial eigenvalue in the ergodic case. 
Let = {i : ^ i < iV + 1} (A^ ^ oo), 6o > 0, 6iv = if iV < oo, 

Ai = inf {D{f) : /.(/) = 0, ^i{f) = l}, (6.1) 

where n{f) = J fd/i, 

D{f)= I^Mfi+i-fi?= I^i^fi-fi-^f (6-2) 

0^i<Ar l<j<iV+l 

with domain ^'^'^'^{D) = {/ G L'^il^) ■ D{f) < oo}. In (6.1), we presume that 

N 



^/ii<oo. (6.3) 



i=0 

Then the Dirichlet form {D, ^™-^{D)) has a trivial eigenvalue Aq = with con- 
stant eigenfunction 1, and here we are working on the next "eigenvalue" Ai of 
(D, ^'^'''^{D)). If (6.3) does not hold, then 1 ^ L'^ilj) and so Ai is not meaningful. 
Moreover, by (1.3) and Proposition 1.3, the Dirichlet form is unique. In this case, 
the corresponding process is explosive, or zero-recurrent, or transient. The decay 
rate is described by Aq which has already been treated in Sections 2 and 3. Hence, 
throughout this section, we assume (6.3). 

Note that condition (6.3) plus (1.2) means that the unique process is ergodic. 
When N = oo and (1.2) fails, the minimal process was treated in Sections 2 
and 3, and in this section, we arc dealing with the maximal process (cf. [10; 
Proposition 6.56]) as in Section 4, it is indeed the unique honest reversible process. 
Denote by Q = {qij) the birth-death Q- matrix. Then under (6.3), the maximal 
process P™^(A) (Laplace transform) can be expressed as 

i^7-(A) = i^7"(A) + ^'^''"''^-^y i.J^E, A>0, 

where for each fixed j, {P™'°(A) : i ^ E} \s the minimal solution to the equations 

\ ^ Qik , ^ij . ^ „ 

Xi = 2_^ T-^—xk + T-— , « e E, 
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and {zi{X) : i E E) is the maximal solution to the equation 

{XI-Q)u = 0, 
^ u ^ 1, 

(cf. [10; Proposition 6.56]). According to a result due to Z.K. Wang (1964) 
(of. Wang and Yang (1992, §6.8, Theorem 2)): if AT = oo and (1.2) fails, then 
every honest process (may be non-symmetric) is crgodic and so is the maximal 
one. Certainly, within the symmetric context, by using (1.4), it is easy to check 
directly the ergodicity of the maximal process. 

Here, we mention a technical point. If (6.3) fails, then as mentioned before, 
by (1.3), there is precisely one symmctrizablc process (Dirichlet form) which is 
nothing but the minimal one. Thus, if (1.2) also fails, then the unique process must 
be explosive and so there is no honest symmetrizable process. This is a different 
point to the reversible case (i.e., (6.3) holds) for which there exists exactly one 
honest reversible process as just mentioned above. 

We use the same notation I, //, J^j, J^j and ^// defined in Section 4 with 
an addition "/o = 0" in the last four sets, but redefine R and y as follows: 

Ri{v) = Oj+i +bi- ai/vi-i - bi^iVi, < N, 

V-i >0 is free and so is fjv-i if N <oo, 
y = {v : Vi > for all i : ^ i < N - 1}. 

The local operator R is modified from R, replacing by := HmO'm/ X^j^rn 
for V with supp (u) = {0, 1, . . . , m — 2} in the set 

7^ 1^1 f o,i+i ai+i + hi — ail Vi-i , 

-y = \v: <Vi< ^ %±_!_^ for i = 0, 1, . . . , m - 2 



m=2 



and V,- = for i > m — 1 



}• 



Theorem 6.1. Under (6.3), the following variational formulas for Ai hold. 

(1) Difference form: 

inf^ sup Ri{v) = Ai = sup inf Ri{v). 

(2) Summation form: 

mf _ sup Ii{f)~^ = Ai = sup inf Ii{f)~'^ , 

where 

^'i = {/ G L^ifJ-) : /o = 0, / is strictly increasing}, 
f = f-T^{f), 7r = /x/Z and Z = ^.^^Hi. 
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The use of / in the last line is based on the property / = / + c for every 
constant c and so we can fix /o to be 0. 

Proof of Theorem 6.1. In the ergodic case under (1.2), the assertion 

inf _ sup Ii{f)~^ ^ Ai 

was proved in [7; Theorem 2.3] (but in the case that A; = oo in the original proof, 
one requires the L^-integrability condition included in as was pointed out in 
proof (c) of Theorem 4.1). The proof remains the same in the present general 
situation with an obvious modification when N < oo. Next, in the ergodic case 
under (1.2), the following result 

Ai = sup inf Ri(v) = sup inf L(f)~^ (6.4) 

is just [3; Theorem 1.1]. In the present general situation, the proof for the second 
equality in (6.4) needs a slight change only (cf. [3; Lemma 2.1]). To prove the 
first equality in (6.4), we claim that 

Ai = inf {D{f) : fi{\f - Tr{f)f) = 1, fi = fi^m for some m e E, m'^ 1} 
=: Ai. (6.5) 

To see this, first it is clear that Ai ^ Ai. Next, the proof of [4; Theorem 3.2] gives 
us 

Ai ^ sup inf /i(/)-\ 

and furthermore, the equality sign with Ai replaced by Ai holds. Once again, the 
key point for the last statement is to show that the eigenfunction of Ai is strictly 
increasing. For this, the original proof needs only a modification replacing Ai 
by Ai (as indicated in proof (b) of Theorem 4.1). Therefore, (6.4) holds in the 
present general situation. 

Now, we need only to show that 

• inf/e#, snpi^i^E h{f)~^ ^ ^^^ver supo<ci<Af Ri{v), and 

• inf^er s^Po=^»<iv Rii^) ^ Ai. 

(a) Prove that inf^g_^^ supi^^g^ /i(/)-^ ^ inf^er suPo^i<iv Rii^). 

As before, write R{u) instead of R{v). Given u with supp (n) = {0, 1, . . . , m — 
1} so that (vi) G 'f, where Vi = Ui+i/ui > for i < m— 1 and Vi = for i ^ m— 1, 
let 

QiUi-i - biUi, i < m, 
amUm-i, i ^ m. 



Since the constraint in Y is equivalent to minj^^_i Ri{v) > 0, it is easy to check 
that _ 

(/i+i - fi)/ui = Ri{u) > for i < m and fi = fi/^m, 
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and so f + boUo E Moreover, since 

AT m-1 N 

fJ-kfk = ^ A*fc/fe + /m ^ /^j 
fc=i fc=i 3=Tn 

m-1 TV 

= X] l^k{akUk-l - hUk) + fm^ fJ'j 
k=i j=m 

N 
j=m 



we get 



and so 



N 



fe=o 



N 



^IJ-kfk = IJ-iO-iUi-i, 

k=i 

N N N 



k=m k=m k=m 

It follows that 



II n - ^— ' 

Hence, 



Ui-l 

Therefore, we have 

max Ri{u)= max Ii{f)~^^ _ inf max Ii{f)~^^ inf sup Ii{f)~] 

0^i<m l^i^m /i=/iAm ^^i<">- fe^i l^ieE 

and then 

inf_ sup Ri{v) ^ inf sup Ii{f)~^. 

Here, we have used the fact that Ri{v) = — oo for i ^ m — 1 if supp (i;) = 
{0, 1, . . . ,m — 2} and in the last step, we have returned to the original notation 
R{v) instead of R{u). 

(b) Prove that inf^^^ supo^j<jv Ri{v) ^ Ai. 
Because of 

W) = 0, l^{f) =1, f^ = f^Am] C {//(/) = 0, = h f^ = /^A(m+l)}, 



SPEED OF STABILITY FOR BIRTH-DEATH PROCESSES 



63 



by (6.5), it is clear that 

A^™) := inf {D{f) : ^i{f) = 0, i^{f) = 1, fi = /.a..} I Xi as m t iV. 

Actually, this is a special case of an approximation result given in [2; Theorem 4.2 
and Corollary 4.3] or [10; Theorem 9.20 and Corollary 9.21]. Note that X^^^ is 
just the first non-trivial eigenvalue of the local Dirichlet form (D, defined 
by (4.17) replacing the Dirichlet boundary at by the Neumann one (having the 
state space {0,1,... ,m}), with Neumann (reflecting) boundary at m. Denote 
by g the first eigenfunction of A^"*"* and extend it to the whole space by setting 
Qi = 9i/\m- Now, if we set ui = ^j+i — gi for i G E, then > for i ^ m — 1, 
Ui = for i m, and furthermore, 

R.(u) = X["^^ > ioT allium- 1. 

Moreover, by the definition of g, we have 

biUi - ttiUi-i = -X^^^gi, i 

Om^m— 1 — '^l 9m- 

Making a difference of this with the one replacing i by i + 1, we get Ri{u) = A^"*^ 
for alH ^ m — 1 (From this, the reason should be clear why in the definition of 
y , we use "■Uj = for z ^ m — 1" rather than "■Uj = ■UjAm")- Thus, 

A^™^ = max Ri{u) 

^ inf _ max Ri{u) 

u: supp(u)={0,l,... ,ra-l}; {vi=Ui+i /ui)er 0^i<m 

^ inf _ sup Ri{u) 

u: supp (■!i)={0,l,... ,n} for some n^O, n<N; {vi=UiJ^i / Uij&y 0^i<N 

= inf_ sup Riiv). 

Here in the last step, we have returned to the original notation R{v) instead of 
R{u). Letting m ^ N, we obtain the required assertion. □ 

With the same rates (oj, 6j) here but endow with the Dirichlet boundary at 0, 
we return to the situation studied in Section 4. The next result, taken from [7; 
Theorem 2.2] and [6; Theorem 3.5], is a comparison of Ai with the quantities Aq, 
5, 6i and 6[ given in Section 4. Sec also Corollary 6.6 below for an improvement. 

Theorem 6.2 (Criterion and basic estimates). Under (6.3), Xi > iff 6 < 
oo. More precisely, we have 

^ ^ ^ Ao ^ Ai ^ Ao^ ^ 1^ ^ f . (6.6) 
4o Oi Oi 

The next two results are mainly taken from [7; Theorem 2.4] with an addition 
on the monotonicity of {r]n} and {rj'^}- 
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Theorem 6.3 (Approximating procedure). Let (6.3) hold and 5 < oo. Write 
(fo = 0,fi = Eo^i^i-i(/^i^i)~^ =• i^[0,«-l] (1 ^ « e E), / = /-7r(/), tt = n/Z, 
and Z = Yli^^j^Hk =■■ l^[0, N] . 

(1) Define fi = /„ = and ?]„ = sup^^.g^ ^i(/n)- ^/len ?7„ 
is decreasing in n and Ai ^ 77"^ ^ (45) ~^ /o'^ all n ^ 1. 

(2) For /ixed mEE: m^l, define 



/^^ = <^.A„., = [^-{//(^-l)] (• Am), n ^ 2, 



p(m) 



p(m) 



and then define 



ri'n 



sup inf liifn), 



sup 



(m) 2\ 



n ^ 1. 



T/ien ry^ is increasing in n and r]'^ ^ ^ fi~^ ^ Xi for all 1. 

The notation ^^fn-iHifn-iY^ used in the theorem may have 0/0 but it should 
not cost any confusion. Note that here we use the same (vj) as in (2.15). In other 
words, when 60 > 0, we use (2.15). But for its dual, it is more convenient to use 
Uj = [fijCLj)'^ as in Section 4 since 60 = 0- This is consistent with the notation 
used in Section 5. 

As a consequence of Theorem 6.3, we have the following improvement of The- 
orem 6.2. 



Corollary 6.4 (Improved estimates). Let (6.3) hold. Then we have 

(4,5)-^ ^ r/r' ^ Ai ^ f7r\ 

where 



(6.7) 



771= sup (V^+ V/^i-l ) 



1 

r]i = sup 



li[i,N] 



N 



1 

sup 



+ 



l^fe^m— 1 

n[m, N] 



l^keE 

1 



fi[0,N] 



2n 



l^k<fk 



l^fc^m— 1 



/x[0,iV] 



ifml^lO, m - 1] - 2 X fik^k 

l^fe^m— 1 



(6.^ 



(6.9) 



Proof of Theorem 6.3. 

Part 1. We prove that {/«} C L^{jx) in three steps. This was missed in the 
original paper [7]. Certainly, we need only to consider the case that N = 00. 
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(a) First, we show that the functions {hn}, 

j J oo 

ho{i) = l, ieE, /i„(z)=^ ^/Xfe/t„_i(A;), z ^ 1, n ^ 1, 

j=l ^^^^ k=j 

are all in L^{p). Clearly, hi (and then /i„ for n ^ 2) may increase to infinity if 
the minimal process is recurrent which is the main problem we need to handle. 
The required assertion says that even though hn can be unbounded but is still in 
L}{ij)- For this, to distinguish with {/„} used in Theorem 6.3, let {/„} be the 
sequence defined in part (1) of Theorem 4.3: 

1/2 ^ 



^ fc=i 



i ^ 1, Vj-i := 



fn{i) = X] '^.J-i X] f^kfn-i{k), i ^ 1, n ^ 2, 

j=l k=j 

/„(0) = 0, n ^ 1. 
From proof (b) of Theorem 4.3, we have seen that 

j=i /^j'^j' 

Because fi{i) ^ /i(l) = a^^^^ ^ ^ !> this gives us 

hi{i) ^4S^/Elfl(i), i^l. 
By induction, it follows that 

This proves that /i„ G L^il-t) for all n ^ 1 since /i G L^ifJ-) as mentioned in proof 
(b) of Theorem 4.3, due to the assumption S < oo. 

(b) Next, we study the relation between {fn} and {/„}. By definition, we have 

i oo 

/2W = ^^j-i^l^kfi{k) = f2(i) - /ii(i)7r(/i), i ^ 1, 
j=i k=j 

i 00 

/3(*) =^^j-l^l^kMk) = /3(i) - /i2('i)^(/l) - hi{i)TT{f2), 1 ^ 1, 

j = l fc=j 

AW = hii) - h3{i)7r{h) - h2{i)TT{h) - hi{i)7r{f3), i ^ 1. 
Successively, we obtain 



n-l 



fn = fn-^ T^{fk)hn-k, U ^ 2. 



k=l 
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(c) Since /i = /i G L^il^) as shown in proof (b) of Theorem 4.3. Now, to show 
that {/„} C L^ilJ.), by (a) and (b), it suffices to prove that {/„} C L^{lj)- This is 
done in proof (b) of Theorem 4.3. 

Part 2. We now prove the monotonicity of {??n} in two steps. Since /„ values 
both positive and negative or even zero, the proportional property used in the 
proof of the monotonicity of {5n} is currently not available. To overcome this 
difficulty, a finer technique is needed. 

(d) Because 



N 



fJ'iaiifnii) - fn{i - 1)] = ^^J'kfn-l{k), n ^ 2, 



k=i 

by the definition of I{f), we obtain 

N I N 

Vn= sup V]Mj/n(j) / y^Atfc/n-iW, n ^ 2. (6.10) 
Since the denominator is positive, the assertion that ?7„ ^ ?7n-i is equivalent to 

N 

XI ^3 [fnU) - r)n-lfn-l{j)] ^0, 1 E E. 

j=i 

That is, 

1 ^ 

»7n-l7r(/„-l) - 7r(/„) ^ X IJ,j [rin-lfn-l{j) - fn{j)] , i^E. (6.11) 

Let us observe the meaning of this inequality: the left-hand side is the infimum 
(attained at i = 0) of the right-hand side. 

The monotonicity of {rjn} now follows once we show that the right-hand side 
of (6.11) is luckily increasing in i, or equivalently, 

N N 

/^[^^] XI ^^j[^n-lfn-l{j) - fnH)] ^ ^[i + I, N]^ Hj[r]n-lfn-l{j) - fn{j)\ ■ 
j=i+l j=i 

By removing the common term 

N 

+ 1,N] X /^J [Vn-lfn-lij) - fnU)] 
j=i+l 

in both sides, it is enough to check that 

N N 

Vn-lYl [fn-l{j)-fn-l{i)] > ^ [/-OV/nW] > i & E, U ^ 2. (6.12) 

j=i+l j=i+l 
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First, let n ^ 3. Then by the definition of /„ and (6.10), we have 

j N 

fn{j) - fn{i) = XI ^s-l^l^kfn-l{k) 
s=i+l k=s 

j N 

s=i+l k=s 
= Vn-1 [/n-l(i) - fn-l{i)] ■ 

This certainly implies (6.12) in the case of ^ 3, regarded as an application of 
the proportional property. Next, let n = 2. Then by the definition of /2 and rji, 
we have 

j N 
f2{j) - f2(i) = ^ Us-l^Hkfl(.k) 
s=i+l k=s 

^miz [his)- his -I)] 

s=i+l 

= rii[h{j)-hii)]- 

This also implies (6.12) in the case of n = 2. We have thus proved that ry„ ^ Vn-i 
for all n ^ 2. 

Part 3. To prove the monotonicity of {?/^}, for each fixed m, as a dual argument 
(exchanging "sup" and with "inf " and , respectively) of the above proofs 
(d) and (e), we have 



Then the assertion follows by making supremum with respect to m. 
Part 4. The proof of fin ^ v'n given in Lemma 6.5 below. □ 

In practice, using fjn rather than rj^ is based on the following result. 

Lemma 6.5. For every non- decreasing, and non-constant function f satisfying 
f € L^i/J-) and D{f) < oo, we have 

^^ffl ^ inf liif). 
D{f) ^ KiGB 

Similarly, for every nonnegative, non- decreasing, and non-zero function f satis- 
fying f G L^in) and D{f) < oo, we have 
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Proof, (a) Since / is not a constant, we have n{f^) > and D{f) > 0. Moreover, 
since / G L^ifi) is also non-decreasing, we claim that 



N 

oo " 

k=i 



> ^ ijLkfk > for all i € i ^ 1. 



Actually, the non-decreasing sequence {fk}, starting at /o < (since / is non- 
trivial) and having mean zero, should be positive for all large enough k. Thus, if 
Yl^=io f^kfk ^ for some io : 1 ^ io <E E, then we would have fi^ < (otherwise 
/i ^ for all i ^ io and then Y.k=io l^^fk > l^jfj > for large enough j). This 
implies that 

and furthermore, 

N 

= /x(/) = X iikfk + X ^kfk ^ X /^fc/fe < 0, 

k—iQ k-^iQ — l k^ii^ — 1 

which is a contradiction. Because of the assertion we have just proved and using 
the convention that 1/0 = oo, it follows that infi^jg^ /,(/) G [0, oo). 
Let 7 = infi^jgE /»(/)• Then we have 

AT 

- XI f^kfk = 'YlJ'kfk^ HJ-iO-iifi- fi-l) 
k^i—1 k=i 

first for those i with fi > and then for all i : 1 ^ i E E. Multiplying both 
sides by /i — ^ 0, we obtain 

-{fi - fi-i) X f^k^k ^ ll^iaiUi - fi-if, ie E, i^l. 

k^i-l 

Making a summation over i from 1 to m, it follows that 

m m 

-^ifi-fi-i) X l^klk^l^l^iaiifi- fi-if- 

i=l k^i—1 i=l 

Noticing that the mean of / equals zero and exchanging the order of the sums, 
the left-hand side is equal to 

mm m 
k=0 i=k+l k=0 



= —fm, X ^fc-^fc + X ^'^fk 

k=0 k=0 
N m 

= X ^^kfkfm + '^fJ'kfk- 
k=m+l k=0 
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As mentioned in the last paragraph, fm> ^ first for some m and then for all large 
enough m since / is non-decreasing, the right-hand side is controlled, for large 
enough m, from above by 

N m 
k=m+l k=0 

With the assumption D{f) < oo in mind, the required assertion now follows 
immediately by passing the limit as m — t- A^. 

(b) For the second assertion, since / G -^^(m) is nonnegative and non-zero, we 
have 

JV 

oo > E Mj/j > foi' all i & E. 

Now, if /o = 0, then there is an io such that /io-i = but /j^ > and so hoif) < 
oo. If /o > and infi^jge = oo, then / should be a positive constant, and 
hence, D{f) = 0. In this case, the assertion is trivial since IJ,{f) > 0. Therefore, 
we may assume that 7 := infi^jg^ Ii{f) < 00. We now have 

N 

^fikfk > llJ'iaiifi - fi-i), ieE,i^l. 

k=i 

Hence, 

m N m 

^ifi - fi-i)^l^kfk > l^l^iaiifi - fi-if. 

i=l k=i i=i 

Exchanging the order of the sums, the left-hand side is equal to 

N kAm N N 

Mfe/fc E ~ /i-l) = E t^kfkifkAm — fo) ^ E l^kfkfkAm ^ 
fc=l i=l k=l k=l 

Combining this with the last inequality, we have obtained the required asser- 
tion. □ 

Having the comparison of f]n ^ fj'^ (Lemma 6.5) in mind, one may expect a 
parallel result for 5'^ and J„ defined in Theorem 4.3. All the examples we have 
ever computed support the conjecture that 5„ ^ S!^, however, there is still no 
proof. In general, wc have 5n+i ^ S'^ only as stated in Theorem 4.3. Note that 

is defined by using 11 {fn) rather than I{fn)- If we redefine 6'^ by using I{fn) 
as in [7; Theorem 2.2], denoted by 5'^ for a moment, then by the second assertion 
of Lemma 6.5, we do have S'^- Besides, by the theorem just quoted, we also 
have S'^^ S'^^ ^n-i- This remark is also meaningful for those S'^ and 5„ defined 
in Section 3. 

Note that the factor of the upper and lower bounds of Ai given in Theorem 6.2 
is 4:Z > 4. The next result has a factor 4 only. A simple comparison of k. below 
and S^^'*^ shows that it is not easy to find such a result. Its proof is delayed to 
the next section. 
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Corollary 6.6 (Criterion and basic estimates). Let (6.3) hold. Then we have 
K~^/4 ^ Ai ^ where 



, ^ \ -In / m-1 

= inf 

:n+i \ 

i=0 



0^n<m<N+l 



j=n 



Furthermore, we have 
where Z = YliLol^i' 



N m N-1 



5l= sup V ^ V /Xj = (5(^-^\ Sr= sup y2l^J^~ir- 

Kn<JV+l ^ Mifli ^ 0^m<JV ^ f^kbk 



In the case that the minimal process is ergodic, since 
1 < Z < oo, — — = oo, 

we have 5r = co and so the second assertion of Corollary 6.6 goes back to Theo- 
rem 6.2. However, the first assertion of Corollary 6.6 is clearly finer. An extension 
of Corollary 6.6 to a more general state space is given in Corollary 7.9 below. 

Most of the examples below are taken from [10; Examples 9.27]. The compu- 
tation of fji, r]i/fji, and k is newly added. 

Example 6.7. Let bi = b{i ^ 0), and Oj = a (i ^ 1), a > b. Then 

Ai = (v^- \/fe)^, 5 = K = a(a - 5)"^, fj-^^ = d[ = {a + b)/{a - bf , 

and rji = X^^ which is sharp. Besides, rji/fji ^ 2, the equality sign holds iff 
b = a. Note that Xi, rj^^ , and fji^ all tend to zero as b ^ a. Furthermore, 

(r/i,??i) C (k,4k). 

Example 6.8. The typical linear model: let bi = /3ii + /5o (/3o > 0, /3i ^ 0), 
and Oj = 7ii (71 > /3i) for i ^ 0. Then Ai = 71 When Pq = 0, we have 

Ar=7i-A. 

Example 6.9. Let bi = b/{i + 1) (6 > 0) for i ^ 0, = a > for i ^ 1. Then 
Ai = a - + 4ab - b) /2. 

Example 6.10. Let bi = b{b > 0) for i ^ 0, Oi = {i A k)a (a > 0) for i ^ 1 and 

some k ^ 2 satisfying 

k-^ ^a/b^k{k-l)-^. 

Then Ai = {Vak - Vb)^ . 
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Example 6.11. Let bo = 1, hi = i, and ai = 2i, i 1. Then Ai ^ Aq = 1 but 
the precise value is unknown. Moreover, 

fji « 0.55, 771 « 0.9986 and rn/fji « 1.82 < 2. 

Besides, k ^ 0.4856 and so (771,771) C {k,4:k). 

The next one is a continuation of [6; Example 3.10]. 

Example 6.12. Let E = {0, 1}. Then Ai = ZXq = fj^^ = and Aq = . 
Hence, the last upper bound in (6.6) and the one in Corollary 6.6 are sharp but 
is not an upper bound of Xi. 

The first lower bound in (6.6) and the one in Corollary 6.6 are sharp for the 
seventh example in Table 6.1 below. 

Examples 6.13. Here are some additional examples, given in Table 6.1, for 
which the quantities 771 ^ X'^^ ^ 771 and k ^ A^^ ^ 4k are compared. For all 
these examples, we have (771,771) C (/c, 4k) and so the estimates given in Corol- 
lary 6.4 are better than the ones in Corollary 6.6. 

Table 6.1 Exact Ai and its estimates for eight examples 



bi (i > 0) 


Oi (i ^ 1) 






m 






i + 1 


2i 


1 


^ 0.8 


^ 1.48 


^ 1.85 


2/3 


i + 1 


2i + 3 
oo = 


1/2 


^ 0.346 


^ 0.638 


^ 1.84 


f5i 0.28 


i + 1 


2i+4+V2 
Oo = 


1/3 


« 0.218 


^ 0.398 


^ 1.83 


0.18 


(z + l)-i 


1 


2(3 - v/5)-i 
« 2.618 


PS 1.92 


« 3.24 


« 1.69 


« 1.6 


1 


i^2 


(V2-l)-2 
5.8284 


3 


fa 5.8284 


^ 1.9 


2 


7 + 2 


i^ 


1/2 


« 0.47 


« 0.85 


f« 1.81 


« 0.47 


60 = 1 




4 


2 




2 


1 


2 + (-l)^ 
, 7-V33 
2 


2[2+(-l)^] 
ao = 


(6 - V33)-i 
^ 3.9 


« 2.11 


« 4.21 


2 


« 1.56 



7. Bilateral absorbing (Dirichlet) boundaries 

This section deals with the fourth case of boundary conditions. It consists 
of two parts. The first one is for the ordinary birth-death processes as studied 
in the previous sections and the second one deals with the bilateral birth-death 
processes with a more general state space. 



72 



MU-FA CHEN 



First, let us consider the processes with state space E = {i: l^i<N + 
1} {N ^ cx)) with Dirichlet boundaries at (ai > 0) and Ar + 1 if iV < oo. Similar 
to Section 2, define 

Ao = inf {D{f) : (.(f) = 1, / G JT}, (7.1) 

where the symmetric measure {/li) is the same as in Section 4, = '^^^e f^kfk, 
and 

fcG-B 

with domain ^™™(D). Clearly, if one changes only the boundary condition at 0, 
then the resulting Aq is bigger or equal to Aq ' . Note that if (1.3) fails, then the 

eigenvalues Aq'*'^^ and Aq^'^-* are different which correspond to the maximal and 

the minimal Dirichlet forms, respectively. However, as mentioned in Section 4, 
once (1.3) holds, Aq^'^"* coincides with Aq^'"'^''. Then there are three cases. The first 
one is that X^jMi < and X^j(Aiiai)~^ = oo. This case is treated in Section 4. 
In this section, we are mainly studying the second case that ^ - /Xj = oo but 

^ 1 

Y < DO. (7.2) 

t^i ^^^^'^ 

The third case is that ^^Hi = oo and = which is treated in the 

next theorem. In this degenerated case, since there is a killing at 1 (i.e., ai > 0), 
the process is transient. Without using duality, by Corollary 7.3 below, we also 
obtain that Aq = 0. See the comments right after Corollary 7.3. 

Theorem 7.1. 

(1) First, let (7.2) hold. Define the dual rates [ai,bi) by (5.1) in the inverse 
way: 

bi = Oj+i, ai = bi, < N + 1, (7.3) 

and denote by Ai the eigenvalue defined in Section 6 for the dual process. 
Then we have Aq = Ai. 

(2) Next, let (7.2) fail and ^^^e l^i ~ Then Aq (as well as Aq'''^'') is equal 
to its dual Ao = Aq^'^^ = 0. 

Proof (a) By (7.3), (7.2) and (5.5), we have 

TV 



^fik <oo. (7.4) 



fc=0 



Clearly, the dual process with rates (a^, 6^) has the state space E = {i : ^ i < 
N + 1}. By exchanging (oj, 6j, Vi) and [ai,bi, Vi) in part (1) of Theorem 6.1, 



Ai = sup inf 

f, 0^i<N 



Oj+i + bi- Oi+iVi 

Vi-1 
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sup inf 



Vi-l 



- bi+iVi 



sup inf 



a,: 1 



the first assertion of Theorem 7.1 now follows from the variational formula given 
on the right-hand side of (9.2) in Section 9. 

(b) Similarly, replacing the use of Theorem 6.1 by Proposition 2.7 (1), we obtain 
the second assertion. In this CclSGj clS already mentioned at the beginning of 
Section 4, we have Aq = Aq'^'^-'. The fact that Aq = Aq^'^^ = comes from (5.5) 
and Theorem 3.1. □ 

By Theorem 7.1 (1), all the results obtained in Section 6 can be transformed 
into the present setup. For instance, by Corollary 6.4, we obtain the following 
result. 

Corollary 7.2. Under (7.2), we have (4(^)~^ ^ 5f ^ ^ Aq ^ S^^ , where 



1 



sup ^i[l,n]{u[n+l,N]+ 



\ /)"[1>^]1{7V«X)} 

V — ~L ' 



5l= sup {^/(p'^ + ^/(p^) ( V'i " V'l 



oi = sup 



l^wbN 

v[i + 1, iV] + (/^Af6Af)~^l|jv<oo} 

^[^^^] + (y"ArbAr)"^l{Af<oo} 



N-l 



^k+lfkAm 7— 31-{iV<oo}- 



k=l 



1^[l,N] + {fiNbN) ^l{Ar<oo} 



N-l 



/ > ^k+lVkAm H T~ ^{N<oo} 



j=i 



Proof. Starting from Corollary 6.4 with its notation, write everything we need in 
its dual. First by (5.4), we have 



fJ-n = ai^n+l, l^n = T-p'n+1, 0^n<N, fljsf = if N < OO. 



Here, recall that f„ = {jJinbn) ^ but z>„ = i^jlnC^n) ^ ■ Then the constant 6 defined 
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in (4.4) becomes S = sup^g^ z/[0,n — 1] iJ,[n,N]. Moreover, we have 

i-i ^ 

^1^^ = —fi[l,i], l^i<N + l, 



IJ,[rn, n\ 



j=0 



ai 



/ij = aiz>[r?7. + 1, n + 1], ^ m ^ n < iV, 



j=m 



IJ,[m,N] = aii)[m+l,N]+HNl[N^^} = ai 



i'[m+l,N] + - 



1 



{N<oo} 



fiNb 



N J 



m<N, 



N-l 



N-l 



ji'NON 



Inserting these quantities into (6.8) and (6.9), making a Uttle simphfication, and 
then ignoring the hat everywhere, we obtain Corollary 7.2. □ 

The next result is a criterion for the positivity of Aq, and is a particular case 
of Corollary 8.4 with B = L}{ij) in the next section. It is not deduced from 

the last section in terms of duality (Theorem 7.1) but conversely, it provides an 
improvement of Theorem 6.2 as shown by the proof of Corollary 6.6 below. 

Corollary 7.3 (Criterion and basic estimates). Without condition (7.2), we 
have K~^/4 ^ Aq ^ k""*^, where 



K 



inf 



E — ) +(E;^ 



Eft 



3=n 



Furthermore, we have 

5l A (5fl ^ AC ^ (l{s=oo} + (ai^)"^) (<5i A 5r) , 

where 

^ 1 1 

S =z] 1 7— l{Ar<oo}j 

^ Hitti unOn 

n ^ N m N ^ 

^R= sup Y^^^jJ^—] 



(7.5) 



% — i 7 — 



Km<iV+l fe=. 



Note that bi, = 6^^'^^ and Sji almost coincides with 5^^'^\ except for there is 
a shift of the state space. The second assertion of Corollary 7.3 means that Aq > 
iff the process goes to either or A/' + 1 exponentially fast. This is intuitively clear 
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by (7.1). Obviously, we have Xq = k -"^ = if Ylit^i — ^ ^'^d ^jifJ'jCij) "^ = 00 
since then Sl = Sr = 00. See also Corollary 8.6 below. 

Proof of Corollary 6.6. For given rates (oj, 6j) in the setup of Section 6, by (5.3), 
we have 

q ^ q ^ q-1 q q q-1 



i=p 



i=p—l 



J=P 



J=P 



j=p-l 



and 

AT 



N-1 



N-1 



tah ■ tiNbN ^0 . bo bo 



Regarding the process studied in Corollary 7.3 as a dual of the one given in the 
last section and then add a hat to each quantity of Corollary 7.3. It follows that 



inf 



inf 



inf 

0^n<m<JV+l 



i=l i=m 
n-1 ^ -1 /AT 

i=0 



+ ( E^^ 

=m 

N 



j=n 
In / m-1 



E^O + E^ 



j=1 

2-1 

E 

=n — 



^3 

■■ j=n-l 
-In /m-1 X -1 



j=n 



K 



Next, we have 

Vi=i MatOat / i=o 

Sl = Sr, and Sr = Sl- Since < 00, by Theorem 7.1 (1), we have Ai = Aq. 

Thus, Corollary 6.6 now follows from Corollary 7.3 immediately except for a slight 

change of the lower bound in the second assertion. For which, since Z < 00, the 
term "A(5i^" is not needed (cf. Proof of Corollary 8.4). □ 

Proof of Theorem 1.5. (a) Condition (1.3) implies that N = 00 and furthermore 
the uniqueness of the symmetric process on L^(/[x) by Proposition 1.3. Now, a* = 
Ai or Aq by [2; Theorem 5.3] or Proposition 1.2, respectively. 

(b) In the case that '^^fM = 00 and — process is zero- 
recurrent and so we have a* = 0. Noting that d^^-^\ k^^-^^\ and k^^-^) are 
all equal to infinity, the conclusions of the theorem become obvious. Hence, in 
what follows, we may assume that only one of Ylil^i Yli{fJ'ibi)~^ is equal to 
infinity. 

(c) Let bo = 0. Then the basic estimate follows from Corollary 7.3. 
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(d) We now prove the first two parts of the theorem under the assumption that 
bo = 0. In the case that Y^- < co but = oo, wc have k^'^-^'> = 
which gives us part (1) of the theorem. Next, if = oo but ^jX/iiO?:)"^ < oo, 
then for 5^ and 6r given in Corollary 7.3, we have 5^ = d'^^'^^ = oo and then 
K*^^'^^ < oo iff < oo. Clearly, (5r < oo iff S^^'^'^ < oo since < ^• 
This gives us part (2) of the theorem. 

(e) Finally, let 6o > 0. This is a dual case of that 6o = treated in (c) and (d). 
By exchanging the measures fi and v, we obtain the remaining conclusions of the 
theorem. 

Actually, part (1) of the theorem is a combination of Theorems 4.2 and 6.2, 
and part (2) is a combination of Theorems 3.1 and 7.1. □ 

We are now ready to prove an extension of Theorem 1.5. 

Theorem 7.4 (Criterion and basic estimates). Without condition (1.3), The- 
orem 1.5 remains true provided 

(1) the process in part (1) is replaced by the maximal one and a* is replaced 
by a™^^: the largest e such that 

\pij (t) -TTjl^ d e-'\ t ^ 0, (7.6) 

J 

for some L^ij^) -locally integrable function Ci depending on i only; and 

(2) the process in part (2) needs no change (i.e., the minimal one). 

Proof. Since Ai is equivalent to Aq'*'^^ (Theorem 6.2) and by duality, Ai = Aq^'^'' 

and Aq^'^^ = it is clear that is equivalent to Ag^'^''. Alternatively, 

one can use Corollary 7.3 to arrive at the same conclusion. Now, part (2) of 
the theorem follows by Proposition 1.2 for which we do not assume (1.3). As 
mentioned in the last proof, part (1) with the original a* also follows by [2] 
provided (1.3) holds. 

Even though in the previous study ([12], for instance), we consider only the 
ergodic processes under (1.2), but Ai can be actually identified with some expo- 
nentially ergodic convergence rate for more general ergodic processes (reversible 
Markov chains, in particular). First, the fact that the L^-exponential conver- 
gence rate is described by Ai does not require the regularity of the Dirichlet form 
(cf. proof of Proposition 1.1, for instance). Next, for a Markov process with state 
space {X, ^,7r) and transition probabilities {Pt{x, •)}, let £i be the largest e such 
that 

||Pt(x,-)-7r||var ^C(x)e-^*, xGX, t^O, (7.7) 

for some L^{X, 7r)-locally integrable function C{x) depending on x only. Then for 
a reversible process having density pt{x,y), we have Ai = ei provided 

P« (•>•) e AYc ^) for some s > 0, (7.8) 

and the set of bounded functions with compact support is dense in L'^{X, tt). The 
outline of the proof is as follows. 

(i) Prove that £i ^ Ai. 
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(ii) Show that \\{Pt — 7r)/||^ ^ C/ e~^^* for bounded / with compact support. 

(iii) Remove the constant C/ in the last line for fixed /. 

(iv) Extend / to L'^{X, it) and then claim that Ai ^ £i. 

By assumption, the last step is obvious. The detailed proof for the first three 
steps is given, respectively, in [12]: (8.6), the last formula in §8.3 replacing ei 
with El, and the proof of Lemma 8.12. Actually, this is a small correction to [12; 
Theorem 8.13 (4)] (i.e., replacing Si by £i) and its proof. It is known that £i > 
iff £i > (as well as £2 > used in the original proof of the cited theorem). Hence, 
the exponential ergodicity is kept but the rates £1 ^ £1 ^ £2 may be different. 
By the way, we mention that the change of topology is necessary in many cases. 
For instance, the pointwise convergence is natural in the discrete case but is not 
in the continuous case. In the ergodic situation, the total variation norm is good 
enough in general but it is meaningless in the non-ergodic case. 

Having this result at hand, part (1) of the theorem follows since we have a™*^ = 
£1 in the present context. □ 

We now introduce an interpretation, similar to Section 5, of the duality used 
in Theorem 7.1. For the ergodic process with Q-matrix, 



62 

-as. 

we have a simpler transformation matrix 



Q 





-60 


bo 









-oi - bi 


bi 


V 





0,2 


-02 - 62 








as 



ai,bi > 0, 



M = 







1^2 


^^3 









^3 










("3 


V 








)«3 



M- 












^ 


1 _ 




















1^2 ^ 


-1^2^ 




^ 











Then 



MQM-^ = 



















\ 


bo 


-ai - 


bo 


ai 










bi 




-a2 - 


bi 


a2 


62/ 










b2 




-as - 


















\ 


di 


-di - 


h 


h 










d2 




-0.2 - 




h 












as 




-ds - 



Wc obtain a process having an absorbing state at and being killed at the state 3. 
The original trivial eigenvalue with non-zero constant eigenfunction is transferred 
into the trivial one with eigenfunction IL{o}- Our dual matrix Q is now obtained 
by eliminating the first row and the first column from the matrix on the right- 
hand side. The elimination is to make the symmetrizability of Q and at the same 
time removes the trivial eigenvalue of the last matrix. Unlike the example given in 
Section 5 where the size of the state space stays the same: {0, 1, 2, 3} — ^ {1, 2, 3, 4} 
with a shift for the dual one, here the size of the state space is reduced by one: 
{0,1,2,3} ^{1,2,3}. 

We are now ready to examine some examples. 
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Examples 7.5. (1) Let N = 1. Then the Q-matrix is degenerated to be a single 

killing — c and so Xo = c = . 
(2) Let A^ = 2. Then 

Ao = (ai + 02 + 6i + 62 - V (ai - 02 + 61 - 62)^ + 40261 )• 



The next two examples are taken from Chen, Zhang and Zhao (2003, Exam- 
ples 2.2 and 2.3) 

Examples 7.6. (1) Let iV = 2, oi = 02 = 1, 61 = 2, and 62 = 3. Then Aq = 2, 
and by Corollary 7.2, we have 

^1 ^ Ao' =0.5^(5i, 

where 

^ 4W3 _ = ^ = 0.46, « L23. 

10 '15 Si 

Next, K ^ Aq ^ ^ 4k with k = 3/7. Obviously, [Si, Si) C (k, 4k). 
(2) Let AT = 2, 61 = 1, 62 = 2, 

ai = ^^, £g[0,x/2), 

and 02 = 1. Then Ao = 2 — e, and we have 

^Ao' = (2-£)-^ ^^1, 



where 



4 + V2 + (2 + ^/2)e-£2 1 {l + e){V2-e 
di = — — ' 



^1 



8 + 2e-3e2 Aq 8 + 2e-3£2 ' 

8 + 6£-£2 1 £2 



Hence, 



16 + 4£-6£2 Ao 2(8 + 2£-3£2)' 



j, ^ ^_ 2(4- V2 )(! + .) 

5l 8 + 6£-£2 



Next, K ^ Aq ^ 4k with 



1 . [ 1 e 

«; = mm 



2 



Ao 1^8 + 2£-3£2' 8-4£2+£3 

Even though it is not so obvious now but we do have [Si, Si) C (k, 4k). 
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Examples 7.7. Because of Theorem 7.1, we can now transfer [10; Examples 
9.27] into the present context, sec Table 7.1, by using (5.1) and (5.7). Here, for 
the sixth example, we need a restriction: 1/k < b/a ^ k/{k — 1)^ (fc ^ 2). 

Table 7.1 Exact Aq for nine examples 



tti {i ^ 1) 


bi (i ^ 1) 




Vi {i ^ 1) 


a 


b(a<b) 


\ V * / 




7i(i-l)+7o 
7o > 0, 7i ^ 


/3ii(^i>7i) 


Pi - 71 


7ii + 70 


i-l + Po 
Po>0 


2(i+l)+/3o 


2 


(z + l)(i + /3o) 
i[2{i+l)+Po] 


i 


2i+4+\/2 


3 


i + 1 r 2(z + x/2) 1 

2/+4+v/lL /(' + 2v^-l)J 


a 
i 


6 


, \/a?+Aab—a 
^ 2 


\/a2+4a6 + a 
26 i 


a 


(i A A;)6 


[Vbk-V^]^ 




i + 1 




2 




ai > 




1 
4 


2z + l 
2{i + 1) 


2 + (-1)^-1 
(^>2) 
7-V33 
2 


2[2+(-l)i 


6- V33 


V33 + (-l)^ 
8 



We now go to the second part of this section. Consider the birth-death pro- 
cesses with a more general state space E = {i : —M — 1 < i < + 1}, M, N ^ oo 
and with Dirichlct boundaries at — M — 1 if M > — oo and at + 1 if A?" < oo. 
Its Q-matrix now is = 5j > 0, = Oj > 0, and Qij = if |i — j| > 1 for 

i,j G E. Fix a reference point 6 E E. Define 



Correspondingly, 
D{f) = 



ae-iae-2 ■ ■ ■ oe+n+i 



bebe-i 
1 



' be+n 



<n^-2, 
1 



bebe-i aebe' 
be+ibe+2 • • • ^e+n-i 
aeo-e+i ■ ■ ■ ae+n ' 



-M - 1 - 
2^n<N+l-e 



(7.9) 



E 

M-l<i^e 



i^idiifi - Si-if + 



E 



l^ihifi+l — fi 



(7.10) 



f e J^, f-M-i = if M < CX3 and /at+i = if A?" < cxd. 
Let us begin with a particular application of Corollary 8.4 to 



lHm)- 
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Corollary 7.8 (Criterion and basic estimates). Let Aq be defined by (7.1) 
with the present state space E. Then we have /A ^ Aq ^ where 



K = inf 



i=-M 



j=m 



(7.11) 



By the way, wc extend Corollary 6.6 to the present general state space. 

Corollary 7.9 (Criterion and basic estimates). Let YlieE A^j < oo and define 
Ai as in (6.1). Then we have k~^/4 ^ Ai ^ k~^, where 



inf 



m \ -1 / N 



i=-M 



n — 1 ^ 



(7.12) 



Proof. When M < oo, the corollary is simply a modification of Corollary 6.6 by 
shifting the left end-point of the state space from to —M. Thus, when M = oo, 
we can choose a sequence {Mp}"^^ such that Mp t oo as p t and then the 
assertion holds if M is replaced by Mp for each p. In which case, the corresponding 

Ai is denoted by A^^^"^ for a moment. Because ^^^^ l^i following the proof 

above (4.2), it follows that 

Ai= inf {£>(/) : m(/^) = 1, //(/) = 0, /»= /(ivm)An for some m,ne E, m< n]. 

Hence, we have A^^"-* | Ai as p\ co. Similarly, replacing M by Mp, we have the 
notation The proof will be done once we show that 



(k(^'')) asp too. 



Obviously, we have 



(/«(*^-)) \ aspt and {k^^^'^) ^ ^ . 

To prove the required assertion, let £ > 0. Then by definition of k there exist 
mo, no G -E, mo < no such that 



mo \ -1 
i=-M 



N 

+ ( 

i=no 



-In / riQ-l ^ 

y 



^ +e. 



Next, since Ylif^i < fo^" fixed mo, no and large enough Mp {—Mp < mo), we 
have 



I- ^ i=—M„ 



Y f^i) + E ^ 



N 



i=no 



no-1 



"U -L -I 



j=mo 



^ + 2£. 
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Combining these facts with the definition of k^^p\ we obtain 



-1 / N X -It / n-1 

inf I / — 



-Mp^m<n<N+l 
^ K"^ +2£. 



7 



Since e is arbitrary, we have proved that (^k^^p^^ ^ ^ k ^ as p ^ oo. □ 

For the remainder of this section, we study a splitting technique. It provides 
a different tool to study the problem having bilateral Dirichlet boundaries. This 
approach is especially meaningful if the duality discussed in Section 5 does not 
work, such as in studying the processes on the whole Z or the Poincare-type 
inequalities given in the next section. We remark that Corollaries 7.8 and 7.9 use 
slightly the splitting idea only (cf. Proof (b) of Theorem 8.2 below). The idea is 
splitting the state space into two parts and then estimating the first (non-trivial) 
eigenvalue in terms of the local ones. Wc have used this technique several times 
before: Chen and Wang (1998) with Dirichlet boundary for the unbounded region, 
Chen, Zhang and Zhao (2003), as well as Mao and Xia (2009), with Neumann 
boundary. The first and the third papers work on a very general setup. Here, we 
follow the second one with some addition. 

To state our result, we need to construct two birth-death processes on the left- 
and the right-hand sides, respectively, for a given birth-death process with rates 
(aj,6j) and state space E. Fix a constant 7 > 1. 

(L) The process on the left-hand side has state space = {i : —M — 1 < 
i ^ 9}, reflects at 9. Its transition structure is the same as the original 
one except ae is replaced by 70^). 

(R) The process on the right-hand side has state space E^~^ = {i : 9 ^ i < 
-1- 1}, reflects at 9. Its transition structure is again the same as the 
original one except be is replaced by 7(7 — l)~^be. 

For the process on the right-hand side, the state is a Neumann boundary. At 
AT -I- 1, it is a Dirichlet boundary if A?" < 00. For this process, the first eigenvalue, 
denoted by Aq"'"'''', has already been studied in Sections 2 and 3. With a change 
of the order of the state space, it follows that the process on the left-hand side 
has the same boundary condition, denote by Aq its first eigenvalue. Note that 
ignoring a finite number of the states does not change the positivity of Aq , 
in the qualitative case, we simply denote them by Aq^\ respectively. In general, 

according to and/or Ei=-M(^i«i)"^ Eile and/or ^L-M^i 

being finite or not, there are eight cases for the processes on Z. For instance, if 
Yl!i=6^t^i^i)~^ = 00, then Aq"'"'' = by Theorem 3.1. Since in this section, we are 
working on bilateral Dirichlet boundaries, it is natural to assume that Aq'^'' > 0. 

The other cases may be treated in a parallel way. For instance, when Aq ^ = 0, 
it is more natural to consider the process on [— M, + 1) with reflecting at some 
finite —M and then pass to the limit as M ^ 00 (cf. the proof of Corollary 7.8). 
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In this case, the eigenfunction should be strictly decreasing once Aq > 0. Hence, 
there is no reason to use the splitting technique. Note that the explicit criterion 
for A^^^ > is given by Theorem 3.1. We can now state the main result of the 
second part of this section as follows. 

Theorem 7.10. 

(1) In general, the Dirichlet eigenvalue Aq of the birth-death process on 
E = {i : -M -1 <i <N + 1} satisfies 

V A^'^) = Ao ^ sup sup (A^-^ A A^'^) , (7.13) 



inf inf (An 

eeE 7>i ^ ^ 



sup 

M-l^es^N+l 7>1 



where on the right-hand, when 9 = — M — 1, define Aq ''^ = oo, and Xq^''^ 
to be the first eigenvalue of the original process (independent ofj) reflected 
at —M if M < oo; when = N + 1, define Xq^''^ = oo, and Aq"''*' to be 
the one reflected at N if N < oo. 
(2) The second equality in (7.13) replacing sup_jvf_i^e^Ar+i sup^gg also 
holds provided Aq^-* > 0, and moreover, 
e N 

Hi = oo if M = oo and /Xj = oo if N = oo. (7.14) 

i=-M i=e 

Theorem 7.10 was proved in Chen, Zhang and Zhao (2003) for the half-space 
(i.e., one of M and A'' is finite), under the hypotheses that J^iif^i^i)'^ < oo and 
Hi <oo which is essentially the case of having a finite state space. 
To prove Theorem 7.10, we need some preparation. First, we couple these two 
processes on a common state space E = {i : — \ < i < N + 2} . Next, separate 
the two processes by shifting the state space E^~^ by one to the right: 1 + E^~^ . 
Denote by (ai,6i) the rates of the connected process. For this, we need to build 
a bridge for the processes on the two sides by adding two more rates 6^ = 7 — 1 
and ae+i = 1. The construction here will become clear once we have a deeper 
understanding about the eigenfunction and it will be explained in Part H of the 
proof of the theorem. Roughly speaking, there are two possible shapes of the 
eigenfunction, the construction enables us to transform one of them to the other 
so that the splitting with Neumann boundaries becomes practical. For which, one 
needs the parameter 7 as shown in Lemma 7.12 below. In detail, we now have 

bi, -M-Ki^e-1, 



1, 



-M - 1 < i ^ 



1, 



i = e + i, 

e + 2^i<N + 2, 



7-1, 



I = 



7 
6,-- 



1 



^ + 1, 



6* + 2 ^ ?; < iV + 2. 

Applying (7.9) to the present setup and removing the factor 65/(1 — 7)"^ (which 
simplifies the notation but does not change the ratio D(/)//x(/^)), we obtain 



li-i 



He 



1 



7 



-M- 1 <i^e-l, 
7-1 



7 



Ate, 



(7.15) 



1) 



+ 2 ^ i < iV + 2. 



SPEED OF STABILITY FOR BIRTH-DEATH PROCESSES 



83 



Then 

_-7-l 

fiitti = fiitti, i^9, fJ-ebe = /iq, Hibi = fii-ibi-i, i^9 + l. (7.16) 

7 

The next two results are basic in using the sphtting technique. 

Lemma 7.11. Given f on E, define f on E as follows: fi = fi for i ^ 9 and 
fi = fi-i for 9 + 1. Then we have fi{p) = ,x{f^) and D{f) = D{f). 

Proof. Clearly, we have fe = fe+i- Then 

-M-Ki^e-1 6l+2^i<Ar+2 
W)= E m^{f^-h-lY+ U^Hf^+l-f^Y (7.17) 

= Y l^i<^iUi - fi-lf + Yj - /j-l)^ 

-M-l<i^6> 6»+l^i<Ar+2 

= D{f) (by (7.10)). □ 

Lemma 7.12. For a given birth-death process with state space E and rates {ai, bi), 
if its eigenfunction g of X satisfies ge-i < ge > 99+i (resp. ge-i > ge < ge+i) for 
some 9 & E (of course, g_j^_i =0 if M < oo, and 5jv+i — ^ if ^ = oo), let 

7 = 1 + ^4^^^^^>1, (7.18) 
aeige - 9e-i) 

and let gi = gi for i ^ 9, gi = gi-i for i ^ 9 + 1. Then for the {ai,bi) -process, 
g is the eigenfunction of X = X having the property go+i = ge. Furtherm,ore, 
g\(^-M-i,e] is the eigenfunction of X of the process on the left-hand side reflecting 
at 9, and similarly 3|[5i_|_i_7v+i) is the eigenfunction of the process on the right-hand 
side reflecting at 9 



Proof. By the construction of {ai,bi) and g, we have 
ng(i) = 



Qg{i) = -Xgi = -Xgi, i < 6* - 1, 

Qg{i - 1) = -A^i-i = -Xgi, i^ 9 + 2. 

Next, by (7.15), we have 

Ug{9) = beige+i - ge) + ae{ge-i - ge) = ae{ge-i - ge) = iae{9e-i - ge), 
ng{9 + 1) = be+i{ge+2 - ge+i) + ae+i{ge - ge+i) 

= be+iige+i - ge) 

7 

^ be{ge+i- ge)- 



7-1 



In the first formula, the term containing bg vanishes. This is the reason why 
we can regard ^ as a reflecting boundary for the process on the left-hand side. 
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Similarly, one can regard ^ + 1 as the one for the process on the right-hand side 
in view of the second formula. By (7.18), the right-hand sides are the same which 
is equal to 



, beige -ge+i) 



aeige-i-ge) = ae{ge-i- ge) + he{ge- ge+i) = -Xge = -^ge- 



aeige - ge-i) 
We have thus proved the lemma. □ 

Proof of Theorem 7.10. Part I. In this part, we prove Theorem 7.10 (1) with the 

first "=" replaced by The proof of this part is relatively easier. Let / G 

/ 7^ 0. Define / as in Lemma 7.11. For fixed 6 G E and 7 > 1, noting that 

re-labeling the state space does not change Aq"*"''*^, by (7.17), we have 



Hence by Lemma 7.11, 

^ Ar-^ A A 



Dif) Dif) ^ , ,e+n 



Making the supremum with respect to 7 and 9, it follows that 
— — ^ sup sup (Aq A Ag^ ■ 

fJ'KJ ) -M-l<e<N+l 7>1 

At the boundaries, say 9 = —M — 1 for instance, by (7.10) and the convention, 
we have 

Dif)> J2 ^^Mf.+i-f^f>C''Kf) = [Ar'^'AA^-]/x(f). 

-M-l<j<JV+l 

Therefore, we indeed have 



^> sup sup(A^^AA^'^). 



Making infimum with respect to /, we obtain 

Ao = inf ^HX^ sup sup (Aq"'^ A Aq"^'^). 

This proves the (second) inequality in (7.13). 

To prove the upper estimate, fix G -E and 7 > 1 again. As we have seen 
from the last part of proof (g) of Theorem 2.4 and Proposition 2.5, if we let 
^s-i-,7,n jjgj-^Q^g ^YiQ local eigenvalue with Neumann boundary at 9 and Dirichlet 
boundary at n + 1, then Aq"'"'''''"' J, Aq"'"'''' as n f 00. Thus, for each £ > 0, we 
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have Aq^''*''"' < Aq^'''' + e for large enough n. By Proposition 2.2, we can assume 
that the corresponding eigenfunction of Aq'*''''''"' satisfies g^^'"'^ = 1 and 

g(+,n) ^ Q i > n{> 6). Similarly, we have A^"'^'"" < A^"'^ + e for 

small enough — m, and moreover, the eigenfunction of Aq '^'"^ satisfies 

g\ = 1 and g\ = for all i < —m{< 0). Let / be defined as above, 
connecting g(~'"^^ and g^~^'"'\ Then / has a finite support, fe = /e+i = 1, and 
moreover by (7.10), 

By Lemmas 7.11 and 7.12, this gives us 

since e is arbitrary. Furthermore, we have 

Ao ^ inf inf (An"''^ V An+''^) 

as required. □ 

The proof of the equalities in Theorem 7.10 is much harder. For which, we need 
once again a deeper understanding of the eigenfunction of Aq. To have a concrete 
impression, we mention that the eigenfunction in Examples 7.6 (2) is go = gs = 0, 
gi = {1 + £)g2- Thus, when e = 0, we have gi = g2- Besides, it is rather easy to 
see the shape of eigenfunction g of the examples given in Table 7.1 since Vi < 1 
iff gi+i < gi for all i. 

Definition 7.13. 

(1) A function f is said to be unimodal if there exists a finite k such that fi 

is strictly increasing for i ^ k and strictly decreasing for i ^ k. 

(2) A function f is said to be a simple echelon if there exists a k such that 
fk = fk+i, fi is strictly increasing for i ^ k and strictly decreasing for 
i^ k + 1. 

Proposition 7.14. Let g be a positive eigenfunction of X > for a birth-death 
process. Then g is strictly monotone, or unimodal, or a simple echelon. 

Proof, (a) Let gk ^ gk+i for some k. We prove that g is strictly decreasing for 
i ^ k + 1. To do so, note that 



bk+i{gk+2 — gk+i) — —Xgk+i — ak+i{gk — gk+i) ^ —Xgu+i < 0. 
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Thus, we have gk+2 < 9k+i- Assume that gn < gn-i for some n ^ k + 2. Then 
the eigenequation shows that 

bn{gn+l - gn) = ->^gn - an{gn-l - gn) < ->^gn < 0. 

By induction, this gives us gn+i < gn for all n ^ A; + 1. 

(b) By symmetry, we can handle with the case that g^ ^ gu+i for some k. One 
starts at 

ak{9k-i - gk) = ->^gk - bk{gk+i - gk) ^ ->^gk < 0. 

We obtain gk-i < gk and then gn-i < gn for all n ^ A; by induction. 

(c) By (a) and (b), it follows that there is no local convex part of g. Otherwise, 
there is a k such that either gk-i > Qk < fffc+i or gk-i > gk = gk+i < gk+2 which 
contradict what we proved in (a) and (b). 

(d) We claim that for every k, say A; = for simplicity, the two cases "g-i ^ 
go" and "(70 ^ gi" cannot happen at the same time. Otherwise, there are four 
situations: 

9-i=go = gi, g-i>go<gi, g-i>go = gi, and g-i=go<gi- 

The first one cannot happen, otherwise we have gi = 0. By (c), the second case 
is impossible. The last two cases are also impossible by (b) and (a), respectively. 

(e) Having these preparations at hand, we are ready to prove the main assertion 
of the proposition. Clearly, we need only to consider the case that g is not strictly 
monotone. Choose a starting point, say for instance. By (d), we have only one 
possibility: either g^i ^ go or go ^ gi- Without loss of generality, assume that 
9o ^ gi- If 9o = gi, then by (a) and (b), 5 is a simple echelon. If go < gi, then 
on the one hand, by (b), gi is strictly increasing for all i ^ 1, and on the other 
hand, we can find a A; ^ 1 such that gi < g2 < ■ ■ ■ < gk ^ 9k+i since g is not 
strictly monotone by assumption. Applying (a) again, it follows that g is strictly 
decreasing for alH ^ A; + 1. Hence, g is either unimodal or a simple echelon. □ 

Proposition 7.15. For the birth-death process on Z, the following assertions 
hold. 

(1) The eigenfunction g of X satisfies the following successive formulas: 



gk+l — 9k + 



fJ-kbk 
1 

gk-1 = gk-\ 

fJ'kO'k 



IJ^eaeige - ge-i) - 



l^eaeige-i - ge) - A jijgj 



i=k 



k^9, 

(7.19) 

k<e. 



(2) If X = 0, then the non-trivial eigenfunction g with ge = 1 for some 9 E 7. 
is given by 

n-l j 

j=e k=9 ^ 

0-1 9-1 , 

g^ =1 - {1 - g,_,)J2 n "<^- 

, . , 1 (^k 
J=n k=j+l 
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In this case, the function g is either the constant function 1 or strictly 
monotone on Z. 
(3) // A > and 

e oo 

^ IJ'i = ^l^i = ^, (7-20) 

i= — oo i=9 

then the non-trivial eigenfunction g of X cannot be monotone. 

Proof, (a) Part (1) of the proposition follows from the eigenequation. 

(b) When A = 0, with Ui := ^j+i — gi{i E Z), the eigenequation = OjUj-i 
gives us 

Uj = (1 - 9e-i) n IT' ^' % = (1 - 9e-i) IT ~' ^ 
k=0 k=j+i ""^ 

It follows that cither gi = 1 oi g is strictly monotone on Z. Now, part (2) of the 
proposition follows by making a summation of j from ^ to n — 1 and from n to 
0—1, respectively. 

(c) Without loss of generality, assume that ge = I for some 6 e Z. Suppose 
that g is non-decreasing, then by the first equation in part (1), we would have 

oo > ^ Hugh ^ 2^ Hk ^ oo as n ^ oo. 

k=0 k=e 

Otherwise, if g is non-increasing, then by the second equation in part (1), we 
would have 

, s 9-1 0-1 

H0ae[g0-i - 90) ^ 
oo > ^ 2^ ^ 2^ A^fc ^ oo as ra ^ -oo. 

k=n k=n 

We have thus proved part (3) of the proposition. □ 

We remark that Proposition 7.15 (2) is different from Proposition 2.2 where 
the eigenfunction of A = must be a constant. Here is a simple example with 
9 = 0: Qi = bi = \i\ if i ^ and ao = bo = 1, then corresponding to A = 0, we 
have a family of linear eigenfunctions {g^^^ = 1 + — i : i E (normalized 
at 0) with one-parameter 7. 

Proposition 7.16. Let (7.14) hold and g be a non-zero eigenfunction of Xq > 0. 

Then g is either positive or negative on E. 

Proof. If one of M or iV is finite, then the conclusion follows from Proposi- 
tion 2.2 (1). From now on in the proof, assume that M = N = 00. 

(a) If the conclusion of the proposition does not hold, then there is a A; (say) 
such that gk ^ and either gk-i > or gk+i > 0. By symmetry, assume that 
9k+i > 0. 
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(b) We now prove that > for alH ^ A; + 1. Given m, n G Z with m ^ n, 
denote by Aq"*'"^' the first eigenvalue of the process restricted on the state space 
{i : m ^ i ^ n} with Dirichlet boundaries at m — 1 and n + 1 in the sense similar 
to (7.1). If the assertion does not hold, then there is a fco ^ > ^ + 1 such that 
Qko ^ 0. Now, let g satisfy gt = 0, gi = gi ioi i = k + 1, . . . ,ko - 1, gko = s for 
some £ > 0, 5j = for i ^ fco + 1. Note that 

( - n 5) (fc + 1) = bk+i {h+i - ak+2) + flfc+i [h+i - 9k) 



Because of Aq > and following proof (b) of Proposition 2.1, we can choose a 
suitable e > such that 



It follows that Aq ' " < Ao- However, it is obvious that Aq ^ Aq ' " and so 
we get a contradiction. We have thus proved that gi > for all i ^ fc" + 1. 

(c) By (7.14) and proof (c) of Proposition 7.15, g cannot be non-decreasing since 
Ao > 0. Hence, there is a ^ fc + 2 such that gk+i < gk+2 < ■ ■ ■ < ge ^ 9e+i- In 
the case that ge > ge+i, by introducing an additional point but keeping the same 
Ao as shown in Lemma 7.12, one can reduce to the case that gg = gg+i- Hence, 
one can split the original process into two as in (L) and (R). Now, starting from 
at which gg > 0, look at the process on the left-hand side in the inverse way, one 
finds the point fc < at which g^ ^ 0. Applying proof (b) above to this process, 
one may get a contradiction. It follows that g > on (—00,^] D (— oo,fc]. 

Therefore, we should have g > on Z. □ 

Proof of Theorem 7.10. Part E. We now prove the equalities in (7.13). By as- 
sumption Aq"^^ > and the second inequality in (7.13), it follows that Ao > 0. If 
one of M and N is finite, then the non-trivial eigenfunction g must be positive 
by Proposition 2.2 (1). In this case, it is helpful to include the boundary into the 
domain of g for understanding its shape. Then by Proposition 7.14, there are only 
two possibilities: 

(i) g is unimodal; 

(ii) (/ is a simple echelon. 

Next, if M = = 00, then by Proposition 7.16, we have g > 0. Moreover, by 
Proposition 7.15, g cannot be monotone. Hence, by Proposition 7.14, g has again 
one of shapes (i) and (ii) as above. 

We now prove the equalities in (7.13) only in the case that M = N = 00. The 
proof for the other case is simpler. 

(a) Case (ii). We use the operator // defined in Section 2: 



= bk+l{9k+l — 9k+2) + 0,k+l{9k+l — 9k) + o,k+igk 
= Ao5fc+i + ak+igk 

^ ^0 9k+l, 




i=k+l 



i=k+l 
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For each / satisfying: /j = fi ioi i ^ 9 and /j = for z ^ ^ + 1 for some / on 
E, by (7.15) and (7.16), we have 



N+l 



3=1 



"7 — — 7 — iJ'0+ife + /^fc-i/fe-i 



k=e+2 



N 



7 



i-1 



k=e+i 

N 



= ^ y -^Yukfk--^ y — 



AT 



+ 1 ^ i < iV + 2. 



(7.21) 



Similarly, we have 



j=-M k=j 



M-l<i< 



(7.22) 



Because gg = ge+i , we can regard as a Neumann boundary of the original process 
restricted on the left-hand side and at the same time, regard ^ as a Neumann 
boundary of the original process restricted on the right-hand side. Because Aq^"* > 
0, by Proposition 2.5 (2), we have g±^ = 0. Hence, by (2.11), (7.21), and (7.22), 
we obtain 



Hi [9) 



^0 V 7/ 9i-i ^ijbj 



E^. 



e + l^i<N + 2, 



^--(,)= 1 E 



-M-l<i<e. 



By Proposition 2.2 (2), we have 



N 



sup 



e^i<N+i 9i l^jbj Ao 



sup -—^ > ^ — . 
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Therefore, by the second inequahty in (7.13) and Theorem 2.4 (3), it follows that 



Ao ^ sup sup 

e'&E 7>i 



' + ,7 



^ sup 

7>1 

^ sup 

7>1 



( inf ^r'"(5)-')Af inf li'^'^i^)-') 



sup inf if-^'^^iq] ^ 

^>l 6l+l<i<iV+2 ^ ' 

sup <^ — + 1 - - sup > — 



Ao. 



We have thus proved in Case (ii) the second equality in (7.13). 

To prove the first equality in (7.13), noting the inequality was proved in Part 
I, we have dually 



An ^ inf inf 
7>i 



A^ V Ao 



'+,7' 



^ inf 

7>1 



Ao"'^ V A^+'^ 



^ inf 

7>1 



( sup lfr^''{^)-'\y{ sup 

inf sup Il\ '^(5) ^ 

7>1 -M-l<i<6» 



7> 



sup 

7>1 L '^O 



7 



sup 



* 1 

E — 



However, there is a problem in the second line of the proof. To apply Theo- 
rem 2.4 (3), one requires that either g £ L'^il^-) or g is local. Hence, an additional 
work is required. Anyhow, the conclusion holds whenever both M and A'^ are 
finite. We will come back to the proof in proof (c) below. 

(b) Case (i). By Lemma 7.12, this case can be reduced to Case (ii). Actually, 
the proof becomes easier now. With 7 given by (7.18), we have 

11'.^'^ {g) = 1-, e+i^i<N + 2, 

^0 

Hence the second equality in (7.13) holds. Moreover, the first equality in (7.13) 
also holds whenever both M and N are finite. 
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(c) To complete the proof for the first equahty in (7.13), we need to overcome 
the unbounded problem. For this, choose Mp,Np f oo as p — oo. Denote by 
^^-,7,p^ ^^e+,J,p Aq^\ respectively, the quantities Aq"'"*, Aq^''^, and Aq when 
M and N are replaced by Mp and Np. Note that for a finite state space, we 
certainly have Aq'^ > 0, its cigenfunction is positive (by Proposition 2.2 (1)) and 
has properties (i) and (ii) mentioned in the above proof (by Proposition 7.14). 
Clearly, for each fixed and 7, we have 

Ag^''*''^ 4, Ag^'''', Ao^"* 4, Ao as p — >■ 00. 

Thus, as proved in (a) and (b), whether we are in Case (i) or (ii), we have for 
each p, 



ee[-Mp,Np] 7>i 



^ inf 

Oe[-Mp,Np] 7> 

> inf inf [An 

eeE 7>i '- " 



inf [An 

7>i '- ^ 



,7,PvA^+'^'P] 



'"^ V A^+'^] 



'^VA^'^]. 



Therefore, by the first inequality in (7.13) proved in Part I, it follows that 

Ao ^ inf inf [A^"'^ V Ag+'^l ^ A^,^' ^ Aq as ^ 00. 
6eE 7>i 



We have thus completed the proof of the theorem. □ 

Here are remarks about the assumption made in part (2) of Theorem 7.10. 
Similar to the upper estimate, we do have 



sup sup 

ee[-Mp,Np] 7>i 



A 



The problem is that Aq '^'^ | Ag '''' as p — > cx) goes to the opposite direction and 
the approximating sequences {Mp} and {Np} depend on 9 and 7. Hence, the proof 
for the upper estimate does not work for the lower one. Next, to prove the second 
equality in (7.13), it seems more natural to assume that Ag"''*^ AAq"'"''*^ > for some 
6 and 7, that is, Aq A Ag"'"'' > 0, rather than Aq V Aq"*"^ > as we made. However, 
if one of them is zero, say Ag ■* = 0, then as mentioned before Theorem 7.10, we 
have a single Dirichlet boundary but not the bilateral Dirichlet ones, and the 
variational formula takes a different form (i.e., the second inequality in (7.13) at 
the boundaries). Condition (7.14) is due to the same reason. In particular, when 
M = —1, for instance, if X^j/^i < 00 and X^j(/Ujaj)~^ = 00, then Ag^'' = by 
Theorem 3.1, and we go back to the case studied in Section 4. In which case, the 
cigenfunction of Ag is strictly increasing. 

To conclude this section, we introduce a complement result to [12; Proposi- 
tion 5.13] about the principal eigenvalue for general Markov chains. 
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Proposition 7.17. Let {qij : i,j G E) be symmetric with respect to {jii) on a 
countable set E, not necessarily conservative (or having killings): 

di --Qi-^ qij ^ 0, Qi := -Qii e [0, oo). 

Define 

D{f) = IY1 - fif + E ^'^difl 

i,jeE ieE 

and 

Ao = inf {D{f) : f has a finite support and /^(/^) = l}- 
Then we have infjgs Qi ^ Xo- 

Proof. Without loss of generality, assume that E = Z_|_ = {0, 1, . . . }. Fix k e E 
and take / = t{k}- Then n{f^) = Hk and 

Dif)= E ^^^<1^3U3-h?+Y.^'^'^^f^ 
i,j: i<j ieE 

= E l^kQkjifk - fjf + E IJ-i^likifi - fkf + E i^i'^ifi 
3>k i<k ieE 

= E t^kQkj + E t^ilik + IJ'kdk- 

By the symmetry of iiiQij, we get 

j>k i<k j¥=k ' 

It follows that 

Ao ^ D{f)/f,{f) = qk. 
The assertion now follows since k & E \s arbitrary. □ 

8. Criteria for Poinc are- type inequalities 

As in [9] for the crgodic case having N < oo or (1.2), the results studied in 
Sections 2, 3 and 7 can be extended to a more general setup, so called Poincare- 
type inequalities. In this way, one obtains various types of stability, not only the 
L^-exponential one studied in the other sections of the paper. Here we consider 
only the criteria and the basic estimates for the inequalities. In other words, we 
extend Theorems 3.1, 4.2, and 6.2 to the general setup with some improvement. 
At the same time, we introduce a criterion for the processes studied in Section 
7 in this setup. To do so, we need a class of normed linear spaces (B, || • ||b, /u) 
consisting of real Borel measurable functions on a measurable space (X, ^,11). 
We now modify the hypotheses on the normed linear spaces given in [12; Chapter 
7] as follows. 

(HI) In the case that fj,{X) = 00, Ik G B for all compact K. Otherwise, 1 G B. 
(H2) If /i G B and |/| ^ h, then / € B. 
(H3) ||/||b = supgg^ l/lyd/x, 
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where to be specified case by case, is a class of nonnegative =^-measurable 
functions. A typical example is ^ = {1} and then B = L^{ji). Throughout this 
section, we assume (H1)-(H3) for (B, || • ||b, //) without mentioning again. 
Before moving further, let us mention the following result. 

Remark 8.1. Without using (1.2), the results in [9] remain true under the condi- 
tion Hi < oo replacing the original process with the maximal one if necessary. 

The key reason is that without condition (1.2), the same conclusion holds in 
Section 4 on which the cited paper is based on. 

In this section, our state space \s E = {i : —M — l<i<A^ + l} (M, N ^ oo) 
as in the second part of Section 7. The next result is the main one in this section; 
it has several corollaries as we have seen in the last section. Note that the factor 
4 in (8.2) below is universal, independent of B. 

Theorem 8.2. Consider the minimal birth-death process with Dirichlet boun- 
daries at —M — 1 if M < CO and at N + 1 if N < oo. Assume that ^ contains 
a locally positive element. Then the optimal constant Ag, in the Poincare-type 
inequality 

\\f^\\^^AgD{f), fe^^'^D), (8.1) 



satisfies 



Bn ^ ^ 4Sb 



(8.2) 



where the isoperimetric constant Bm can be expressed as follows: 



Bts 



inf 



i=-M 



1 



+ 



y— 



-In 



111 



m,n\ I 



(8.3) 



In particular, when B = U'/'^{ijl)(jp ^ 2) (then (8.1) is called the Sobolev-type 
inequality), we have 



b: 



inf 



m ^ 

E — 



y— 

t—'n ' 



2/p 

Eh) ■ (8-4) 

j=m 



Proof, (a) First consider the transient case, in particular when one of M or is 
finite. We use the proof of [11; Corollary 4.1] or [12; Corollary 7.5] with a slight 
modification. In proof (a) there, it was shown that one can replace ^ 1" by 

= 1" in computing the capacity Cap{K) for compact K. Without loss of 
generality, assume that / ^ 0. Otherwise, replace / with |/|. In the proof just 
mentioned, the condition "X^^ < oo" was used so that 1 G S'{D). We cannot 
use this assumption now, but for a given nonnegative / G ^'"^(L') n ^c{E), 
where ^c{E) is the set of continuous functions with compact support, we can 
simply choose a nonnegative smooth h G '^c{E) such that /i|supp(/) = 1- Then 
h G ^"'"(L*), / A G ^'"•"(L'), and so one can use f A h G ^™'"(D) instead of 
/ A l to arrive at the same conclusion D{f) ^ D{f A h) as in the original proof 
(a). 
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The first step in the original proof (b) shows that one can replace a finite num- 
ber of disjointed finite intervals {Ki} by the connected one [minUjiTj, maxUj-fCj]. 
This part of the proof needs no change. 

Note that in the original proof, the state space is {1,2,...} with Dirichlet 
boundary at 0. The main body of the original proof (b) is to find a minimizer 
(actually unique) / G '^dE) for D(f) having the properties /o = and fix = 1- 
Replacing N with q for the consistence with the notation used here and let K = 
{m, m + 1, . . . ,n} (1 ^ m ^ n, here m and n are exchanged from the original 
proof). Now, within the class of /: /q = 0, /|x = 1 and supp(/) = {1, . . . ,q} 
(n ^ q < N + 1), the minimal solution is 

-(/)-(Ei)""H|:^)"" (-) 

To handle with the general state space, one needs to move the original left-end 
point 1 of the state space to somewhere, say p > —M — 1. In detail, replace the 
condition m ^ 1 used in defining the compact set K by m > — M — 1. At the 
same time, replace {!,... , q} hy {p,p + 1, .. . ,q} with — M — 1 < p ^ m for the 
supp(/). Then the last formula reads as follows: 

^ i=p ' ^ t=n ' 

In the original proof, the ergodic condition and (1.2) are mainly used here to 
remove the second term on the right-hand side. We now keep it. Since the right- 
hand side is increasing in p and decreasing in by making the infimum with 
respect to /, it follows that 

Cap(iC) := inf {£>(/) : / G n ^,{E) and /|x ^ 1} 

/ m ^ ^ -1 . N ^ .-\ 

The assertion of the theorem now follows by using 

n III-RtIIb ||l[m,n] ||b 

ijB := sup — — = SUp 



K Cap(ir) _M-i<m^n<Ar+i Cap([m,n]) 

and applying [11; Theorem 1.1] or [12; Theorem 7.2]. The last result is an exten- 
sion of Fukushima and Uemura (2003, Theorem 3.1). 

(b) Next, consider the recurrent case: both X]i<e(/^i'^j)~^ ^"^^ J2i>eil^i^i)~^ 
are diverged. Here is actually a direct proof of the lower estimate in (8.2). Without 
loss of generality, assume that the reference point ^ = 0. Fix m' ^ m ^ and 
n' ^ n ^ 0. Based on the knowledge about the eigenfunction given in the last 
section, and similar to proof (b) of Theorem 3.1, define 



fi={ 



fc=iVn ^'^ " 
iA(— m) ^ 

7 Y] ^{i^-m'}, ^^0, 

k=—m' 
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where 

7 := 7(m',m,n,n') = V — 7-/ . 

fe=n k=—m' 

Here, 7 is chosen to make / be a constant on [— m,n]. By (7.10), we have 

— m n 

i=—m' i=n 
—m 1 n' -j 



-In 



Moreover, 



Hence, 



p2|| > || ^|2 



-m,n\ 



An > ^^(^y ^ II l[-m,n] I 



/ 1 \^ 



k=—m' 



Prom this, we obtain the lower estimate in (8.2). Since contains a locally 
positive element, we have ||l[-m,n]l|B > for large enough m and n. Letting 
m',n' 00, by the recurrent assumption, it follows that = 00. Besides, it is 

obvious that Ba = 00 in this case and so the first and then the second assertion 
of the theorem becomes trivial in the recurrent case. □ 

Proof (b) above indicates an easy improvement of the lower bound of ^b- Use 
the same / as above, and define 



(m,m' ,n,n') 



—m 1 / — m 

E E — 

,±tl,f^kak/ ,^^,f^kak} 



-m,—m- 



+ 



Then a simple computation shows that 
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Hence 



D{f) 



^WM-mM+h^""'"" '''''' \ { E ^) 



k=—m 



j=n 



Noting that the right-hand side is increasing in m' and n', and making a change 
of the variable —m—^m,we obtain 



^ sup 



Denote by Cb the right-hand side. Then the conclusion of Theorem 8.2 can be 
restated as Bb ^ Cb ^ ^ 4Bb. Certainly, this remark is meaningful in other 
cases but we will not mention again. 

The next result is an easier consequence of Theorem 8.2. 

Corollary 8.3. Everything in the premise is the same as in Theorem 8.2. Then 
(1) we have ^ -Bl A Br, where 

1 ^ 1 

Bl = sup N ||l[n,JV+l)||B, -Bit = sup \ r-||l(_M-l,nll|B- 

neE ._ J.. fJ'iO.i neE - ^ fJ'iOi 

The equality sign holds once 

N 



+ 



1 



i=-M 



l^iO-i UNbN 



-1 



{N<oo) = OO. 



(2) Next, we have Bn ^ (Bl A Br)1^s=oo} + S ^B, where 



B = sup 

m,n&E: m^n L 



i=-M 



k=n 



Proof. Clearly, by (8.3), we have 



i=-M 



meE V ' II n - 



i=-M 



and so B^ ^ Bl. The equality sign holds once Yl!i=eil''i^i) "'^=00. Similarly, 
we have i?B =^ Br. The equality sign holds once 'Ylfi=-Mil^i^i)~ ~ Hence, 
B-a ^ Bl /\ Br and the equality sign holds once S = 00. 
Next, when 5 < 00, we have 



[m,N+l)\ 



Bi^ ^ S inf 



i=-M 



N 
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We have thus proved the corollary. □ 

Of course, one can decompose the constant B in Corollary 8.3 (2). For instance, 
for fixed mo, we have 

mo -, 



^ 1 



mo,n\ 



5 > ( V I sup 

The last factor is close to Bji when mo is negative enough. However, when 
mo — )• —M, the first term tends to zero since S < oo, unless M < oo. This 
indicates that bounding B^ by B^ and Bji is rather rough, especially in the case 
that E = Z (cf. Example 8.9 below). This is a particularly different point of the 
processes on the whole Z or on the half space Z+ as shown by Corollary 8.4 below. 

We now specify Theorem 8.2 and Corollary 8.3 to the half space: either M or 
N is finite. This corresponds to the processes studied in the first part of Section 
7 (see Corollary 7.3). 

Corollary 8.4. In Theorem 8.2, let M = —1. Then we have 

-1 / N ^ . -1-, 



Btb 



inf 



111 



Furthermore, we have 

BlABr^Bm^ (1{S=oo} + {aiS)-^) {Bl A Br), 

where 



(8.6) 



1 

Bl= sup nil 



N 



n,N+l)\ 



B 



R 



s 



E- . 



+ 



-1 



{iV<oo}- 



(8.7) 



Proof. The first assertion follows from Theorem 8.2 with M = — 1 and an ex- 
change of m and n again. The second one follows from Corollary 8.3 except the 
last estimate. When S = oo, we have B^ = Bl. While when 5 < oo, we have 



B^^ ^ S inf 



^ aiS inf 

l^m<Ar+l 



K=m 



n,m\ 



aiSB;^\ 



Therefore, 

Bm ^ Bl1{s=oo} + {aiSy^BR ^ {l{s=oo} + (ai-^)"') {Bl A Br) 
as required. □ 

When one of M or is finite and its Dirichlet boundary is replaced by the 
Neumann one, the solution becomes simpler. The next result corresponds to the 
processes studied in Sections 2 and 3. 
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Theorem 8.5. Let M = be the Neumann boundary and assume that ^ contains 
a locally positive element. Then the isoperimetric constant := sup^ || l^f Us/Cap (if) 
can be expressed as 



^ 1 

Bm= sup V ||l[0,n]||B- (8.8) 

0^n<A'+l fMOi 
^ i=n 

In particular, for the Sobolev-type inequality, we have 

N I / ^ \2/p 
0^n<N+lf^lJ'ibi\'^^ J 

Proof. The proof is nearly the same as that of Theorem 8.2 except one point. 
In proof (b) of [11; Corollary 4.1] or [12; Corollary 7.5], to find a minimizer / 
for D(f), since the constraint /o = and /n = 1, / cannot be a constant on 
{0, 1, . . . ,n}. Now, without the constraint /o = 0, the minimizer should satisfy 
fj = 1 for all j : ^ j ^ n. Thus, instead of (8.5), the minimal solution becomes 



Then the necessary change of the proof of Theorem 8.2 after (8.5) should be 
clear. □ 

Applying Theorem 8.5 to IB = L^(//), we return to Theorem 3.1. Actually, in 
parallel to [9], one may extend the results in Sections 2 and 3, Theorem 3.1 in 
particular, to the present setup of normed linear spaces and then deduce Theo- 
rem 8.5. The next result is obvious, it says that for a null-recurrent process, the 
LP {p ^ l)-Sobolev inequality is still not weak enough. 

Corollary 8.6. Consider a birth-death process on Z_|_. IfYli^if^i ~ °° ^'^^ 
= DO, then B^'^^ = B^'^"* = oo for all p ^ 2. 

Remark 8.7. We now compare (8.6) and (8.8) in the particular case that J2i 
= oo. Then the constant Bm given in (8.6) becomes 

B^-^^ = sup ||l[l,m]l|B- 

Rewrite the constant B-g. given in (8.8) as 

N , X / N 



^ (si "'^°^'"") V (r,S.r£i "''°'"'"")- 
By (H3), we have 

l|l[l,n]l|B ^ l|l[0,n]l|B ^ ||1L{o}]||b + ||1[i,„]||b. 

Next, by (HI), we have ||1{o}||b < oo. It follows that B^'^'' < oo iff B^'^'' < oo. 

We conclude this section by a simple example to show the role of the Poincare- 
type inequalities. 
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Example 8.8. Consider a birth-death process on Z_|_ with fii = {i + 1)'^ (7 > 1) 
and hi = 1. Then ai = i'^{i + 1)"''' and 



n>0 



n -1 2/p 00 ^ 

^ i=0 



j=n 



U + 1)^ 



P>2. 



Hence, Bp^'^^ < 00 iff 



P>2{1 + 



7-1 



However, S'^^'^^ = ^2^'^^ = 00 for all -f > 1. 

Example 8.9. Let E = Z, bi = 1, Hi = e*^ and B = Then for the 

quantities given in Corollary 8.3, we have Bl = Br = 00 but B^ < 00. 

Proof. Obviously, Bl = Br = 00. To show that Bb < 00, since 

xVy^x + y^ 2{x V y), 

it suffices to prove that 



sup 



m -|X /CO ^ \ -I n 



k=n 3=m 



jJLj < 00. 



By symmetry, without loss of generality, it is enough to show that 



sup 



I— —00 ' j=—m 



or 



/ —m \ m 



m>0 



00. 



Equivalently, 



sup 



Uiai 

i=—oo ' j =—m 



00 ^ \ m 

Ei E". 



^0 \ Mi^*! 



< 00. 



3=0 



The assertion now follows by using Conte's inequality: 



X 



and Gautschi's estimate: 
1 r 



{xP + 2f/P-x Ke""" / e-y^dy^Cj 



x^ + 



Cr 



1/p 



(8.10) 



£e^^^g(l-e-^^), .^0 



— a; 



a; ^ 0, 



Cp = r(l + l/p)^/^^ '\ p>l; C2=7r/4. 



Alternatively, one may check directly that the function under supremum on the 
left-hand side of (8.10) is decreasing in m 1) and then (8.10) follows easily. □ 
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9. General killing 

In Sections 4 and 7, we have studied the special case having a killing at 1 only. 
We now study the process with general killing, as described by (2.1) with state 
space shifted byl: E = {i : 1 ^ i < N + 1}. We use the same symmetric measure 
as in Section 4. 

The next preliminary result is quite useful. To which it is more convenient to 
use oi + ci and + cjv for the killing rates at boundaries 1 and (if < oo), 
respectively, rather than ci and cjy used in Proposition 2.1. Note that the killing 
rates in the next proposition are allowed to be zero identically. 

Proposition 9.1. Let (a^) and {hi) he positive hut ai ^ 0, hj\[ ^ if N < oo, 
and let (cj) he nonnegative on E. Define Aq = Xo{ai,bi,Ci) as follows: 

Ao = inf {£>(/) :M(r) = l,/eJr}, 

where 

-^(/) = IJ-Aifi+i - fif + fJ-iaifi + '^l^iCifi, /jv+i = if N < oo. 

ieE ieE 

Write Aq = Ao(oi,&i,0) for simplicity. Then we have 

(1) Ao(aj,5j,c-) ^ \o{ai,hi,Ci) if c[ ^ Cj for alii £ E. 

(2) Xo{ai,bi, Cj + c) = Xo{ai,bi,Ci) + c for constant c ^ 0. 

(3) Aq + sup^g^ Ci ^ Aq ^ Aq + inf jgfi Cj and the equalities hold if is a 
constant on E. 

Proof. Since a change of {cj}^i makes no influence to {/ijj^i, part (1) is simply 
a comparison of the Dirichlct forms on the same space i^^(/x) with common core 
J(f . Similarly, one can prove the other assertions. □ 

Note that Proposition 9.1 makes a comparison for the killing rates only. Actu- 
ally, a more general comparison is available in view of [3; Theorem 3.1]. Next, if 
(1.3) holds, then by Proposition 1.3 and the remark below (4.3), the Dirichlet is 
unique, and so the condition f £ J(f can be ignored in defining Aq. 

It is worthy to mention that the principal eigenvalue Aq studied here can be 
extended to a more general class of Schrodinger operators. That is, we may 
replace the nonnegative potential (cj) with the one bounded below by a constant: 
infj Ci ^ — M > — oo. Then we have q + M ^ for all i and 

\o{ai,hi,Ci) = Xo{ai,bi,Ci + M) - M ^ -M. 

Having Proposition 9.1 at hand, all the examples for Aq given in Sections 3, 
5 and 7, can be translated into the case of Aq with constant killing rate. For 
instance, we have the following example which already shows the complexity of 
the problem studied in this section. 
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Example 9.2. Let Oj = a > for i ^ 2, bi = b > and Cj = c ^ for 1. 

(1) // ai = a, or ai = but still a ^ b, then Aq = {s/a — V^)^ + c. 

(2) // ai = and a > b, then Aq = c. 

Proof. By Proposition 9.1, we need only to study Aq. In the last case, since the 
process is ergodic, we have Aq = 0. Next, we have Aq = {^/a — Vb) according 
different cases by 

(i) Example 5.3 if ai = a and a ^ b, 

(ii) Example 7.7 if ai = a and a ^ b, 

(iii) Example 3.4 if ai = and a ^b. □ 

From now on, we return to a convention made in Section 2, the rates ai and 
bN are combined into Ci and cat if AT < oo. Thus, in Theorem 7.1 for instance, we 

have ai = and ci > 0, and moreover, 6jv = and c^r > if A'^ < oo. In general, 
we assume that Cj ^ 0. Otherwise, we will return to what we treated in Sections 
2 and 3. Define the operator R: 



Ri{v) = 0,(1 - v^\) + hi{l - Vi) + Cj, 

for -y in the set y = > : 1 i < N}. 
N = oo, define 



i E E, vo = oo, vn = if N < oo. 
Next, define Y = Y ii N < oo. When 



1^ = 1^ < Vi : Vi > for i < m, Vi = ioi i ^ m> 

m=l ^ ^ 

\^ <v : Vi > on E, the function /: /i = 1, /» = JJ t'fe (^ ^ 2) is in L^{iJ,) 
^ k=l 

and satisfies O/// ^ rj on E for some constant rj^ . (9.1) 

For V e y with finite support, R,{v) is also well defined by setting 1/0 = oo. 

Theorem 9.3. Assume that Ci ^ 0. For Aq defined by (2.2) with state space 
E = {i : 1 ^ i < N + 1}, the following variational formulas hold: 

ini_ sup Ri{v) = Aq = sup inf Ri{v). (9.2) 

ver ieE ver '^-^ 



Proof, (a) First, we study the lower estimate. In the case that ^^^e l^k < oo, as 
a particular consequence of [5; Theorem 1.1], we have 



Ao ^ sup inf (i) = sup inf 

g>o g g>o ieE 



ai 1 



9i-i 



+ bAl 



9i+i 



+ Ci 



(9.3) 



where g^ := and gN+i = if < oo. The proof remains true when XlfceB l^k = 
oo, simply using E^ = {1,2,... ,m} (m < iV + 1) instead of the original one. 
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Actually, the conclusion holds in a very general setup (cf. Shiozawa and Takeda 
(2005) and its extension to the unbounded test functions by Zhang (2007)). 

Suppose that Aq > for a moment. Then by Proposition 2.1 (with a shift by 1 
of the state space), the eigenfunction g of Aq is positive. It follows that the first 
equality sign in (9.3) can be attained and so does the last equality in (9.2) with 
Vi = gi^i/gi > 0(1 ^ i < N). Next, if Aq = 0, then N = oo since and bi arc 
positive for i : 2 ^ i < iV, and Ci ^ (in the case of Theorem 7.1, we have ci > 
and also cat > if TV < oo). By setting = 1 for i G E, we get 



inf 

ieE 



ai(l — ) + bi{l - Vi) + Ci 

Vi—l 



^ inf Ci ^ 0. (9.4) 



Hence, the last term of (9.2) is nonnegative. Therefore, the last equality in (9.2) 
is trivial if Aq = 0, in view of (9.3). 

(b) Next, we study the upper estimate. We consider only the case that N = oo. 
Otherwise, the proof is easier. Given v e y,let j = 7(1;) = supj^j^^^ Ri{v) and 
as in the definition of y, set 

/o = 0, /i = 1, /i = n vk, i ^ 2. (9.5) 

k=l 

First, suppose that supp (v) = {1, 2, . . . , m— 1} for a finite m. Then supp (/) = 
{1,2,... , m\ and 

—$1/ 

{i) = Ri{v) ^-i, i = l,2,...,m. 



/ 

Hence, 



fe=i fc=i 

m+l TO m 



l^l^kfl >^l^kfk{-^f){k) 

k=l 
m+l 

= ^2 f^kakfk-lifk-l — fk) — ^2 f^kakfk{fk-l — fk) + ^ fJ-kCkfk 
k=2 k=l k=l 

m 

= ^ fJ-k [akifk-l — fkT + Cfc/I] + fJ-m+iam+lfmifm — fm+l) 

k=l 
m+l 

= ^l^k [ak{fk-i - fkf + Ckfl] ■ 

k=l 

We have not only 7^0 (actually 7 > when m is large enough since Cj ^ 0) but 
also 

^0 ^ ^ ^ liv) (9.6) 
for all V & y with finite support. 
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(c) Next, we are going to prove (9.6) in the case that v & Y with > for 
six i ^ \. In this case, the positivity condition of v is not enough for the first 
equahty in (9.2), as mentioned in Section 2 (above the proofs of Theorem 2.4 and 
Proposition 2.5). See also the specific situation given in the proof of Example 9.17 
below. This explains why two additional conditions are included in the second 
union of the definition of Y . The condition "/ € L^(/i)" is essential but not the 
one ^^^f / f ^ T]" since the eigenfunction g of Aq satisfies ^^^g/g = — Aq". To prove 
(9.6), without loss of generality, assume that 7 = 7(f) < 00. Otherwise, (9.6) is 
trivial. Clearly, 7 ^ Ri{v) = 6i(l — vi) + ci > —00. Note that by assumptions, 
the function / possesses the following properties: 

(i) / > on E. 

(ii) / G L'^ifJ') and then Ptf € L'^{^J.), where = {Pij(t)) is the minimal 
semigroup determined by the Dirichlet form. 

(iii) \nf(i)\ = I X;^. Qijfjl ^ max{|r/|, \j\}fi for all i G E. 

Here, property (iii) comes from 

-vf ^ -nf ^ 7/. 

Since ipij{t)) satisfies the forward Kolmogorov equation: 

Pijit) = Sij + / y'pifc(s)5feids 
•^0 k 

and (i), it follows that 

Ptfii) = /i + XI / '^Pik{s)qkjfjds. 
3 k 

By (ii), Ptf{i) < 00 and is continuous in t. Because of this and (iii), the order of 
the last two sums and also the integration are exchangeable. This leads to 

Ptf{i)^ fi-1 I y^Pik{s)fkds = h-^ j PJ{i)ds, ieE, (9.7) 
Jo J. Jo 

since by assumption fif ^ —7/. Therefore, we obtain 

< D{f) = lim !(/, / - PJ) ^ Urn ^ j\f, PJ)ds = j{v)\\ff < 00. 

Here, the first limit is due to (ii) and the first equality in (1.10), the last inequality 
comes from (i) and (9.7). We have thus proved that not only 7 > but also 
/ G i^(D) and so we have returned to (9.6). In other words, (9.6) holds for all 
V . By making infimum with respect to G 1^, we obtain 

Ao ^ inf_ sup Ri (y) . 
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(d) To prove the equality sign in the last formula holds, in view of proof (b), 
we have actually proved that for every finite m, 



A^™) := inf{D(/) : /o = 0, /, = for alH ^ m + 1, ||/|| = 1} 



^ inf sup Ri{v), 



(9.8) 



where _ 

Yrn = {vi ■■ Vi > ior i < m, Vm = 0}. 

Actually, there is a v G 1^ such that Ri{v) = Aq"*^ > for alH (1 ^ z ^ m) since 
m < oo and then the equality sign in (9.8) holds. Therefore, the first equality in 
(9.2) holds since Aq 4- Ao as tti 'I" oo. n 

We now begin to study the estimate of Aq. First, by Proposition 7.17, we have 
a simple upper bound: 

Ao ^ inf (ai + bi + Ci). 

Hence, Aq = whenever lim„^oo((in + bn + c„) = 0. The next result provides us 
a finer upper bound. It is motivated from Theorem 3.1. 

Proposition 9.4. Zei Cj = Cj — inf j Cj . Then 

-1 r / m -, \ —1 m 



Ao ^ inf Ci + ini ( } Hi] inf 
ieE eeE\^^ J E3m^e 



^ inf Ci + inf ( > 
ieE eeE \ ^ 



i=\ 



I 



I I V ^ ^ _L III/ 



l^iCi 



i=l 



i=l 



(9.9) 
(9.10) 



Proof. By Proposition 9.1, it is enough to consider the case that Cj = Cj, i.e., 
inf,- d =0. Fix £ < m and define 



Then 



'Pi = 



(£,m) 



m ^ 



E. 



Hence, 



i=l i=e+l i=l 

m £ m m 



k=l 



i=l i=l+l 



i=l 



t X -1 
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Because (^(^'"^) G JT, it follows that 
An ^ inf inf "^^^'^ 



t N -1 



^ inf ( > iia\ inf 



inf 

E3m^e 



m ^ \ —1 m 



^ inf ( y^/^ 
eeE V 



i=l 



(=1 



i=l 



□ 



As an immediate consequence of (9.10), we obtain the following result. 
Corollary 9.5. IfYl'iZil^i ~ °° ^^'^ 



then 



Mm flrnbm ( Mi ) = 0, 



1=1 



Ao ^ inf Cj + lim } UiCi / } Ui ^ inf + lim c. 



m— >00 ,_i 



i=l 



Proof. Without loss of generality, assume that Ci = Ci. 
By assumptions, it follows that 



lim ( ^//i 

i=l 



1^1 



i=l 



^ ^^^^ i m ^0 + lim ( y^/x, 



i=l 



i=l 



The first inequality now follows from (9.10). 

To prove the second inequality, let 7 = lim„_).oo c„ G [0, 00]. Then for every 
£ > 0, we have supj^.>„ ^ 7 + e for large enough n. Hence, 



X /^«*^' = X '^^^^ + X ^^"^^ ^ X + (7 + e) X ^ 

i=n+l i=l i=n+l 



n. 



We have thus obtained 

/ e .-1 e 

— ( X '"m X ^ /™ ( X 



Y HiCi + {j + e) ^ 

i=l i=n+l 



= 7 + £ 
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as required. □ 

To study the lower estimate of Aq, we observe that not every positive sequence 
{vi) is useful for the lower estimate given in (9.2) since one may have inf j Ri{v) < 0. 
In order for inf, ^ 0, it is necessary that 



1 

0<Vi^ — 
h 



Ci + ai + hi 



Vi-l 



From this, we obtain the following necessary condition: 



<Vi^Xi- 



Xi-l - 



Vi-l 



Xi-2 



yi-2 



Xs - 



y2 

X2 

Xl 



where 



Xi 



Cj + Oj + hi 



yi = -r- 
hi 



However, the condition is clearly not practical. Because of this reason, we are 

now going to introduce an alternative variational formula for the lower estimates. 

For a given sequence (rj), define an operator 11^ = of "double sum" on 

the set of positive functions (fi) as follows: 



i-l 



i-1 



k=l j=l j=l k=j 

Write //(/) = II^{f). For a fixed sequence (cj), let Ci = Ci — infj Cj and define 

^ = {f>0:fi<fi + Ilfif) for allE3i^2}. (9.11) 

Clearly, if Ci > 0, then every positive constant function belongs to Otherwise, 
every / > with /j < /i for all E 3 i ^ 2 belongs to ^. 

Theorem 9.6. Let , (cj) and ^ he defined as ahove. Next, for each fixed 
f E define 



E3i^2 IIi{f) 



N = oo, 



mf ^^^^ A^^lv N<oo, 



CiriJ) 



inf 

E3i^2,ci<ri 



Ci + 



iv - Ci)fi 



V, 



[0,^]. 



{E3i^2: Ci<vi} ^ 
{E 3i^2: Ci <r]} = 



(9.12) 



(9.13) 
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Then we have 

Ao ^ inf c, + C(??,/) and vXivJ), f e ^, r, e [0,^]. (9.14) 

Moreover, for fixed f, C(?7, /) is increasing in r) and furthermore, 

Ao = inf Q+ supC(?,/). (9.15) 

Remark 9.7. To indicate the dependence on (q), rewrite CiVi f) '^^ C(ci, "?/,/)• 
Similarly, we have ^{ci, /). Then for each / G ^ and constant ^ 0, we have a 
shift property as follows: 

?(ci + 7,/) = 7 + C(Q,/), C{ci + 7,V + 7J) = l + C{ci,vJ)- (9.16) 

Hence, the use of inii^E Ci in Theorem 9.6 is not essential but only for simplifying 
the computations. The same property holds for (9.10) but not for (9.9). 

As will be illustrated later by Examples 9.17 and 9.19, it is not unusual that 
> Ao for some / G In that case, we certainly have ^/ > C(0' /)• This means 

that may not be a lower bound of Aq and so the use of (^{rj, f) in Theorem 9.6 

is necessary. 

Proof of Theorem 9.6. By Proposition 9.1, for simplicity, we assume that q = Cj. 

(a) First, we prove "Aq ^" in (9.14). Fix / G ^. Then ^ = ^ 0. Without 
loss of generality, assume that (^ ^) ?? > 0. Otherwise, the assertion is trivial. Let 

hi = f,+Iir''if), ieE, rje{0,^]. 

Since by (9.12), 

fi-fi + Ii:if)^Vmf)>0 

foTE3i^2 and /ii = /i > 0, we have h > 0. Next, define Vi = 
{vo ■= oo and vn = if N < oo). Then for i :2 ^i < N, since 

hi - = iir\f) - iit+^u) = -V E(^ - c,)n,f,, 



we have 



ai(l - v^\) +bi{l- Vi) 



^ [ai{hi - hi-i) + bi{hi - /li+i)] 



1 

hi 



i-l 



XI ~ Cj)/xj/j + X(^ - 



lii-ibi-i ^ fiA .^^ 



iv - Ci)fi 
hi 
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This also holds when i = 1 (noting that ai =0): 



6i(l - Vi) = — (/ti - h2) = 7-^ — =r]-ci. 

Ill III 



11 N < oo, then at i = A^, by assumption 

TV 



TV 



j = l I 3 = 1 



we get 

aiv(l - v'^_^ + hj^{\ - vn) = - 



TV-l 



hNli-N-l^N-l 



(rj - CN)fN 



h 



N 



Combining these facts together, we arrive at 



Ri{v) = Cj + Oj (l - v-\) +bi{l- Vi) ^ Ci + 



jri - Ci)fi 
hi 



i e E. (9.17) 



We now show that the right-hand side of (9.17) is nonnegative for all i and so 
we have ruled out the useless case that infj Ri{v) < 0. Since /i > 0, the assertion 
is equivalent to 

Cihi ^ (cj - ry)/i, i G E, 



or 



This is trivial if CiIIi{f) ^ fi (in particular if i = 1) since fi — fi + Hiif) ^ for 
all i? 9 ^ ^ 2 and / G Otherwise, by the definition of ^ and r], we have 

fi-fi + nt{f) > mif) ^ Vlliif) > ri [Iliif) - fi/ci] , E3i^2. (9.18) 

We have thus proved the required assertion. 
By Theorem 9.3 and (9.17), we obtain 



Ao ^ sup inf 



Ci + 



= sup < ?7 A inf 



iv - Ci)fi 

f,+iir\f). 

iv - Ci)fi 



(9.19) 



Ci + 



h+iinf) 



Here, the last line is due to the fact that IlKf) = 0. 

(b) To prove the first assertion of the theorem, we show that for each i: 2 ^ 
i G E, 

{V - Ci)fi 



Ci + 



^77 iff Ci^rj. 
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Clearly, the inequality is equivalent to 

{ri-c,)f-^{rj-c,)[h + Iir\f)]. 
The required assertion then follows since by (9.18), we already have 

f.^f.+Iir^if). 

As a consequence of the assertion, we have rj ^ C(^)/)- Now, from (9.19), it 
follows that 

Ao ^ sup r/ A C(??, /) = sup Civ, /)• 

This gives us the first assertion of the theorem. 

(c) To prove the monotonicity of ({r], /) in rj, let ??i < r/2 ^ ^. If {E 3 i ^ 2 : 
Ci < 112} = 0, then {E 3 i 2 : Ci < rji} = and so 

({v2,f) = v2>vi = CiviJ)- 

If 3 i ^ 2 : Q < 771} 7^ 0, since {E 3 i ^ 2 : Ci < iji} C {E 3 i ^ 2 : Ci < 772}, 
we need only to show that 

JUlZ^^Jni^f^ on {E3i^2:c,<r„}^t 

Actually, this is enough even ii {E 3 i ^ 2 : Ci < rji} = $ in view of (b). Now, 
the required conclusion is trivial on the set {E 9 i ^ 2 : 771 ^ Cj < 772}. Hence, it 
suffices to show that 

^ on \E 3 1^2: Ci <m\. 

A simple computation shows that this is equivalent to 

h+IIt{f)^CiIh{f), 

which holds on {E 3 i ^ 2 : Ci < 771} in view of (9.12) and ^ > 771. 

(d) To prove (9.15), it suffices to show that the equality in (9.19) holds for 77 = ^. 
Noting that the right-hand side of (9.19) is nonnegative, without loss of generality, 
we may assume that Aq > 0. Then, by Proposition 2.1, the eigenfunction g > 
of Ao satisfies 

k 

i^'khigk - gk+i) = Xl^-^o ~ (^j)H9j, ke E, gN+i = if at < 00. 

Hence, 

N 

gi-gi = II^°-%g), i e E, ^(Aq - c,)/x,-5,- = if AT < 00 
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and furthermore, g E It follows that 



iiiig) 



= Ao if AT < OO, 



Ci + T^^zy— = Ao, ^ G E. 



(Aq - Ci)gi 
g^^Iir'^^ig) 

Therefore, .^g = Aq, and furthermore, the equality sign in (9.19) is attained at 

(/,r/) = (5,Ao). □ 

We now make a rough comparison of Theorems 9.6 and 9.3 for the lower esti- 
mate. See also the comment below the proof of Corollary 9.9. 

Remark 9.8. For a given positive sequence {vi) such that inf jg^; Ri{v) := 7„ ^ 0, 
corresponding to the sequence (/j) and^f defined by (9.5) and (9.12), respectively, 

we have ^ 7„. 

Proof. Prom the assumption 

Riiv) = Ci + ai(l - v^\) + bi{l - Vi) ^ 7^ =: 7, i e E, 
it follows that 

1 ^ 

fk - fk+l ^ ~ Cj)/Xj/j, 

and then 

i-l k 

- ^ E 7^ - = i e E. 

k=i ^ j=i 

To prove our assertion, without loss of generality, assume that 7 > 0. Then it is 
clear not only that / G =^ but also ^/ ^ 7. □ 

As a complement to Remark 9.8, it would be nice if we could show that 

(0 (i')Ji ^ ^j^g set {E 3 i ^ 2 : Ci < Ef}. 

This holds obviously on the subset {7^ ^ Cj < ^/}, but is not clear on the subset 
{E3i^2:c,< 7^}. 

The next result is a particular application of Theorem 9.6. It is a complement of 
Corollary 9.5. The combination of Proposition 9.4 and Corollary 9.5 with Corol- 
lary 9.9 below indicates that when Ao(aj,6j,0) = 0, the condition lim„^oo c„ > 
is crucial for Ao(aj, 6j, Cj) > 0. This is more or less clear in terms of the Peynman- 
Kac formula: 

where {i^t^jt^o is the minimal semigroup generalized by the operator with rates 
{ai,bi,Ci), {Xt}o<^t<T is the minimal process with rates {ai,bi), and r is the life 
time of {Xt}. Note that Ao(ai,fei,0) > 0, and hence, XQ{ai,bi,Ci) > if the 
uniqueness condition (1.2) fails. Otherwise, r = 00. 
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Corollary 9.9. Let e G (0, 1). Define 

\ — £ + ciZi + exi 



inf 



Zi + evi 



N = oo, 



. 1-e + ciZi+exi » ci + e 2^ '^j^j 
mt /\ , N < oo, 



(9.20) 



inf 



Ci + 



1 + ciZi + exi - ^s{zi + £yi)_ 



(9.21) 



where 

Xi= ^ Cj yi= ^ fiji/[j,i-l], Zi = u[l,i-1], 

and viij] = Ylif^k^jil^kh)'^- Then we have Xq ^ infj^B Q + sup£g(o,i) Ce- The 
same conclusion holds if in (9.21) is replaced by rj E [0,^£]. In particular, if 
lim , Cn > 0, then An > 0. 

Proof, (a) The main assertion of the corohary is an apphcation of Theorem 9.6 
to the specific / G /i = 1, /j = e G (0, 1) (£^ 3 i ^ 2), for which we have 

//[(/) = 0, _L.+,^r,,,Y: ^ 

Then (9.20) and (9.21) follows from (9.12) and (9.13), respectively. 

We now prove the particular assertion for which N = oo. 

(b) If (1.2) does not hold, then Ao(oj, 6j, 0) > by Theorem 3.1, and so Aq > 
by part (3) of Proposition 9.1. Similarly, if inf^ Cj > 0, then we have again Ao > 0. 
Thus, without loss of generality, assume that 

irif Cj = and (1.2) holds. 



(c) With the test function / given in (a), by (9.12), we have 

l-e + IIf{f) 



Ce = inf 

i>2 



iiiU) 



By assumption, there exist 7 > and m ^ 2 such that Cj > 7 for all i ^ m. 
Certainly, we have 



Ce^ inf 



1-e + IIfif) 



/\inf 
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For i > m, we have 



By assumption (1.2), the right-hand side goes to £7 > as i — > oo. It fohows 
that infj>^ //?(/) ///j(/) > 0, and furthermore, there exists rj G (0,7) such that 

(d) Noting that the set {i ^ 2 : Cj < r/} C {i : 2 ^ z < m} is finite, by (9.13), 
we have 



^ mm > 0. 



Now, by Theorem 9.6 or proof (a) above, we conclude that Aq > 0. □ 
From proof (c) above, we have seen that when N = 00, 

Ce>0 iff inf //f(l)//7i(l) > for all m ^ 2. (9.22) 

i>m 



Note that 

i>m //j(l) 



j=i j=i 



and the right-hand side is positive iff 

m f m 

liSl Y.f'-^^->/ E/^i >0. (9.23) 

j=l I j=l 

Thus, Corohary 9.9 is quaUtatively consistent with Corollary 9.5. 

In view of Remark 9.8, it is not obvious that Theorem 9.6 improves Theo- 
rem 9.3. An easier way to see the improvement is as follows. Recall that the last 
assertion of Corollary 9.9 is deduced in terms of the test function / used in its 
proof (a). For which, the corresponding sequence (wj) is vi = 1/2 and Vi = 1 for 
all i ^ 2. Inserting this into R{v), we get 

inf Ri{v) = (ci -I- 61/2) A (c2 - 02) A inf Cj. 

Thus, for infill i?i(T;) > 0, it is necessary that inf^^s q > 0, which is clearly much 
stronger than the last condition lim ^^^ c„ > used in Corollary 9.9. 
As Proposition 9.4, the next result is also motivated from Theorem 3.1. 
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Corollary 9.10. An explicit lower estimate can be obtained by Theorem 9.6 using 
the specific test function /^"^^ ; 

/ m , X 1/2 

where m may be optimized over {m G £" : m ^ 2} (or over E if ci > 0). 

We now show that some special killing (or Schrodinger) case can be regarded 
as a perturbation of the one without killing. To do so, fix constants /3, 7 > 0, and 
define 

Oj = bi-i, 2 ^i < N + 1, bi = Oj+i, 1 ^i < N, 
di = /3, bN = ^ ii N < 00; 
Oj = ttj+i, 0^i<N, bi = bi, l^i<N + l. 

bo = P, aN=l UN <oo. (9.24) 

Note that (dj, 6j) and (dj, hi) are dual each other in the sense of Section 5 but they 
are clearly different from (aj,6j). Recall that ai = and = by convention. 
Next, let (cj) satisfy 

' Oj+i - Oj - 6j + 6j_i, 2^i < N, 
Ci^ < a2-bi+ P, i = l, (9.25) 

j — ttN + bN-i, i = N<oo. 

Note that the right-hand side of (9.25) can be negative. Conversely, for given 
rates (di, bi) , the inverse transform is as follows: 

ttj = Sj_i, 2 ^i < N + 1, bi = dj+i, l^i < N, 
Ci ^ bi - bi-i - ai+i + cii, 1 ^i < N + 1 {or i e E). (9.26) 

Proposition 9.11. Suppose that the given rates {ai,bi,Ci : i £ E) satisfy (9.25). 
Define Xo{ai,bi,Ci) as in Proposition 9.1 without preassuming that Ci ^ for all 
i E E. Next, define (dj,6j) and (di,6j) by (9.24). 

(1) If Y^f=2l^i^T-i ~ then Xo{ai,bi,Ci) ^ Aq, where Aq is defined by (4-1) 

with rates (d^, bi) . 

(2) Otherwise, Ao(aj,6j,Cj) ^ Ai, where Ai is defined by (6.1) with rates 
{ai,bi). 

(3) The equality sign of the conclusions in parts (1) and (2) holds provided it 
does in (9.25). 

Proof, (a) As an application of Proposition 9.1, without loss of generality, we may 
and will assume that the equality sign for q in (9.26) holds. Then, we prove that 
the equality sign of the conclusions in parts (1) and (2) holds. 
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Clearly, we have 

fii = 1, /iidi = /3, fii = fiiUi = 2 < N + 1. (9.27) 

Mi 

(b) Recall the operators: 

^/W = bi{fi+i - fi) + ai{fi-i - fi) - Cifi, 

^f(i) = Hfi+i - fi) + ai(/i_i - fi), / G JT, /o = 0, fN+i = if AT < oo. 

Clearly, Ao(ai,6i,Cj) is the principal eigenvalue of Q, and the idea is describing it 
in terms of the first eigenvalue Amin of Vt. Let U be the diagonal matrix with 
diagonal elements {^i : i & E). Then is simply the diagonal matrix with 
diagonal elements {jli : i E E). For each function h with = and /ijv+i = if 
AT < oo, by (9.27), (9.24) and (9.26), we have 

(nU~^h){i) = bi{fii+ihi+i - fiihi) + ai[jli_ihi_i - jlihi) - Cifiihi 

= ~ Ai^j) + bi-i[{j,i-\hi-i — jlihi) 

- [bi - bi-i - aj+i + ai)jlihi 

= {cLi+ifli+ihi+i - bifiihi) + [bi-ifti-ihi-i - Uifiihi) 

= fiibi{hi+i - hi) + fiiai{hi-i - hi) 

= fii^ h{i) 

= {U-^hh){i), 2f^i<N. 

It is easy to check that the identity holds also for i = 1 and i = N, and then for 
all i & E. Multiplying U from the left on the both sides, we obtain 

UnU-'^ = n. (9.28) 

Furthermore, we get 

{f,ng)^ = {u-'f, iunu-')Ug)^ = {Uf,{Um-^)Ug)j, = {f,ng)f, 

for all f,g € J(f^, where the mapping /—)■/:=[// is an isometry from LF'{^) to 
L'^{fi)- Since / G Jf iff / G it follows that the operators and O with the 
same core J(f are isospectral. In particular, Ao(aj,6j,Cj) = Amin- 

(c) For assertion (1), since {fiibi) ^ = oo by assumption, it follows that 
N = oo and the Dirichlet form corresponding to is regular by Proposition 1.3. 
Hence, the minimal and the maximal domains of the Dirichlet form are coincided. 
Therefore, Amin is equal to Aq^'"*^"* replacing the original rates {ai,bi) by (di,6i). 

For assertion (2), since oi > and bN > 0, we come to the setup of Section 7: 
Amin = Aq'''''^^ with (a,, bi) replaced by (dj, hi) . Next, because of {fiibi) ^ < oo, 
by Theorem 7.1, it turns out Amin = Ai in terms of the dual rates (di,^i) of 
{ai,bi). □ 

We now summarize our main qualitative result about Aq. The three parts 
given below are obtained by Corollary 9.9, Proposition 9.1 plus Proposition 1.3, 
and Corollary 9.5, respectively. 
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Summary 9.12. We have Aq > whenever N < oo. Next, let N = oo. Then 

(1) Ao > if lim ^_^^ c„ > (in particular, j/infj Ci > 0). 

(2) Ao = Xo{ai,bi,Ci) > Ao(aj, 6j, Cjl{i}) > which can be checked case 
by case by 

(i) Theorem 3.1 when Ci = and 



5:J-<oo; (9.29) 



(ii) Theorems 7.1 and 6.2 when Ci > and (9.29) holds; 
(Hi) Theorem 4-2 when ci > but (9.29) fails. 
(3) Xo = if 

lim ^ ixid / '^iJ,k = 0, lim Umbm ( X] Mi ) =0 ^ /Xj = oo 



m-i-oo 



fc=l ^ i=l 



Open problem 9.13 (Explicit criterion for Aq > 0). As will be seen soon in 
Example 9.18 below, for Aq > 0, it can happen that lim ^_^^ c„ = 0. Hence, the 
simple condition "lim^^^ c„ > 0" in part (1) is sufficient only but not necessary. 
Naturally, this condition becomes necessary for the first one in part (3) for which 
a sufficient condition is lim„_).oo c„ = 0. Thus, it is more or less satisfactory 
whenever (c„) has a limit. Otherwise, there is a gap. In contrast with the first 
condition in part (3), condition (9.23) is sufficient for > but there is still a 
distance to deduce the positivity of Aq in view of Corollary 9.9. 

Next, since we are dealing with the minimal Dirichlet form, a general criterion 
for Hardy-type inequalities (cf. [12; Theorems 7.1 and 7.2]) which was successfully 
used in Section 8, is also available in the present situation, hence, there is already 
a criterion for Aq > in terms of capacity which is unfortunately not explicit. 
More seriously, the technique to produce an explicit result used in [12; pages 134- 
136] does not work at the beginning (replacing a finite number of disjointed finite 
intervals {Ki} by the connected one [minUii^^i, maxUii^i]) in the present setup. 
Thus, it is still an unsolved problem to figure out an explicit criterion for Aq > 
in the present setup. 

It is our position to illustrate by examples the application of the results ob- 
tained in this section. First, by using Proposition 9.11 and (9.26), it is easy to 
transfer the examples given in Sections 3 and 6 to the present context. However, 
most of the resulting killing rates arc rather simple. We are now going to con- 
struct some new examples, all of them are out of the scope of Proposition 9.11. 
In the most cases, we use simple (aj,6j) and pay more attention on (cj). Let us 
begin with the following simplest case. 

Example 9.14. Let 

-bi - ci ^1 ^ 

02 -02 - C2 j ' 
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Then as in Examples 7.5 (2), we have 

Ao = ^ {0.2 + 61 + ci + C2 - \/ (02 + C2 - 61 - ci)2 + 40261 ^ , 

with eigenvector 



202 



02 



+ C2 - 61 - Ci + \/ (a2 + C2 - 61 - Ci)2 + 402^1 



Even in such a simple case, the role for Aq played by the parameters a^, hi, and 
Cj is ambiguous. For instance, since ci — £2 = ci — C2 (cfc := — ci A C2), one 
can separate out the constant Ci A C2 from the above expression of Ao- However, 
this obvious separation property becomes completely mazed for the next example 
having three states only. 



Example 9.15. Let 



Q 



Then 



-hi — ci hi 

02 — a2 — 62 — 02 62 
as -as - cs , 



1 o -U 
Ao = ~3^i + ~3~ 



1 / V 

- arc cos 

3 V 2 



-U 



-3/2 



27r 



with eigenvector 



^ 61(03 + C3 - Aq) ^ Cs - Aq ^~ 

03(61 + ci - Ao) ' as ' " 



where 



(9.30) 



= 72 - 7iV3, = 73 - 7172/3 + 2(7i/3)3, 

and + 7iA^ + 72A + 7s the eigenpolynomial of —Q with coefficients: 

71 = -0.2 - as - 61 - 62 - ci - C2 - Cs, 

72 = a2as + 61 as + cias + C2as + 6162 + a2Ci + 62C1 + 61 C2 + C1C2 + a2Cs 

+ 61 Cs + 62CS + CiCs + C2Cs, 

73 = -a203Cl - 0361 C2 - a3CiC2 - 6162C3 - a2CiCs - 62C1CS - 61C2CS - C1C2C3. 



Proof. Since the eigenvalues of — Q are all real, it is easier to write them down. 
By using the notation given above, the eigenvalues of —Q can be expressed as 



1 n -U 

--71 + 2\ cos 

3 ' V 3 



1 / V 

-arccos( -- 



-U 



-3/2 



+ 



2k'K 



A; = 0,1, 2. 
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Among them, the minimal one is Aq given in (9.30). Clearly, the solution is indeed 
rather complicated in view of the coefficients of the eigenpolynomial. □ 

To see the role played by the killing rate (q), in the following examples, we 
restrict ourselves to the case that Ao(ai,^i,0) = (and then N = oo). The 
examples are arranged according the increasing order of the polynomial rates (cj) 
and (bi). Actually, all the examples in the paper are either standard or constructed 
by using simple rates and simple eigenfunctions. They are used first as a guidance 
of the study and then to justify the power of the theoretic results. 

In contract to the explosive case (cf. Theorem 3.1), Aq can still be zero for the 
process having positive killing rate, as shown by the following example. 

Example 9.16. Let h\ = 1, ai = hi = 1 for i ^ 2, and let (ci) satisfy lim„_>.oo c„ = 
0. Then we have Aq = 0, even though Ci can be very large locally. 

Proof. Apply Corollary 9.5. □ 

Example 9.17. Let at = h = I for i ^ 2 and Ci = p-^^S - 1)^ (^ > 0) for i ^ 1. 
Then for every ai ^ and bi > 0, we have Aq = — 1)^. 

Proof. Since Xo(ai,bi,0) = and (cj) is a constant, this is a consequence of part 
(3) of Proposition 9.1. 

Note that the lower estimates given by Proposition 9.1, Theorem 9.3, and (9.4) 
are all sharp for this example. To see this, simply choose Vi = 1 in (9.2) and 
(9.4). We now consider a more specific situation: ai = 0, bi = 1 — f3 and 
(3 <E (0,1). If we set Vi = /3~^, then it is easy to check that Ri{v) = and 
so inft,>o supj>]^ i?i(t;) = 0. This shows that the truncating procedure used in 
Theorem 9.3 for the upper estimate is necessary in the case that the function / 
defined by (9.5) does not belong to L'^{fi), even though Ri{v) is a constant. In 
the present case, Vi > 1 for all i and so the corresponding function / is strictly 
increasing. Since /li is a constant for i ^ 2, it is clear that /Xj = oo and then 



Example 9.18. Let oi = 0, 6i = 5/2, = 2 and 6, = 1 for i ^ 2, Cj = for odd 
i and Ci = 13/6 for even i. Then Aq = 5/6. The upper bound provided by (9.9) is 
approximately 1.03. For the lower estimate, Proposition 9.11 is available but not 
Corollary 9.9. 

Proof, (a) Let Vi = \ + (— 1)Y3. Then it is easy to check that Ri{v) = 5/6. Next, 
define 



We claim that g G L'^{fi) by using Kummer's test. To do so, note that to study the 
convergence/divergence of the series Ylinl^n9ni the constant k defined by (3.13) 
takes a simpler form as follows: 



fiL\^,). □ 



n-l 



k=l 



(9.31) 




(9.32) 
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Now, because g G L'^{lj) and 

-^g/g = R{v) = 5/6, 

we have Aq = 5/6 by Theorem 9.3. Clearly, this eigenfunction g of Aq is not 

monotone since gij^i/gi = Vi =2/3 for odd i and = 4/3 for even i. 

(b) Next, we study the upper estimates of Aq. First, we have 

5 5 
/lii = 1, //i = — , i ^ 2; i^ibi = — , i > 1. 

The upper bound provided by (9.9) is approximately 1.03. 

(c) For a lower estimate, we apply Proposition 9.11 (2). The modified rates 
are as follows: = 2 (i ^ 1), bi = 5/2, and 6^ = 1 (i ^ 2). However, (c,) does 
not satisfy (9.25) at i = 2. We now replace (q) by (q := q + 1/6) and choose 
bo = /3 = 2/3. Then (cj) satisfy (9.25). With the modified (q), we are in the 
ergodic case, and moreover, Ai = (\/2 — l) with eigenfunction g: % = —1 and 



a. - J_2*/2 
~ 20 



101 + 6OV2 + (^41 - 25V2) i 



i ^ 1. 



Therefore, by Proposition 9.11 (2), we obtain Ao(ai, bi,Ci) ^ (\/2 — l)^. Returning 
to the original (cj) by Proposition 9.1 (2), we get a rough lower bound as follows: 

Ao = Xo{ai,bi,Ci) - ^ ^ ^ - 2V2 «i 0.005. 

6 D 

Before moving further, let us remark that if only 60 is changed from 2/3 to 

1/2, then for the (a^, 6^) -process, we still have Ai = (\/2 — l) with a similar 
eigenfunction g: go = —1 and 



— 2*/2 
10 



67 + 42V2 + (27 - 17V2) 



i > 1. 



Now, as an application of Proposition 9.11 (2) with the original (c^) replacing 
C2 = 4/3 by C2 = 3/2 only, the resulting Aq has a lower bound (a/2 — l) « 0.17. 

(d) To apply Corollary 9.9, we write q = 13(1 + (-1)*)/12 and use (9.20) and 
(9.21): 

^^ = inf i^£±£^, = inf -, r]G[0,C,], (9.33) 

Zi + eyi ^ odd i^3l + exi-ri{zi+eyi)' 

Note that the numerator of given in (9.33) is independent of i but in the 
denominator, Xj, yi and Zi all tend to infinity as i — )• 00. To avoid the trivial 
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estimate, one needs to cancel the leading term in i of exi — r]{zi + eyi) in the 
denominator. This leads to the following solution: 

65 e 



V 



9(2 + 5 e) 



Inserting this into exi — r]{zi + eyi), it follows that the denominator of in (9.33) 
becomes 

65 (-2 i + {-ly + 3) £2 + 2 (78 i + 13(-1)* - 245) e - 144 
^ 72(5£ + 2) ■ 

Now, in order to remove the leading term in i, the only solution is 

e = 78/65 > 1, 

which does not belong to the domain of e G (0, 1). Therefore, the test function 
used in Corollary 9.9 does not provide enough freedom to cover this example. 

Note that without the killing rate, the process with rates (oj) and (6j) is expo- 
nentially ergodic and so Ao(aj,6j,0) =0. □ 

For the following examples, we assume that = 6j for i ^ 2. Then 

fj,i = 1, fii = 6ia,"\ i ^ 2; /Xjbj = h, 1. 

The quantities ^£ and defined in (9.20) and (9.21), respectively, are now deter- 
mined by 

i-j ~ sr^ i-j ^ - 1 



Example 9.19. Let ai =0,bi = 2/3(1 -/3)(1 -2/3)-i (/3 G (0,1/2)), a, = b,=^i 
for i ^ 2, Cj = (1 — /3)^(i — 1) for i ^ 1. Then Xq = 2^(1 — P). In the special case 
that P = we have Ao = 3/8. The upper and lower bounds provided by (9.9) 
and Corollary 9.9 are 3/4 and approximately 0.274, respectively. 

Proof. Let Vi = /3(1 + z"^) for i > 1. Then Ri{v) = 2^(1 - 13). By Rummer's 
test (cf. (9.32)), the corresponding function g defined by (9.31) belongs to L'^{fi). 
Hence, the assertion follows from Theorem 9.3. Note that Vi < 1 for all i, and g is 
strictly decreasing even though ci = < Aq and q > Aq for alH > (1+/3)(1 — /3)~^ 
(compare with (2.5)). 

As an application of (9.9) with (£, m) = (1, 1) or (3,4), we obtain 

23 - 40/3 + 23/32 



Ao ^ (1 - P) 



\ 1 - 2/3 / \ 



2(8-11/3) 

To study the lower bound, for simplicity, wc let /3 = 1/4. Then Aq = 3/8 and the 
upper bound in the last formula is 3/4. Choose e = (\/409 — 5)/24 so that the 
infimum = (29 — \/409)/32 0.274 is attained simultaneously at z = 2 and 
i = 3. Since < C2, the set {i ^ 2 : q < S,^} is empty. Therefore, the lower 
bound provided by Corollary 9.9 is approximately 0.274. □ 

For the following two examples, without the killing rate, the process is expo- 
nentially ergodic and so Ao(aj,6j,0) = 0. 
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Example 9.20. Let ai =0, bi = 4/5, Qi = bi = for 2, and 

"8 2 



+ 5 



i > 1. 



3i-8 3i-4 

Then Aq = 4. The upper and lower bounds provided by (9.9) and Corollary 9.9 
are 14/3 and approximately 2.82, respectively. 

Proof. The proof is similar as before using 

Vi = l- „."*". , i ^ 1. 
3z — 4 

Note that Cj has minimum at i = 2. The upper bound provided by (9.9) with 
{£,m) = (2,2) is 14/3. The lower bound produced by Corollary 9.9 with e = 1 
is 48/17 ~ 2.82. Since ci > 0, the parameter e = 1 is allowed. Then = 48/11 
is attained at i = 3, and Ce = 48/17 is attained at i = 2 (noting that the set 
{i ^ 2 : Ci < is a singleton {2}). □ 

Example 9.21. Let ai = 0, bi = 3/2, ci = 15, 

ai = bi=i{i- 4-^)(12i2 - 31i + 27), i ^ 2, 
4 1 o 301 227 

Then Aq = 119/8 = 14.875. The upper and lower bounds provided by (9.9) and 
Corollary 9.9 are approximately 15.42 and 13.18, respectively. 

Proof. Note that Cj is convex and has its minimum at i = 2. For 

3 2 7 

Vi = -- - + — — -, I ^ 1, 

4 I 4z — 1 

we have Ri(y) = 119/8. Note that vi> 1 and Vi <1 for all i ^ 2. The function g 
defined by (9.31) is not monotone but is bounded. Next, since jii ~ we have 
g G L^(/u). The assertion now follows from Theorem 9.3. 

Clearly, infill Ci = 11/4. The upper bound provided by (9.9) with {i,m) = 
(2,4) is approximately 15.42. To get a lower estimate, we apply Corollary 9.9. 
Because ci > 0, we can choose e = 1. Then = 354679/29504 is attained at 
i = 4. Next, since the set {i ^ 2 : q < ^e} is a singleton {2}, we need only to 
compute at i = 2: Q « 10.43. Thus, the lower bound produced by Corollary 9.9 
is approximately 13.18. □ 

To conclude this section, we return to the uniqueness problem for birth-death 
processes with killing of the Dirichlet form as discussed at the end of Section 1. 
Certainly, the problem is meaningful only if = oo. Recall that for a given 
Q-matrix, not necessarily conservative (i.e., may have killing), the exit space 
is the set of the solutions (uj) to the following equation: 

(A7 - Q)u = 0, 
^ ' A>0. 

Note that the dimension of is independent of A > 0. By (2.5) replacing A 
with —A, it follows that the non-trivial exit solution, if it exists, is unique and is 
strictly increasing. 
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Theorem 9.22 (Uniqueness of the Dirichlet form). Let N = oo. 

(1) The Dirichlet form satisfying the Kolmogorov's equations is unique if^\ = 
{0}. Equivalently, 

cx) ^ n 

^— -^/Xfc(H-Cfc) = oo. (9.34) 

n=l fc=l 

(2) Let Y^i^El^i < Then the Dirichlet form is unique iff 

/XjCj < oo and — — = oo. 

(3) Let 5^jgEA*i = oo. Then the Dirichlet form is unique if 

either inf ^ P^'"^{\) > or ^ mci < oo 
ieE jeE ieE 

holds, where E = {i & E : Ci > 0}. 

Here are some comments about the theorem. 

(i) Suppose that only a finite number of q are non-zero. 

Then condition (9.34) is equivalent to (1.2) [Certainly in this item, we 
are using the modified (1.2) and (1.3) by removing from the state space]: 

5Z — jr ^i^k^^ — 7- XI '"'^(^ + cfc) ^ c X — — X nk, 

n=l ^""^ fc=l n=l fc=l n=l ^""^ fc=l 

where C = maXi:c->o(l + q) < oo. Hence, condition (9.34) is stronger 
than (1.3). In this case, condition (9.34) is even not needed in part (3) 
where the last two conditions are automatic. 

When J2i l^i < (1-3) is equivalent to (1.2) which coincides with 
(9.34). When J2il^i — both (1.3) and part (3) hold. Therefore, if only 
a finite number of Ci are non-zero, then we have 

criterion (1.3) <^=^ one of parts (2) and (3) holds. 

(ii) When 7^ for infinite number of z, except condition (9.34), an additional 

condition on the killing rates (q) is required. The condition means that 
if Cj increases very fast, then there exist some Dirichlet forms that do not 
satisfy the Kolmogorov equations. 

(iii) The second condition in part (3) is the same as the one in part (2). For 
the first condition in part (3), it is easy to write down some more explicit 
sufficient conditions. This is due to the following fact. Since for each fixed 
J, {Pj™'"(A) : z G is the minimal solution to the equations 

\ ^ Qik , ^ij . ^ „ 

Xi = 2_^ T-^—xk + T-— , « e E, 

f-f. X + qi X + qi 
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by the linear combination theorem, { J2jeE -^ij^W '■ i & is the mini- 
mal solution to the equations 



EQik ,1 -^77 



This minimal solution (a:;|) can be obtained in the following way. Let 



x^^^ = ^ i G E 



...X + Qi^ A + , 



Then x\"^ t < as n ^ oo for every i G -B (cf. [10; §2.1 and §2.2]). 
Hence, for each n ^ 1, xf^^ is a lower bound of J2jeE ^ij^'^W- 
this discussion, it is clear that the first condition in part (3) is also a 
restriction on the growing of the killing rates (q). This is consistent with 
the second condition there. It is regretted that we do not know at the 
moment whether the conditions in part (3) are necessary or not. 

Proof of Theorem 9.22. Part (1) follows from [10; Theorem 3.2] and Chen et. 
al. (2005) [1] with a fictitious state 0. The last cited result is an application of 
the single birth processes. Noting that if J2i^E f^i ^ °° ^i^E l^i'^i ^ 
then (9.34) holds iff YliieEil^i^i)~^ ~ hence, part (2) is a special case of [10; 
Theorem 6.42]. Next, noting that the unique exit solution is strictly increasing, 
when = oo, we have H L^i/j) = {0}. Hence, part (3) is a particular 

application of [10; Theorem 6.41]. □ 



10. Notes 

10.1 Open problems and basic estimates for diffusions. 

Having seen such a long paper, the reader may feel strange if we claim that 
the story is still incomplete even in the context of birth-death processes. Unfor- 
tunately, it is the case. 

All of the examples we have done so far show that the following facts hold. 

(1) The ratio of the improved upper and lower bounds belongs to [1,2]. 

(2) The sequence is increasing in n and (5„ ^ 5^ for all n. 

(3) The sequences {6^}, {^n}^ and {Sn} all converge to Aq ^ as n ^ oo. 

(4) The relation {fji,r]i) C {k,4k) discussed in Section 6 holds. 

However, there is still no analytic proof for them. The difficulty for the first 

question is that the maximum/minimum of 5i and 5i may locate in different 
places. In the case that (2) would be true, then the story could be simplified since 
we need the first sequence only. For Questions (2) and (3), the assertions are 
numerically justified for almost all of the examples in the paper but the results 
are not included. We have not worked on Question (3) hardly enough since one 
can go ahead only in a finite number of steps in the symbol computation but 



SPEED OF STABILITY FOR BIRTH-DEATH PROCESSES 



123 



the question is certainly meaningful and in the numerical computation, only in 

a few steps one achieves the eigenvalue. For the sequences {fjn} and {rjn}, we 
have similar questions as (1) and (3) about, but the corresponding question (2) 
is answered by Lemma 6.5. 

There is a parallel story for the one-dimensional diffusions. In many cases, one 
can easily guess what the result should be, even though there may exist a new 
difficulty in its proof. For instance, as a combination of the proofs of Theorem 8.2 
and [12; Corollary 7.6], one may prove the following result. 

Theorem 10.1. Consider the minimal diffusion on (—M,N) (M,N ^ oo) with 
operator 

d^ d / b(x) 

L = a(x)——:r + b(x)— ( a(x) > 0, — — is locally integrable 
dx-^ ax \ a[x) 

and Dirichlet boundaries at —M if M < oo, and at N if N < oo. Let C{x) = 
J0 b/a for some fixed reference point 6 G {—M,N) and assume additionally that 
ja is also locally integrable. Denote by A-g, the optimal constant in Poincare- 
type inequality (8.1) with Dirichlet form 

/N 
/'V, fe^^{-M,N), 
-M 

Then satisfies ^ ^ 45b ; where 

^(x,y)\\n^- (10-1) 



Bts 



inf 

-M<x<y<N 



-C 



M 



+ 



iV 



-C 



- 1 



By the way, we prove a dual result of Theorem 10.1 for ergodic diffusions. As 
discussed in the proof of Theorem 7.5, the exponentially ergodic rate often coin- 
cides with the first non-trivial eigenvalue Ai defined below. Consider a diffusion 
process with operator L as in Theorem 10.1, with state space {—M,N) {M,N ^ 
00) and reflecting boundaries at —M if M < 00, and at iV if iV < 00. For 
convenience, we define two measures as follows: 

Scale measure: i'{dx) = e~'-^^^^dx, C{x) := -, 

Je 0, 

where 6 G (— M, N) is a fixed reference point. 
Speed measure: ii{dx) = — dx. 



With these measures, the operator L takes a compact form: 
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Next, suppose that iJ,{—M, N) < oo, and denote by tt the normalized probability 
measure of /x. Set 

= {/ : / is absolutely continuous in (— M, iV)}, 

and define 

Ai = mf{D{f) : f G L'{fi) n ^, 7r(/) = 0, ||/|| = 1}, 

where 

/N 
af'^d^i, / G ^. 
-M 

Clearly, in the definition of Ai, only those / in the set {/ G L'^ijj) H : D{f) < 
oo} are useful. In other words, we are here using the maximal Dirichlet form, as 
in Section 6. 

Theorem 10.2. Let a > 0, a and b be continuous on [—M,N] (or {—M,N] if 
M = oo, for instance). Assume that fi{—M,N) < oo. Then for Xi, we have the 
basic estimate: /A ^ Ai ^ , where 



K ^ = inf 



-M<x<y<N 



1 / rN \ -In 



dfi] + / d/X 
M J \Jy 



](/:-)"• 



(10.3) 



Proof, (a) First we show that for the basic estimate, it suffices to consider the 
finite M and N with smooth a and b. Since a and b are continuous, if M = = oo 
for instance, we may choose Mp, Np ^ oo as p ^ oo such that 9 G {—Mp,Np) 
for all p. Then, by Chen and Wang (1997, Lemma 5.1), we have X^^'^^'^p^ | Ai 
as p — > oo (This is parallel to the localizing procedure used in Section 6). At the 
same time, the isoperimetric constants k^^p^ ^ as p — t- oo (cf. proof of 
Corollary 7.9). Hence, in what follows, we may assume that M, < oo. Next, by 
using the continuity of a and b again, and using a standard smoothing procedure, 
we can choose smooth Op and bp such that — > a and bp b (as p ^ oo) 
uniformly on finite intervals, and furthermore, we can assume that > on 
each fixed closed finite interval. Clearly, the corresponding Kp converges to k as 
p ^ oo. Therefore, without loss of generality, we assume, unless otherwise stated, 
that not only M, N < oo but also a and b are smooth with a > on [—M,N]. 
(b) Recall the following differential form of variational formula for Ai: 



Ai = sup inf 

/g^ xe{-M,N) 



^, af" + {a> + b)f' 



(10.4) 



where 

,^ = {f (^^\-M,N]r\^^{-M,N): f{-M) = f{N) = Q, /|(-m,jv) > O}. 



This is an analog of the variational formula for the lower estimate in Theo- 
rem 6.1 (1). In the original study by Chen and Wang (1997, (2.3)), the state 
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space is the half-line, not finite, but this is not essential. It works also for finite 

state spaces. Besides, it was stated as in (10.4) only. For "=", one simply 
chooses f = g' , where g is the eigenfunction of Ai. This gives us the boundary 
condition: f{—M) = f{N) = since g'{—M) = g'{N) = by assumption. Here, 
one requires that g G '^^(— M, A'^) which is satisfied since we arc now in a finite in- 
terval having smooth a and h. Alternatively, instead of the original coupling proof, 
one may use the analytic one which leads to (6.4) for birth-death processes. 

(c) We are now going to handle with a more general situation: M, iV ^ oo, 
a,b £ 'rf^{~M,N) and a > on { — M,N). Let us define a dual operator L of 
L. In view of the Karlin and McGregor's construction, the dual of a birth-death 
process is simply an exchange of the scale and speed measures jl = u and u = ii 
up to a constant (cf. (5.3)). Thus, in view of (10.2), the dual operator L, as was 
introduced by Cox and Rosier (1983), should be given by 

Again, the speed and scale measures fi and u oi L should be expressed as 

a 



dfl = ^dx, du = e ^dx 



in terms of the coefficients a and 6 of L to be determined now. Because ft = v 
and = jj,, we have 

d/t dz> du dfl 
dx dx dx dx 

It follows that a = a. Then using the equation jl = v, we get 

C = -C + \oga = -C + log a. 

Thus, from 

a a a 

we get b = a' — b. Therefore, the dual operator L has the following expression: 

For the dual process, the Dirichlet boundary is endowed at —M and N (cf. proof 
(e) below). Clearly, the dual operator L is symmetric on L^(e~^dx). We remark 
that the assumption on a and b can be weakened in this paragraph, 
(d) Define a Schrodinger operator as follows: 

d^ d 

Ls = a{x)^ + {a'{x) + b{x))— + b'{x) 
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with Dirichlet boundaries at — M and N provided they are finite. Clearly, Ls is 
symmetric on L^(^e'^dx). Denote by the principal eigenvalue of Ls: 

In the setup of (b), formula (10.4) becomes 

Ai = sup inf — ^(x). 
/ejr xe(-M,JV) / 

This leads to the study on A5. 

(e) An elementary computation shows that 

e^Lse-^ = L. (10.8) 

Note that 

/N /-N /-N 

-M J-M J-M 

where the mapping f ^ f '■= e*^/ is an isomctry from L^(e'-'dx) to (e~'-'da:) , 
and that / G "rf^i-M, iV) iff / G '^i{-M, N). Since "tfH-M, N) is also a common 
core of Ls and L by the assumption on the coefficients a and b, it follows that the 
operators Ls and L with the same core '^q°{—M, N) are isospectral. In particular, 
we have A5 = Aq- When M, N < OC', this means that L has Dirichlet boundaries 
at —M and N since so does Ls- Now, the basic estimates for Xs can be obtained 
in terms of the ones for Aq, as will be shown in part (f) below. 
To go back to Ai, noting that by (10.8) again, we also have 

-Lsf _ -{e^Lse-^){e^f) _ -Lf 
f / ■ 

By (a), we can assume that M,N < 00 and a > on [—M,N]. Prom Shiozawa 
and Takeda (2005) and X. Zhang (2007), it is known that 

Ao ^ sup inf — ^(x) 
/e^ xe{-M,N) f 

(i.e., Barta's inequality). To see that the equality sign holds, simply choose / 
to be the eigenfunction g of Aq. The fact that g e ^ is guaranteed by the 
assumptions that M,N < 00, a and b are continuous, and a > on [—M,N]. 
This is a standard (regular) Sturm-Liouville eigenvalue problem. The property 
g\(-M,N) > is due to the fact that Aq is the minimal eigenvalue. We have thus 
returned to Ai from Aq through A5. 
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(f) For the dual operator L defined in part (c), applying Theorem 10.1 to 
B = L^{fi), we obtain k~^/4 ^ Aq ^ k~^, where 

Now, the theorem fohows by the dual transform fi = v and u = ^. 
Finally, the proof of Theorem 10.2 can be summarized as follows: 

Ai for general M, N and continuous a, h 
— >■ Ai for finite M, N and smooth a, h 

(by approximating and smoothing procedure) 
— > \s (by coupling method leading to the Schrodinger operator) 

—7- Aq (by isometry in terms of the dual operator) 

—7- basic estimate of Aq (by capacitary method: Theorem 10.1) 

-> basic estimate of Ai (by duality). □ 

Actually, we have also proved the following result (cf. parts (c)-(f) in the last 
proof) which is parallel to Proposition 9.11. 

Proposition 10.3. Let M, N ^ oo, a,b e '^^(-M,iV) and a > on (-M,iV). 
Then for the Schrodinger operator Lg on L^(e^dx) having the form (10.7) with 

Dirichlet boundaries at —M if M < oo, and at N if N < oo, we have Xs = Aq, 
and furthermore, k~^/4 ^ As ^ k~^, where k is defined by (10.3). 

The following simplified estimate of 

^(10.3) ig 

helpful in practice. Recall that by 
assumption, ii{—M,N) < oo. Let m{f^) be the median of fj, (i.e., iJ,{—M,m{fj,)) = 
IJL{m{ii) , N)) . Given x G {—M,m{ij)), let y = y{x) be the unique solution to the 
equation: fi{y, N) = fi{—M, x). The A-G inequality a + ^ ^ 1\foL^ suggests the 
use of y(a;), which then leads to a simpler bound: 

^(10.3)^2-1 sup ii{-M,x)v{x,y{x)). 

xe{-M,m{n)) 

We remark that the equality sign here holds in some cases, but the inequality sign 
can happen in general. Anyhow, this provides us a guidance in seeking for the 
infimum in (10.3). Certainly, the similar discussion is meaningful for 

Having Theorems 10.1 and 10.2 at hand, the basic estimates in the other cases 
(A^^ and A^'^) mentioned in Section 1 should be clear. 

The study on the one-dimensional case provides a comparison tool for the study 
on the higher dimensional situation, as we did a lot before. Hence, there is no 
doubt for the development in the higher dimensional context. 

10.2 h- Transform. 

In an earlier draft of this paper (roughly speaking, up to Theorem 7.1 plus 
a part of Theorem 9.3), the author mentioned an open question: how to handle 



inf 

-M<x<y<N 



dz) 



M 



+ 



d7> 
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the case that (1.3) fails? Then two answers have appeared. The first one is the 

use of so-called /i-transform by Wang (2008a) where the transient case studied in 
Section 7 is transferred into the one studied in Section 4. Next, with the help 
of the duality given in Theorem 7.1, the ergodic case studied in Section 6 can 
be also transferred into the one studied in Section 4. In this way, with a Tise of 
Theorem 4.2, Wang obtains a criterion for Ai (Section 6) with a factor 4. To have 
a taste of this technique, let us quote a particular result here. 

Theorem 10.4 (Wang (2008a, Theorem 1.2)). Set hi = YljLil^j- Then we have 
d~^/4 ^ Ai ^ 6~^, where 



sup (^-l-)S2^hl (10.9) 

l^i<7V+l \hi ho J HiUi 

J — * 



Comparing this result with Corollary 6.6, the factor 4 is in common but the 
isoperimetric constants are quite different. The advantage here is that only 
one variable is required in the supremum, but in Corollary 6.6 two variables are 
needed. The price one has to pay to (10.9) is involving a new quantity h. The 
natural extension of Corollary 6.6 to the whole line (Corollary 7.9) exhibits an 
interesting symmetry of the left and the right half-lines. Such an extension of 
Theorem 10.4 with the same factor 4 is unclear to the author. Along the same 
line and using [9], Wang then extends the results to Poincare-type inequalities as 
well as functional inequalities, see Wang (2008b, c). Clearly, Wang's papers show 
that the ^-transform is a powerful tool and may be useful in other cases. 

While the author's solution to the above open question is the use of the max- 
imal process as included into this version of the paper. As shown in the paper, 
Corollary 6.6 comes from the author's previous general result without using the 
^-transform. An interesting question in mind is to use the variational formulas 
in Section 6 to derive Corollary 6.6 directly. Besides, a direct generalization of 
Sections 2,3, and 7 to the Poincare-type inequalities is still meaningful in practice 
since the formulas are quite different (in view of Theorem 10.4) and some of them 
may be more practicable. 

10.3 Remark on some knov^n results. 

As mentioned in Section 5, duality (5.1) goes back to Karlin and McGregor 
(1957b). The author learned this technique mainly from van Doorn (1981; 1985) 

based on which the proof of the basic result Ai = a* was done, cf. [2]. It is now 
known that such a result holds in a very general setup as indicated in the proof 
of Theorem 7.4. 

We now discuss the situation that (1.3) holds. Then there are three cases. 

(1) Mi = c)o and J2iit^ibi)~^ < oo. 

(2) J2il^i <oo and = oo. 

(3) T.ifJ'i = J2iiP'ibi)~'^ = oo. 

First, let bo > 0. In cases (1) or (3), by Theorem 2.4 (1) and Proposition 2.7 (1), 
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X^'^^ is equal to 



sup inf [oj+i + bi- ai/vi-i - bi+iVi]. 
y := {v : v-i is free, Vi > for all i ^ 0}. 



(10.10) 



In case (2), by Theorem 6.1 (1), Ai can be expressed by (10.10). Thus, in view 
of Proposition 1.2 and [2; Theorem 5.3], the convergence rate a* can be also 
expressed by (10.10). 

Next, let bo = 0. Then in case (2), by Corollary 5.2, Proposition 2.7(1), and 
using (5.8) in an inverse way, it follows that Aq'*'^^ is equal to 



sup inf 



a,: 1 



1 



Vi-l 



+ biil 



(10.11) 



y := {v : vq = oo, Vi > for all i ^ 1}. 



In case (1), Aq is equal to Aq ' . By Theorem 7.1 (1), in terms of Theo- 
rem 6.1 (1) and using (5.8) in an inverse way, we obtain the same expression 
(10.11) for Aq^'^''. Finally, in the degenerated case (3), we indeed have Aq^'^^ = 
X'o ' = which can be expressed as (10.11) by Theorem 7.1 (2). Hence, by 
Proposition 1.2, the convergence rate a* can also be expressed by (10.11). We 
have thus obtained the following result. 

Theorem 10.5 (van Doorn (2002)). Let (1.3) hold. Then the exponential con- 
vergence rate a* is given by (10.10) or (10.11), respectively, according to bo > 
or bo = 0. 

With a slightly different expression, this result was given in van Doorn (2002) 
by the analysis on the extreme zeros of orthogonal polynomials in Karlin and 
McGregor's representation, and was actually implied in van Doom's earlier papers 
(1985; 1987) as mentioned in the paper just cited or in [3]. In the last paper, this 
result was rediscovered in the study on Ai, using the coupling methods. The lower 
estimate was also obtained by Zeifman (1991) using a different method in the case 
that the rates of the processes are bounded, with a missing of the equality. 

A progress made in the paper is removing Condition (1.3) and even (1.2). In 
particular, the situation having finite state spaces is included. This is meaningful 
not only theoretically but also in practice since the infinite situation can be ap- 
proximated by the finite ones. Besides, when 6o > and < oo, the duality 
given by (5.9) is essentially different from (5.8) (cf. Remark 2.8). Prom the other 

point of view, the dual of this case goes to X^'^^ rather than Ag^'^^ However, 
we then have to use the maximal process in Section 6, as we did in Theorem 7.1, 
rather than the minimal one used in Sections 2 and 3, except using (1.3) (which 
is equivalent to (1.2) if J2il^i < °°)- Prom analytical point of view, the use of 
the maximal process is natural since one looks for the inequality to be held for 
the largest class of functions, as illustrated by the weighted Hardy inequality in 
Section 4.1. 
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In van Doorn (2002) , some variational formulas of difference form for the upper 

bound of a* are also presented but wc do not use them here. As far as we know, 
the criterion for a* > (Theorem 1.5) has been open for quite a long time; it 
was answered in the ergodic case only till [6] in terms of the study on the first 
non-trivial eigenvalue Ai. For which, the criterion was obtained independently by 
Miclo (1999) based on the weighted Hardy's inequality. Criterion 3.1 follows from 
the variational formulas of single summation form (part (2) of Theorem 2.4), but 
it is not obvious at all to deduce the criterion from (10.10) (or dually from (10.11)) 
directly. More clearly, the variational formula of the difference form for the lower 
bound given in (9.3) which is closely related to (10.11) was known for some years 
and works in a more general setup, but an explicit criterion for the killing case 
is still open (Open Problem 9.13). Anyhow, having the duality (Corollary 5.2 
and Theorem 7.1) at hand. Theorem 1.3 is essentially known from [6], except 
the basic estimates in the ergodic case as well as in the setting of Section 7 is 
presented here for the first time. The techniqiie adopted in this paper depends 
heavily on the spectral theory, potential theory, and harmonic analysis. In the 
transient continuous context. Criterion 3.1 was obtained by Maz'ja (1985, §1.3), 
as a straightforward consequence of Muckenhoupt (1972). The discrete version 
was proved by Mao (2002, Proposition A. 2). In these quoted papers, the problem 
in a more general (L', L^')-setup was done. 

In the continuous context, the Hardy- type or Sobolev-type inequalities (cf. 
Theorem 10.1) were studied first by P. Gurka and then by Opic and Kufncr 
(1990, Theorem 8.3). Instead of ^™™(D), they considered the following class 
of functions: the absolutely continuous functions vanishing at — M and N. This 
seems not essential in view of Aq^'^'' = With a different but equivalent 

isoperimetric constant (i.e., replacing the sum in (10.1) by maximum "V"), they 
obtained upper and lower bounds with ratio 2 « 22, where uj = (\/5 + l) /2 is 
the gold section number. By the way, we mention that the use of weight functions 
w and V in (8.6) in the cited book is formally more general than our setup. One 
can first assume that w and v are positive, otherwise replace them by w + l/n 
and V + 1/ra, respectively, and then pass to the limit as n ^ oo. Next, it is easy 
to rewrite w and v as e*-^/a and e'^ for some functions C and a > 0. Note that 
only C and a (without using h) are needed to deduce the basic estimates in our 
proof. Again, in the continuous context, the splitting technique was also used in 
Theorem 8.8 of the book just quoted where some basic estimates were derived in 
terms of an isoperimetric constant, up to a factor 8. Their isoperimetric constant 
is parallel to the right-hand side of (7.13) replacing Ag^ by the corresponding 
j(3-i)± depending on 6 (certainly, without using the parameter 7). Our Example 
8.9 is an analog of Examples 6.13 and 8.16 in the quoted book. In contrast with 
our probabilistic-analytic proof here, their proof is direct, analytic, and works 
in a more general (L"?, LP)-setup. We have not seen the discrete analog of their 
results in the literature. In the (L^, LP)-sense {p ^ 1), the variational formulas 
in the continuous context were obtained in Jin (2006) but it remains open for 
the more general (L'', L^')-sctup. Even though it is a typical Sturm-Liouville 
eigenvalue problem having richer literature, we are unable to find an analog of 
Theorem 10.2. 



SPEED OF STABILITY FOR BIRTH-DEATH PROCESSES 



131 



Finally, in computing the examples in the paper, the author uses the software 
Mathematica. All the examples were checked by Ling-Di Wang and Chi Zhang 
using MatLab. Most of the author's papers cited here can be found in [8]. 
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